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Preface to the Second Edition 


This second edition of the Theory and Practice of the Testing 
of Gas Wells has resulted from a detailed review of the first or 
Discussion Draft carried out during 1964 and 1965. The Board reques- 
ted the view of the Alberta oil and gas industry on the Discussion 
Draft through the Alberta Division of the Canadian Petroleum Associa- 
tion. The Association established a special committee, under the 
chairmanship of Mr. R. L. Taylor of Shell Canada Limited, to review 
the document and the procedures proposed in it. The other committee 
members were G. C. Whittaker, Pacific Petroleums Ltd., V. P. Peters, 
Pan American Petroleum Corporation, J. Hnatiuk, The British American Oil 
Company Limited, H. S. Simpson, Imperial Oil Limited, and S. M. Thorne, 
Socony Mobil Oil of Canada, Ltd. This group spent many hours in a de- 
tailed study of the Discussion Draft and made many useful suggestions 
which have been incorporated in the Second Edition. The Board 
acknowledges this valuable contribution. 

In addition, Dr. D. L. Katz of the University of Michigan, 
Mr. K. Aziz of Rice University and the late Mr. F. K. Beach, Petroleum 
Consultant, were good enough to send constructive comments for incor= 
poration into the Second Edition. These also are gratefully acknow- 
ledged. 

Mr. G. J. DeSorcy, Chief Gas Engineer for the Board, accepted 
the major responsibility for preparing the revisions with Mr. A. S. 
Telford of the Board's Gas Department giving able assistance. 

In addition to introducing the suggestions of the aforementioned, 
Mr. DeSorcy has made a number of revisions of his own and as suggested 


by other members of the Board organization. 
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.The Board has been pleased with the acceptance of the 
Discussion Draft by industry and believes that this Second Edition, 
incorporating the suggestions received from industry, should serve 
the objective originally set - that of improving engineering practice 


in the testing of gas wells and aiding in conservation. 


G. W. Govier 
Chairman 
Oil and Gas Conservation Board 


June, 1965 


Preface to February 1964 Edition 


This manual on the Theory and Practice of the Testing of Gas Wells 
has been prepared to serve a growing need in the Province of Alberta, It is 
in what might be called "discussion draft" form, It is hoped that sugges- 
tions for its improvement will be received from industry during 1964 and it 
is planned that a revised and improved version will be published in 1965. 

Many advances have been made in the understanding both of the 
reservoir and the well bore flow of natural gas since the publication of 
the well known "Monograph 7" of the U.S. Bureau of Mines in 1937. Few of 
these developments, however, have been incorporated into routine well 
testing procedures and the engineer who is not a specialist in the field 
may not be familiar with the technical literature describing them, 

The manuals published by the Texas Railroad Commission, the Kansas 
State Corporation Commission, and recently by the Interstate Oil Compact 
Commission have contributed greatly to the precision of definition of 
certain of the tests and to a standardization of methods of calculation. 
These publications, however, were not intended to include reviews or 
criticism of the pertinent theory and except for the last mentioned they 
deal only with "stabilized flow" testing. 

One of the most significant improvements in the scientific under- 
standing of gas well testing is in connection with the unsteady state or 
unstabilized flow behavior. It seems appropriate to take advantage of this 
and other developments and to design well testing procedures accordingly. 

The preparation of the first draft of this manual was carried out 
in 1958-59 under my guidance as a class project in the graduate course, 
Advanced Natural Gas Engineering, at the University of Alberta, Edmonton, 


by graduate students: D. Batcheller, T. Fekete, I. Nielsen, C. Winter and 


0 Vall Ee 


S, Qayum. Mr. Fo Werth, Instructor in Petroleum Engineering, assisted in 


the preparation of certain of the tables and figures. Valuable assistance 


also was given at this time by Mro Je Ge Stabback, then Chic? Gas) se aaa 


of the Oil and Gas Conservation Board, and by representatives of the 
Alberta gas industry, including The British American Oil Company Limited, 
Northwestern Utilities Limited, Pacific Petroleums Limited, Shell Canada 
Limited and Imperial Oil Limited, Some of these representatives attended 
class discussions and many made data from company files available. 

The first draft was reworked and revised chiefly by Messrs. J. 
Pletcher, Assistant Gas Engineer and Mr. Se A. Qayum, Temporary Assistant 
Gas Engineer of the staff of the Oil and Gas Conservation Board. Mr. K. 
Aziz, Assistant Professor of Petroleum Engineering, University of Alberta, 
also gave assistance at this time, This work was carried out inter- 
mittently as time permitted during 1961 and 1962. Finally, in mid-1963, 
Mr. Go J. DeSorcy, Assistant Chief Gas Engineer of the Board was assigned 
the full time task of preparing a draft suitable for publication and 
distribution to industry. Mr. DeSorcy and I altered the emphasis of the 
early draft from one directed to flow tests and related bottom hole 
pressure calculations to the broader one of the present treatment. Others 
of the staff of the Board have participated and shared in the checking of 
calculations. 

The manual is being released in its present form not because 
it is thought to be a finished product but because the Board believes that 
it will even now be of some value to industry in improving the practice of 
testing gas wells and, therefore, serve the interests of conservation. 
Also, in this way, the Board hopes that it will receive suggestions from 


industry for improvement of the manual, 


The manual suggests certain procedures, those which the Board 


at abal 


believes best, but it is not intended as a directive with respect to 
procedures, Following its revision in 1965, however, certain sections 
of it may be considered suitable for incorporation, by reference, into 


regulations issued under The Oil and Gas Conservation Act, 


G. W. Govier, P. Eng. 
Chairman 
Oil and Gas Conservation Board 


February, 1964 
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NOMEN CLATURE 


6 
1.418. x 10 u 27, ln (r,/r,) 


a = Th (by definition) 


A = cross sectional area of flow, square cm, 


AOF = absolute open flow. 
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b = Th oa Ra pe Fr, (by definition) 
C = back pressure equation coefficient, 
d = inside diameter of flow pipe, inches. 
Dirt = effective flow atanater of annular space, 
feet = D,, ~ Ds. 
Dy = outside diameter of tubing, feet, 
D,, = inside diameter of casing, feet. 
D = inside diameter of flow pipe, feet. 
D, = average grain diameter, feet, 
e = base of naperian logarithm, 
f = friction factor. 
fox = modified friction factor (Cornell and Katz), 
fopp ™ modified friction factor (Fancher, Lewis and Barnes). 
pe a aes (by definition) 
da 
FY = turbulence contribution factor, 
dF = energy losses due to irreversibilities, 
a = dimension conversion factor, 
G = specific gravity of gas, 
G = specific gravity of separator condensate, 


= XV 


specific gravity of gas resulting from combination 


of separator products. 
specific gravity of separator gas. 
pay thickness, feet. 
permeability, millidarcys. 
permeability, darcys. 
vapor - liquid equilibrium ratio for component i. 
length, feet. 
dimensionless flow rate. 
molecular weight of gas, pounds per pound mole. 
reciprocal slope of the back pressure curve. 
Reynolds number. 
modified Reynolds number (Cornell and Katz). 
modified Reynolds number (Fancher, Lewis and Barnes). 
pressure, psia. 
pressure drop, lbs. per sq. ft. 
average pressure, psia. 
critical pressure, psia. 

p 2 
dimensionless pressure = = (by 

£. 

shut-in formation pressure, psia, 
pressure measured during a build-up or drawdown 
test, psia. 
standard pressure = 14.65 psia in Alberta, 
pressure at time of previous survey, psia. 
reduced pressure, 
reduced average pressure. 


reduced sandface pressure. 


definition) 


Peo! dae ae 


reduced wellhead pressure. 

flowing sandface pressure, psia. 
dimensionless drawdown. 

wellhead pressure, psia. 

false value of final static pressure, psia. 


volume flow rate, cu. cm. per sec. 


gas flow rate, millions of cubic feet per day 


at standard conditions. 


radius, feet. 


apparent or steady-state radius of drainage, feet. 


‘ < ; r 
dimensionless radius = me 
8 


exterior boundary radius, feet. 

effective radius of well bore, feet. 

gas law constant. 

gas-condensate ratio, cubic feet per barrel. 


2GL 
53.34 Tog 


time, hours. 

dimensionless time, 

time of flow, hours. 

time to stabilization, hours. 
shut-in time, hours, 
temperature, °R. 

average temperature, °R. 
critical temperature, °R. 
standard temperature = 520°R. 


reduced temperature. 


(by definition) 


(by definition) 


=e xVd = 


reduced average temperature. 

reduced sandface temperature. 

reduced wellhead temperature. 

sandface or formation temperature, °R. 

wellhead temperature, °R. 

variable related to reservoir and flow characteristics. 
specific volume of gas, cubic feet per pound. 

velocity of flow, feet per second. 

apparent velocity, feet per second. 

condensate vaporizing volume ratio, cubic feet per Darrel. 
mole fraction of component i in liquid phase. 

mole fraction of component i in vapor phase. 

a function of u, 

gas compressibility factor. 

average compressibility factor. 

sandface compressibility factor. 

wellhead compressibility factor. 

turbulence factor, Gentimeterss. 

density of gas, lbs-mass per cubic feet. 

gas filled porosity, fraction. 

absolute roughness of pipe. 

relative roughness of pipe. 

viscosity, centipoise (at any pressure and temperature). 
viscosity, centipoise (at one atmosphere pressure and any 
average viscosity, centipoise, 


viscosity of gas, lbs-mass per foot second. 


temperature! 


1. INTRODUCTION 


The ability to analyse the performance and forecast the 
productivity of gas wells with a reasonable degree of accuracy is of 
utmost importance in today's natural gas industry. One of the most 
useful aids in analysing gas well performance is the flowing well 
test. The results of such tests are often used by regulatory bodies 
in setting maximum gas withdrawal rates, They are also employed by 
producing and transportation companies in projecting gas well deliveries. 

A complete analysis and understanding of the results of an 
appropriate well test enables one to determine the rate at which a well 
will flow against a particular pipeline "back" pressure, and also to pre- 
dict the manner in which the flow rate will decline with depletion and 
the resulting drop in reservoir pressure, This type of projection is 
necessary in the preparation of field development programs, the design 
of gathering and pipeline facilities as well as processing plants and 
is often used in the negotiation of gas sales contracts. 

Other important applications of the well test and information 
gathered during testing are in the estimation of gas reserves associated 
with a well or group of wells and in making various types of special 


reservoir studies. 


History 


The maximum theoretical delivery rate of a gas well has been 
described as the absolute "open flow" potential, which in a limited way 
is still a useful concept. The absolute open flow potential is defined 


as the ability of a well to produce against zero sandface back pressure, 


eee 


In the early days of the industry the maximum potential was directly 
measured by gauging the flow of gas from a well which had been opened 
to the atmosphere. Opening a well to atmosphere, however, results in 
a flow rate which is not precisely equivalent to the absolute open flow 
potential, i. e. the flow that would take place against a zero sandface 
back pressure. This maximum flow rate to the atmosphere at the surface 
is often termed the practical open flow potential and for shallow wells 
is very close to the absolute open flow potential. On the other hand, 
for deep and high capacity wells (where the frictional resistance to 
the flow of the gas is considerable) the difference between these flow 
potentials is often very significant. 

The field gauging of open flow potentials often did not give 
reliable results. This was because the flow rates being measured were 
not the true absolute potentials, and also because in the case of wells 
of high capacity the accurate measurement of the flow rates was impossible. 
Even in cases where the measured potential was relatively accurate, the 
reniite were not particularly useful in that they provided no information 
regarding the ability of a well to produce against various pipeline back 
pressures, nor any information regarding the delivery capacity of the well 
at declining reservoir pressures. In addition to this, the practice of 
flowing a well wide open often resulted in damage to the well by the coning 
of water into the well bore or by the "sloughing off" and production of 
sand particles from the formation. Often too, serious wastage of gage re 
sulted. All things considered, the practice of testing a well "wide open" 
has been recognized as undesirable, and methods have been developed for 
assessing productive capacity by conducting well tests at reasonable and 


controlled rates of flow. 


The Back prosalire Test 

The basic work towards arriving at 4 practical test wus 
carried out by the U. S. Bureau of Mines during the period 1929 to 1935. 
This work was documented in two Reports of Investigation by Pierce and 
Rawlins (66)(67) and culminated with the publication of the well-known 
and widely used Monograph 7 (72) of Rawlins and Schellhardt. The early 
work of Pierce and Rawlins actually first proposed the back pressure method 
of testing gas wells. Monograph 7 repeated the same basic material and 
equations and in addition included sections on testing equipment, g48 
measurement and tables which compiled the pressure differences from sand- 
face to wellhead for various static and flowing conditions, These reports 
presented the empirical justification of the back pressure test, described 
the field conduct of the test, and analysed the results of such a test 
based on data obtained on 582 wells throughout the United States. Briefly, 
{t was shown that if a gas well were appropriately tested by flowing it at 
several rates, @ plot of the difference between the square of the static 


reservoir pressure and the square of the flowing bottom hole pressure versus 


the corresponding rate of flow should yield a straight line on log-log 
paper. It was then demonstrated how the plot could be employed to determine 
the well capacity at any flowing sandface pressure including absolute zero 
(corresponding to absolute open flow conditions), and how it could be used 
to predict the behavior of a well as the reservoir pressure declines. 

Figure l-l presents 4 plot of the relationship proposed by 
Pierce and Rawlins for a typical Alberta gas well. The plot shows how the 
difference in pressures squared versus the corresponding flow rates for 
four test conditions result in a straight line relationship. This relation- 


ship was expressed by the following equation 


-1 
2)n (i=d) 


2 
Q ##C (P,. - P. 


where 


Q Fond 

Pe s 

p = 
8 

( = 

n = 


| Reflects A Zero Gandface Pressure 





rate of flow, millions of cubic feet 

per day at standard conditions. 

shut-in formation pressure, psia. 

flowing sandface pressure, psia. 

a coefficient which is determined mainly 
by the reservoir characteristics at the 
particular well. 

an exponent which depends primarily upon 
the type of flow in the reservoir at the 
particular well. The value of 'n' is the 


reciprocal slope of the back pressure Curve. 
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| — CONVENTIONAL BACK PRESSURE TEST PLOT. 


Figure l-1 shgws how the plot may be used to determine the ab- 
solute open flow potentiul ond the capability of a well to flow against 
@ particular pipeline pressure. 

The back pressure test proposed by the U. S. Bureau of Mines 
was widely accepted and used by industry and regulatory authorities for 
@ period of some 20 years. It was normal practice to produce wells be= 
low a maximum gas well allowable based on some constant percentage 
(usually 20 to 35 per cent) of the open flow potential. In 1954, the 
Alberta Oil and Gas Conservation Board recognized that the use of a 
fixed percentage of the open flow over the life of the pool is equiva- 
lent to a steady reduction in the sandface drawdown pressure over the 
entire pool life (97). Aware that an increasing safety factor over 
the pool life was unnecessary, the Board altered its policy and used 
the results of back pressure tests as a guide in setting allowables 
more nearly to reflect a uniform sandface drawdown at any reservoir 
pressure - i. e. allowables equivalent to an increasing percentage 
of the open flow potential as the reservoir pressure declines. 

Pierce and Rawlins recognized that it is not always possible 
to measure the pressures directly down hole and presented a theory and 
method for converting wellhead pressure mcasurements to bottom hole 
or sandface conditions by taking into account the weight of the column 
of gas, and in the case of the flowing well, the pressure drop duc to 
the friction. Their approsch was based on observations at wells which 
were producing from relatively low pressure reservoirs with little or 
no liquid production. With the discovery of deep high-pressure and 
condensate type reservoirs, problems were encountered in determining 
sandface pressures from top hole measurements. Vitter (93) presented 


a method for calculating the sundface pressure for condensate wells, 


ea 


The method combines the gas and liquid products of a separator in deter- 
mining the weight of the gas column for the case where the condensate 
and the gas recovered from the separator exist ina single phase in the 
well bere. Subsequent to the publication of Vitter's paper, Shell O11 
Company (105) pointed out that a slip had been made in the original 
paper with respect to the flowing case, ana further modified the cal- 
culation. The "Modified Vitter's” method of calculating sandface pres~ 
sures has received widespread acceptance and is in common use today. 

More recently, a number of other investigatore (24)(25)(68)(83) 
have proposed improved methods for calculating sandface pressures (particu- 
larly in the flowing case) in the single phase gas well. The best of the 
methods take advantage of recent developments (16)(29)(77)(78)(79) in the 
field of fluid mechanics in determining the friction resistance factor. 

In addition to this, much work has been done on the problem of determining 
the sandface pressure when both a gas and a liquid phase exist in the well 
pore (6)(11)(32)(33)(34)(37)(38)(54) (56) (69)(75) (87). 

An improved understanding of the flow of fluids through porous 
media has made it clear that the empirical straight line relationship pro- 
posed by Pierce and Rawlins only results under conditions of stabilized 
and laminar flow in the reservoir, and in this case the reciprocal slope 
of the back pressure test curve 'n' should cqual unity. 

Work of Fancher, Lewis and Barnes (28) defined turbulent flow 
in porous media at an early date. It has since been shown (7)(27) that 
under conditions of high flow rates in the reservoir, where turbulence 
occurs, the back pressure "line" even for stabilized flow would curve 
upwards or reflect a value of 'n' decreasing from the value of unity 
expected at low flow rates. In most practical cases the variation in 


t,¢ } 
n' is not too great and it may be represented by a constant average 


SN 


value at something less than unity. These considerations point out 
the danger of excessive extrapolation of the back pressure curve. 

An even more important consideration in gas well testing is 
the problem of stabilization. It has been pointed out that the back 
pressure plot as proposed by Pierce and Rawlins will result in a 
straight line only if the flow rates and pressures are essentially 
stabilized. This condition is seldom attained, particularly in low 
permeability reservoirs where it is now recognized that the convention- 
al back pressure test gives meaningless results, unless wells are 


flowed for prohibitive periods of time at each flow rate, 


New Approaches to Testing 
Gas Wells 


Cullender (23) has shown that a reasonably accurate back 
pressure curve may be obtained where the individual flow rates have 
not stabilized, provided, a certain testing procedure is adhered to, 
This procedure involves flowing the well at three or four flow rates 
for periods of equal duration, with each flow period beginning from 
comparable shut-in conditions, The shut-in conditions must be close 
enough to fully stabilized conditions so that any pressure build-up 
still occurring will not affect the recorded pressures during subsequent 
flow tests. (A rule of thumb useful in establishing whether a well 
pressure is sufficiently stabilized is that the pressure should be 
building up at a rate of less than about one-tenth of one per cent of the 
previously recorded pressure in 30 minutes). The resulting data points, 
when plotted in the conventional back pressure test manner, yield a 
straight line and the correct value of 'n'’, An additional point is used 


from a run with an extended flow period approximating stabilized conditions, 
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A line drawn through this point and with the correct :'n', represents the 
true stabilized back pressure test line, The Cullender method, known 
as the "isochronal test’, is consistent with the present day theory of 
gas well testing and is receiving increasing consideration from industry 
and regulatory bodies. 

Another form of isochronal test has been presented by Katz 
et al (49). This "modified isochronal test" is receiving increasing 
support in industry. The essential modification is that each closed= 
in time between flow periods is also of the same duration as the flow 
period, and the new (unstabilized) closed-in pressure is used for cal- 
culating the difference of squares for the next flow point. While this 
procedure is not as readily justified on theoretical grounds as the true 
isochronal test, the results are very nearly the same. 

In addition to these advances in the understanding, conduct 
and interpretation of flow tests, there has been increasing awareness 
of the importance of an accurate determination of the true static 
reservoir pressure. Many investigators (2)(3)(18)(26)(40)(43)(44) (52) 
(58)(59)(61)(62)(65)(89) have discussed the characteristics of pressure 
build-up in reservoirs based on the theory of fluid flow, and several 
methods have been developed for calculating the stan i ived reservoir 
pressure from pressure-time data taken over a period of shut-in time 
substantially less than would be required to attain true stabilization. 
Out of an improved theoretical understanding of pressure build-up and 
drawdown data, has arisen the possibility of using such tests to supple- 


ment or even replace the results of back pressure tests. 


Natural Gas Properties 


In the early days, the only knowledge of gas properties needed 
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in conducting well tests were those required to calculate the flow rates, 
the gas gravity and compressibility factor. Moreover, they were needed 
only at the temperature and pressure of metering. With further develop- 
ment and refinement of the back pressure test, 4 knowledge of the effects 
of temperature and pressure on the compressibility factor has become 
important, and the viscosity of the gas gained significance. In the 

case of condensate wells, where the separator gas and liquid exists as a 
single phase in the well pore, the gravity, compressibility and viscosity 
of the combined well stream are important. These properties are best 
determined by laboratory measurement on a sample of the gas but they may 
also be calculated from the analysis of the gas. 

A detailed discussion of the important properties of natural 
gas along with examples illustrating how they may be calculated or esti- 
mated is presented in Appendix A. Also included in the Appendix is a 
tabulation listing those physical properties of the constituent components 


of natural gas that are required in carrying out well test calculations. 


Measurement of Flow Rates 

The measurement of gas flow rates during testing operations 
{is normally carried out by means of an orifice meter or a critical flow 
prover. In rare instances, an orifice well tester, a pitot tube or 4 
choke nipple might be used. 

When a well is tested directly into a pipeline, the gas measure- 
ment is normally carried out with an orifice meter. In addition to this, 
in many instances even though gas is being flared, a portable testing set- 
up including an orifice type meter is used. The main advantages of this 


means of measurement are that it is one of the most accurate methods 


available, and that it permits testing directly into a pipeline. Proper 
installation and maintenance, and an intelligent appraisal of the data 
gathered is critical in the use of an orifice meter, A number of pub= 
lications are available which outline methods for calculating gas flow 
rates from basic orifice meter data and which include the necessary 
coefficients to correct these flow rates to certain temperature and 
pressure conditions. Among these are the American Gas Association 
Committee Report No. 3 (99), the California Natural Gasoline Association 
Bulletin No. TS-402 (102), and the back pressure testing manuals of the 
Interstate Oil Compact Commission (103) and the Kansas State Corporation 
Commission (104), The methods and coefficients published in these manuals, 
differ only slightly and in most cases have been based on work done by the 
American Gas Association. The Board accepts their Measurement Committee 
Report No. 3 as the standard reference for the installation and use of 
orifice meters. 

The measurement of gas with the critical flow prover is based 
‘on the fundamentals of gas flow through an orifice under "critical con- 
ditions" or conditions where the velocity of flow has reached a maximum 
and remains constant. These conditions are normally assumed to occur, 
and the use of the critical flow prover is considered valid, when the 
ratio of the pressure downstream of the orifice to the upstream pressure 
is less than approximately 0.5. The main disadvantage in the use of the 
ube ch flow prover is that the produced gas must be blown to atmosphere, 
It is, however, a satisfactory means of measurement and its use involves 
a minimum of equipment. The U. S. Bureau of Mines Monograph 7 first pub- 
lished the method and provided appropriate coefficients for the use of 
the critical flow prover. However, the formula proposed in Monograph 7 
was based on ideal gas behavior and it is now recognized that deviations 


from such behavior should be accounted for, Subsequently, the Texas 


Railroad Commission and others have published appropriate modifications 
to the formula through use of the gas compressibility factor. Monograph 7 
is still accepted as a standard in the design and field use of critical 
flow provers, but unless appropriate modifications are made, other publi- 
cations such as the testing manuals of the Interstate Oil Compact Come 
mission, the Kansas State Corporation Commission, or the Texas Railroad 
Commission (106) should be referred to for appropriate formulae and 
Coellicients, 

The pitot tube and orifice well tester are no longer used 
extensively in gas flow measurements. However, they do serve a useful 
purpose in gas measurement under certain conditions, particularly where 
the well produces from a low pressure reservoir, Choke nipples are 
sometimes used for measuring gas flow rates during testing operations 
when the gas must be flared to atmosphere and where the use of a critical 
flow prover or orifice meter is precluded by excessive pressure or other 
reasons. 

Each of these lesser important methods of gas measurement is 
discussed in some detail in Monograph 7, however, it should be noted 
that the basic equations included therein do not incorporate the com- 


pressibility factor, and should be modified accordingly. 


Measurement of Pressures 

The accurate measurement of pressures corresponding to the pre- 
viously discussed flow rates is also of great importance in gas well test- 
ing. Since interpretations of back pressure test results must be based 
on the theory of flow in the reservoir, it follows that the important 
pressure in interpreting the test is the reservoir pressure, Ideally, 


this pressure should be measured directly through use of an accurate 


carefully calibrated bottom hole pressure gauge. There are many types 
of such gauges available today, all of which when used properly are 
quite adequate for obtaining accurate sandface pressures directly. A 
detailed discussion of the general working principles and use of common 
gauges is available in the Alberta Oil and Gas Conservation Board's 
Manual For The Conducting and Reporting of Bottom-hole Pressure Tests 
(98). 

In some instances, due to mechanical difficulties or for other 
reasons, it is not practical to use a bottom hole gauge. In addition to 
this, for shallow wells the accuracy of conversion from the top hole to 
bottom hole pressures is acceptable, so bottom hole gauge measurements 
are unnecessary. in these situations, wellhead pressures are usually 
measured and converted to sandface conditions, as has already been men= 
tioned, The highest possible accuracy in wellhead measurement is im- 
portant, and for best results these pressures should be taken with a 
dead-weight gauge. The dead-weight gauge is an instrument which utilizes 
the balancing effect of an unknown pressure against a set of calibrated 
weights, A discussion of the operations of a dead-weight gauge is con- 
tained in the "Gas Measurement Manual" published by the American Gas 
Association (100). 

There are many methods of conversion from wellhead to bottom 
hole pressures, each having advantages and disadvantages. While the method 
originally proposed by Pierce and Rawlins and detailed in Monograph 7, along 
with its equivalents in the form of the Texas Railroad Commission method, the 
Vitter's method and the Modified Vitter’s method, is suitable under certain 
circumstances, other methods developed more recently are superior in most 
instances, A brief discussion of some of these methods is included in Appen- 


dix B, for the case of a single phase gas under shut-in and flowing conditions. 
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A less detailed discussion of the determination of bottom hole pressures from 
wellhead measurements is also included for the case where both gas and liquid 
phases exist in the well bore, 

In the case where only the gas phase is present in the well bore, 
the Board concludes that for shallow, sweet gas wells several methods appear 
to give comparable results in terms of accuracy of bottom hole calculations. 
However, if a standard method for relatively correct pressure calculations 
regardless of well depth and type of gas is to be selected, the method of 
Cullender and Smith (24) is probably superior. This method has therefore 
been adopted by the Board for converting wellhead pressures to sandface 


conditions. 


Experience indicates that regardless of the method used, calcu- 
lations of sandface pressures at high flow rates (greater than 12-15 MMcfd) 
are often subject to excessive errors, The reason for this may be the 
difficulty in accurately estimating friction effects at very high flow 
velocities. For relatively deep, high capacity wells flowing at high rates, 
the sandface pressure drawdown is normally small whereas the losses due to 
friction are very high. Consequently, a small error in determining the 
weight of gas and the friction loss can render the calculated sandface 
pressure meaningless. For these reasons, in cases where wells will be 
tested at high flow rates, serious consideration should be given to 
directly measuring the sandface pressures wherever physically possible. 

The Board also concludes that where two phase flow occurs in the 
well bore, the problems arising from the calculation of bottom hole pres- 
sures are virtually insoluble at this time, It is, therefore, recommended 
that wherever an appreciable amount of liquid phase is suspected in the well 


bore, the bottom hole pressure be directly measured, 
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Relation of Theory to Practical Well Tests 


The discussion of the deficiencies in the conventional back 
pressure test, and the examination of the more advanced methods of testing, 
demonstrate the need of a thorough understanding of the theory of gas 
flow in the reservoir in order to obtain maximum utility from well tests. 
For these reasons discussion of the theory of both steady and unsteady 
erate flow in the reservoir is presented in Appendix C. Included are 
reviews of the flow equations under both turbulent and laminar conditions, 
along with example problems to {llustrate some practical uses of these 


equations. 


2. THE FLOW TESTING OF GAS WELLS 


As noted in part 1, Pierce and Rawlins were the first to pro- 
pose and set out a method for testing gas wells by gauging the ability 
of the well to flow against various back pressures. This type of flow 
test has been designated the "conventional back pressure test" in this 
manual, 

The conventional back pressure test involves the determin- 
ation of the shut-in reservoir pressure. The well is then flowed at 
three or four different conditions, and the flow rate and the correspond- 
ing flowing sandface pressure is measured at each of these conditions. 
The results are plotted (as illustrated in Figure l-l) as the difference 
in the squares of the formation and the flowing sandface pressures versus 
the corresponding flow rates. 

This method of testing and the interpretation of the data is 
relatively simple, and the method has been considered the basic acceptable 
standard for testing gas wells for many years. For this reason, and be= 
cause an understanding of some of the problems associated with the con- 
ventional back pressure test justifies and supports the development of 
more practical type flow tests, the theory of the conventional type test 


is reviewed, 


Theory of the Conventional Back Pressure Test 


At the outset it should be understood that the original develop- 
ment of the back pressure relationship by Pierce and Rawlins was based on 
empirical methods, In the period since development of the basic relation- 


ship, much has been learned regarding the flow of fluids in porous media, 
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This now permits a theoretical analysis and better understanding of, 12 
not a complete justification for, the basic back pressure relationship. 

The relationship is commonly expressed in the form of equation 
(1-1): 


Chips ispiaa te Ga) 


&O 
"i 


where 
Q = rate of flow, millions of cubic feet per day 
at standard conditions. 
P = ghut-in formation pressure, psia. 
Pp = flowing sandface pressure, psia. 
C = a coefficient which describes the positioning 
of the back pressure test curve. 
n = an exponent which describes the slope of the 
pack pressure test curve. 
If we compare this equation to the radial laminar flow equation, 


described in Appendix C, 


0.7054 x 10°° kh toa - P.*) (c-8) 
Q = 
Upear, Mm(r,/r,) 
where 
Q = gas flow rate, millions of cubic feet per day 


at 14,65 psia and 60°F. 

k = permeability, millidarcys. 

h = thickness, feet. 

Pe = formation pressure (at exterior drainage 
boundary), psia. 

tn = flowing sandface pressure, psia. 


u = average viscosity, cp. 
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Z = average compressibility factor, 


& 
eS = average flowing temperature, cate 
PS exterior boundary radius, ft. 
os etfective well bore radius, ft. 


it is apparent that equations (1-1) and (C-8) are similar. In fact, if 


the coefficient 'C' is a constant, at a value of 
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and if 'n' is equal to unity, the equations are identical, 

If the flow in a reservoir is steady-state and if the radius 
of drainage, Tins is known to be a constant, the value of 'C' as described 
by equation (2-1) is in fact a constant for a particular reservoir 
pressure. Although flow in a reservoir is not a steady-state phenomenon, 
it has been determined that when the radius of drainage reaches the 
limits of a reservoir or when the drainage areas of two adjacent wells 
meet, "stabilized flow" occurs and although the flow is of an unsteady- 
state nature the steady-state equations approximate actual flow conditions. 
This matter is discussed in greater detail in Appendices C and D, 

As a result, it is concluded that where flow rates and pressures 
represent stabilized conditions, the back pressure relationship as de- 
veloped by Pierce and Rawlins, with 'n' equal to unity, is analogous to 
the steady-state laminar flow equation. Since flow in a reservoir is not 
normally completely laminar, one function of the exponent may be considered 
as a correction for the actual reservoir flow conditions. For completely 
stabilized flow, its value is largely a reflection of the type of flow 
(laminar, transitional, or turbulent) in the reservoir at a particular well. 

One of the major problems related to the conventional back pressure 


test is the matter of stabilization. (This problem is discussed in detail 
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in Appendix D.) It has been pointed out that the pack pressure test will 

result in the correct straight line relationship only if the values of the 

coefficient 'C' (equation (2-1)) reflected at each flow rate in a multi- 

point test are constant, Theoretically, this will occur if the effective 

radius of drainage, Lies ig the same for each flow rate, The easiest way 

to approximate this, is to flow the well at each rate to "stabilization". 

(The effective radius of drainage will then be the exterior boundary of the 

reservoir or the point at which the pressure gradients from two adjacent wells 

meet, In any case, it will be essentially the same for each flow point. ) 
In a reservoir of very high permeability, the time required to 

obtain stabilized flow rates and flowing pressures, 4S well as a stabi- 

lized shut-in formation pressure is usually not excessive. In this type 

of reservoir, a properly stabilized conventional back pressure test may 

be conducted in a reasonable period of time. On the other hand, in low 

permeability reservoirs the time required to even approximate stabilized 

flow conditions may be very long. In this situation, it is not practi- 

cal to conduct a completely stabilized test, and since the results of 

an unstabilized test can be very misleading, the more recently developed 


methods of testing should be used to predict well behavior, 


Isochronal Flow Test 

The conventional back pressure test carried out under stabilized 
conditions, qualifies as an acceptable approach to attaining the straight 
line relationship which is essential in the proper interpretation of tests, 
because it extends each flow rate over a period of time sufficient to permit 
the radius of drainage to reach the outer edge of the reservoir or the 


point of interference between neighbouring wells. This ensures that the 
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radius of drainage, Tos and the back pressure coefficient 'C' are con- 
stants, and except for the effects of turbulence, the back pressure re- 
lationship is a straight line representing stabilized flow conditions in 
the reservoir, 

If each flow rate of a multi-point test extends for any other 
period of time insufficient for stabilization, but such that the effective 
radius of drainage is the same for each point, then the coefficient 'c' 
should also be a constant and the back pressure relationship a straight 
line. 

Cullender (23) suggested that a series of flow tests at different 
rates for equal periods of time will result in a straight line relation- 
ship, and demonstrated that such a performance curve will have a value of 
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the exponent 'n’ essentially the same as that established under stabilized 
flow conditions, This method of testing is commonly referred to as the 
isochronal flow test and may be used in conjunction with only one stabi- 
lized flow point to replace a fully stabilized back pressure test. Briefly, 
the isochronal test consists of alternately closing in the well until a 
stabilized, or very nearly stabilized, pressure is reached and flowing the 
well at different rates for a set period of time, One flow test is con- 
ducted for a time period long enough to attain stabilized conditions, The 
observations reflecting the equal time flow periods are used to construct a 
back pressure line of the proper slope, and the final stabilized conditions 
are used to estimate the stabilized coefficient 'C' and thus the position 
of the curve, 

The test proposed by Cullender, is based on the principle that 
the effective radius of drainage in a given reservoir is a function only 


of dimensionless time, and is independent of the flow rate, This assumption 


is discussed in Appendix E, 
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The true isochronal test requires that the pressure between 
each flow rate be allowed to return to the stabilized shut-in reservoir 
pressure. In very tight reservoirs, it is not always practical to attain 
a& completely equalized reservoir pressure before the initial flow period, 
nor is it always practical during the test to close the reservoir in un- 
til the original pressure is regained, As a result, the true isochronal 
test proves impractical as a means of testing many wells, 

Katz (49) has suggested that a modified isochronal test con- 
ducted with a closed-in period equal to the flow period may give satis- 
factory results provided the associated unstabilized closed-in pressure 
is used for Pe in calculating the difference of pressures squared for the 
next flow rate, This method has been used for testing many wells, and 
indeed has given results which appear quite satisfactory. The method, 
to be referred to as the "modified isochronal flow test", does not yield 
a true isochronal curve but closely approximates the true curve, This is 
illustrated in Appendix E, 

More recently, several methods have been published for estimat- 
ing the stabilized performance of gas wells on the basis of theoretical 
considerations and data obtained from short term flow tests. One such 
method is the "two flow method" as proposed by Carter, Miller and Riley, (14) 
This approach involves closing in a well until the pressure is equalized, 
flowing the well at a constant flow rate for a short period of time, again 
closing the well in until static conditions are attained, and finally flow- 
ing the well at a different and constant flow rate for a time equal to the 
first flow period. Calculations are then made using basic flow theory of the 
skin factor and the turbulence factor, These factors are then used to cal- 
culate the stabilized relationship between pressure drawdown and flow rates, 
This relationship can be plotted in the form of the conventional back pres-— 


sure test curve, Winestock and Colpitts (96) have Suggested a method for 
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interpreting the two flow rate form of test which does not require constant 
flow rates. This modification reduces greatly the expense and inconvenience 
of running such a test, 

Methods such as the one described above are particularly useful 
when testing low permeability gas wells which are not connected to pipelines 
in that they estimate the stabilized performance of a well without wasting 
large amounts of gas, These more recent approaches to testing are being used 
more frequently by industry and undoubtedly will become more important in the 
future. 

The modified isochronal test is probably the best all purpose 
multi-point flow test available for gauging a well's ability to produce 
gas. For this reason, it will be discussed in some detail. It is recog- 
nized that most of the matters discussed in connection with the modified 
isochronal test are of importance in the conduct and interpretation of 
any flow test, and is therefore understood that where applicable these 
matters will be given consideration regardless of the type of test being 
conducted. (Significant differences in the theory or conduct of other 
type tests will be referred to and compared with the modified isochronal 
test.) 

Logically, the modified isochronal test (or the isochronal test) 
comprises three distinct parts. These are the establishment of the stabi- 
lized closed in reservoir pressure, the determination of the proper value 
of the exponent 'n' and the positioning of the back pressure curve to stabi- 
lized conditions, The discussion of the test will be handled therefore 


under three general headings. 


Determination of Stabilized Reservoir Pressure 
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The stabilized formation pressure is normally determined by 
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closing in the well for an appropriate period of time prior to the initial 
flow period, In many reservoirs, due to the tight nature of the formation, 
it is not practical to close the well in until the stabilized pressure can 
be measured directly. In these instances, the well may be shut-in only for 
a period of time long enough to establish a high degree of pressure equal- 
ization in the area immediately surrounding the well bore. At the time when 
this initial closed-in formation pressure is no longer building up rapidly 
(say, less than 1/10 of one per cent of the previously recorded pressure 
in 30 minutes) the modified isochronal test method should give acceptable 
results with respect to the slope of the back pressure line, If the shut-in 
period is of sufficient length, and if the appropriate observations are made, 
the true static reservoir pressure (which is required for positioning the 
line), may be calculated from the pressure build-up data. Many methods have 
been developed for determining this static pressure from pressure-time data. 
Several of these methods are discussed in some detail in Part 4, "Interpret- 
ation of Build-up and Drawdown Data", 

If the initial shut-in period is not long enough to permit either 
ene measurement or calculation of the true static formation pressure (or 
if the necessary data have not been recorded), the pressure can probably 
be determined subsequent to the flow portion of the test with little or 
no effect on accuracy, provided, this is done within a reasonable period 
of time, 

In the case of a new well in a reservoir or one which has pro} 
duced only limited amounts of gas, the discovery pressure for the well 


may be used as the stabilized formation pressure, 


Determination of Slope of Back Pressure Curve 


The most commonly used method for determining the appropriate 


os 


Slope of the back pressure relationship or the value of the exponent 'n! 
is to produce the well at a series of flow rates. In the case of either 
the modified or the true isochronal test, the flow periods are normally of 
equal duration, or at least observations are made at equal time intervals, 
(For a conventional back pressure test, each of the flow periods must be 
of sufficient length to ensure stabilized conditions, so theoretically 
they also should be of equal duration.) 

Theoretically the flow periods for an isochronal test may be 
as short as a matter of minutes (10 minute flow periods have been plotted 
and give accurate results in some pools), but from a practical point of 
view should be long enough to obtain steady conditions with respect to 
the well bore, measurement facilities and other surface equipment. One 
hour flow periods are commonly used and in most reservoirs should give 
accurate results. For the modified isochronal test, the closed=in period 
between flow rates should be of approximately the same duration as the 
flow period, and in any case should be long enough that the pressure is 
no longer building up at a rate that will affect subsequent flow tests. 
The pressure immediately prior to each flow rate must be recorded and used 
as the formation pressure with respect to that rate. (For the isochronal 
test, the closed-in period is of sufficient duration to permit the pres- 
sure to return to the initial closed-in pressure, and for the conventional 
test, the closed-in period is held to a minimum.) The results of the iso- 
chronal test are plotted in the normal manner as the difference in the form- 
ation and the flowing pressures squared versus the corresponding flow rates. 
This results in a straight line relationship with a slope equal to the re- 
ciprocal of the exponent 'n'. 

One of the problems that may arise in connection with back 


pressure testing (regardless of the method used), is the departure from 
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approximate laminar flow and the effects of the onset of turbulence as 

the gas flow rates are increased from point to point during a multi-point 
test. Binckley (7) has shown that for steady-state isothermal flow, the 
value of the exponent of the back pressure curve 'n' is 0.5 for completely 
turbulent flow, Since 'n' equals unity for completely laminar flow, the 
range of slopes to be expected for back pressure curves is 0.5 to 1.0. 

The flow rates associated with virtually all back pressure tests 
are neither in the range of complete laminar nor complete turbulent flow. 
Thus, the degree of turbulence and its effect on 'n' will increase for 
each increasing flow rate, This effect will tend to reduce 'n' and if 
the range of rates of flow for a particular test is wide enough, the change 
in 'n' from near unity towards 0.5 will result in an upwards curvature of 
the back pressure line, Binckley (7), Elenbaas and Katz (27) and others 
have shown that this curvature does exist. 

In most practical cases the curvature will be so slight that 
the back pressure test curve may be represented by a best fit straight 
line with a constant slope. However, inability to predict the degree of 
this curvature makes the interpolation of a back pressure curve to high 
flow rates very difficult. Govier (35) has presented a method for inter- 
preting back pressure tests which attempts to account for the erfect of 
turbulence, Whenever a test is to be extrapolated to flow rates well beyond 
those covered during the conduct of the test, it is recommended that cone 
sideration be given to correcting the test for the turbulence effect, 

It follows that another approach to the determination of the 
appropriate slope of the back pressure curve (although one which is not, 
part of any of the three standard type flow tests) 


» might require only 


one accurate flow point, This point is then plotted as log (P = P = 
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versus log Q, and a line reflecting laminar flow (slope of unity) is 


drawn through the point, The method of Govier (or any other appropriate 


approach), is then used to adjust this line for turbulence at several 
different flow rates. Subsequently, a best fit straight line may be 
drawn through the points thus obtained. This approach is not feasible 
unless a considerable amount of information regarding certain character= 
istics of the formation can be determined with a reasonable degree of 
accuracy, Since this information must be obtained from build-up or draw- 
down tests, and since the method itself is based on the more recently de- 
veloped theory, it is treated in more detail in Part 5, "The Estimation 


of Flow Behavior From Theory and Limited Data", 


Positioning The Back Pressure Curve 


For a properly conducted conventional back pressure test, the 
curve is automatically positioned to stabilized conditions because each 
of the flow rates is stabilized. However, for an isochronal type test, 
the final step involves positioning the curve to reflect stabilized con-= 
ditions, Prerequisite to this step is a knowledge of the static formation 
pressure, the determination of which has been discussed earlier. Normally, 
the final flow rate of an isochronal test is extended until reasonably 
stabilized flow conditions are attained. The stabilized sandface pressure 
and the corresponding flow rate are then used to establish a stabilized 
point on the back pressure plot. A line through this point parallel to 
the isochronal line represents the stabilized back pressure relationship. 

The main problem with respect to this method is establishing 
a suitable time for which the well should be flowed, or a time to approach 
a satisfactory degree of stabilization. It can be seen from Appendix D, 
that any established criterion is somewhat arbitrary, nevertheless, a 
standard is desirable, 


If it is accepted that a suitable degree of stabilization will 
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have been attained when the radius of drainage has reached the outer edge 
of @ one section spacing unit regardless of the actual density of develop~ 
ment, several of the formulae presented and discussed in Appendix D can 

be used to calculate directly the time necessary to reach this condition, 
The Appendix shows that the time required to satisfy this criterion is 
Surprisingly long. (The calculated time to stabilization may be several 
months for very tight reservoirs,) As a result, an attempt has been made 
to select an acceptable alternative and avoid having to position a back 
pressure curve on the basis of exceptionally long flow periods, Bearing 
in mind the use to be made of the flow test results, the Board concludes 
that @ maximum time to Stabilization of 15 days is Satisfactory. In other 
words, if the calculated time for the pressure-radius profile to move to 
the edge of the one section Spacing unit is longer than 15 days, the iso- 
chronal line may be positioned to reflect a "15 day Stabilization", and 
this curve should be useful in deliverability and other related Studies, 
Experience may be relied upon in extending or shortening the calculated 
stabilization time period, if previous tests in the same reservoir have 
Shown that the defined condition of stabilization can be very closely 
approximated in some time other than that calculated, 

The back pressure curve may also be positioned to represent 
stabilized conditions through calculation rather than through a 15 day 
flow test. Various methods (14)(19)(39)(70) (80) (84) (85) ee which 
permit the adjustment of the isochronal line to reflect the degree of 
stabilization corresponding to any length flow period, including that 
defined as adequately stabilized, They require a knowledge of certain 
reservoir characteristics normally determined from either a build-up or 
a drawdown test. Since the calculations are not considered part of the 


isochronal flow test, further details are included in Barbas 
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In summary, the Board believes that a proper back pressure test 
should reflect a degree of stabilization which coincides with the lesser 
of the calculated time to stabilization (adjusted on the basis of exper- 
idence) or 15 days, but recognizes that the adjustment to stabilization 
may be accomplished through either a lengthy flow test or a shorter flow 


test and gagoctated calculations. 


Frequency of Testing 
With respect to the frequency of testing, it is clear that if 


ty? 


the coefficient 'C' and the exponent 'n' of a back pressure relationship 
are constants over the life of a well, then only one back pressure test 
would be necessary to predict the behavior of the well over its life. 
Current knowledge of flow in porous media indicates that the radius of 


drainage, r is a function of the time period over which a well has been 


f? 
flowing, and that the radius and thus the value of the coefficient 'C' 
becomes fixed when it has reached the outer edge of the reservoir or the 
areas to be drained by the well in question. (The flow has then become 
stabilized.) Even when this condition has been satisfied, a study of 
equation (2-1) leads to the conclusion that changes in the average gas 
compressibility and viscosity associated with the decline in pressure as 
the pool is depleted, will result in slight changes in the coefficient 'C'. 
At the same time, it may be assumed that the effect of the depletion of 
the pool on the exponent 'n' at a particular flow rate, should be only 
through the changing gas viscosity as the pressure declines. Since the 
variation in viscosity with pressure is very small, it follows that the 
slope of a back pressure curve should change very little throughout the 


producing life of a well. This fact, and the knowledge that in the absence 


of artificial stimulation, the coefficient 'C' changes very little from 


year to year, supports the theoretical concept that only one back pressure 


test is required over the life of a well, Because of shortcomings in the 


theory and changes that may occur in reservoir properties, especially near 


the sandface, a complete test every five or six years would seem appropriate. 


The actual spacing of these tests will be influenced by many factors. Some 


of the most important ones, and their effect on a properly planned testing 


program are as follows: 


(i) 


‘eEy) 


ei) 


(iv) 


A gas well should be tested following any workover 
which may significantly change the ability of the 
well to deliver gas or if well performance indi- 

cates that an appreciable change has taken place 

in producing characteristics, 

The number of wells drilled to and producing from 

the same pool will effect the testing progran, 
Scheduling should be such that sufficient inform- 
ation is available (at least on an annual basis), 

to accurately trace the pressure history for the pool, 
The possibility that damage may occur to a reservoir 
as a result of excessive producing rates also effects 
the timing of flow tests, If for a particular 
reservoir, strict control of maximum producing rates 
is necessary to prevent well damage, the frequency 

of testing must be increased to ensure that adequate 
information regarding the behavior of the wells is 
available at all times, 

The ultimate market to which the @as will be delivered 
is an important factor, If the market being served 


requires current and highly accurate information 
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regarding the deliverability of the well, then a 
greater number of tests are necessary, 
Acceptance of the premise that the exponent 'n' will remain 

essentially constant over the life of the well leads to the use of the 
80 called "single point test". This procedure involves completion of 
a multi-point stabilized test to determine the exponent 'n' and the stabi- 
lized coefficient 'C'. Thereafter, if there is reason to believe that 
'C' may have changed, only one stabilized flow point is obtained and the 
new back pressure curve is constructed by drawing a line with a reciprocal 
slope equal to the known 'n', but through the newly established point. 
The most desirable testing program, planned to yield the maximum amount 
of technical information in the most economical manner, should include a 
series of both multi-point and single point tests, spaced effectively 


over the life of a well. 
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3. FIELD CONDUCT AND REPORTING OF DATA 
FOR BACK PRESSURE TESTS 


Earlier sections and corresponding appendices have indicated 
the Board's views with respect to the type of back pressure test which 
is recommended, the frequency of testing, and the suggested method for cone 
verting surface recorded pressures to sub-surface pressures. This section 
is intended as an appropriate conclusion regarding each of these matters, 
The discussion of field procedure will deal primarily with the conduct of 
a modified isochronal test, however, mention will be made where the pro- 
cedure must be varied for other types of test. In a similar manner, the 
discussion of the calculations related to the test assumes use of the 
Cullender and Smith method (24) for converting pressures from top hole 
to bottom hole conditions. The standard forms presented for the record-= 
ing of field notes and the calculating of results have been prepared 
specifically for use with these methods, although the form for recording 
field notes may be used for any type test. This is not true with respect 
to the calculation sheet. 

Complete results of a typical modified isochronal test are shown 
on Figures 3-1 to 3-4 inclusive, to facilitate an understanding of the 


forms, 


Field Conduct of the Back Pressure Test 

The conventional back pressure test is satisfactory only in 
reservoirs where stabilized conditions can be attained in a relatively 
short period of time. Where more practical, a form of the isochronal 
test should be conducted, 

In prescribing field conduct of the test, considerable use has 


been made of the recommended form of rules of procedure included in the 


Secs 


Interstate Oil Compact Commission back pressure test manuel (103). Details 


with respect to proper field conduct are listed below in point form. 
Reference is made to Figure 3-1, a suggested form for the recording of 


field notes. 


A. Shutein Pressure 

is The well to be tested shall be produced prior to the initial 
shut-in period for a time sufficient to clear the well bore of any liquid 
hydrocarbons or water that may have accumulated. In the case where the well 
is tied into a pipeline the flow rate shall be large enough to clear the 
well bore of liquids, If the well cannot be cleared of liquids while pro- 
ducing into a pipeline, it shall be blown to atmosphere to remove the 
liquids’. 

Where the well is not connected to a pipeline, it shall 

be blown to atmosphere until it is obvious the well bore is cleared of 
liquids, Particulars of blowing the well shall be recorded in detail as 
shown in Figure 3-1. 

26 The well shall be shut-in for a period long enough to ensure 
stable pressure conditions. If the static pressure is to be determined 
from this portion of the test, the period must be long enough to permit 
the gathering of sufficient pressure-time data to calculate the static 
pressure by one of the methods presented in Part 4 or any other acceptable 
procedure and in all likelihood need not be in excess of six days, (It 
may be recalled from the previous section that the modified isochronal 
test gives satisfactory results even when the initial shut-in pressure 
is not completely stabilized, provided the rate of pressure build-up is 
very small when the flow test begins. A rule of thumb ig that the rate 
of build-up should be less than about 1/10 of one per cent of the 


pre- 


viously recorded pressure in 30 minutes, ) 


wags 


If there is reason to believe that an appreciable accumu- 
lation of liquids may occur in the well bore during the shut-in period, 
the sandface pressures shall be directly measured with a bottom hole 


pressure gauge. 


B. Flow Periods 

1. The well shall be produced for a series of four different 
flow rates. If four significantly different flow rates cannot be obtained, 
three may suffice, (In the case of a conventional test, each point must be 
stabilized so theoretically should be of equal duration, For an isochronal 
type of test, the flow periods need not be of equal duration but at least 
one observation must be taken at the same time interval of each flow period, 
It therefore is normally convenient to make the flow periods of equal du- 
ration for all but the long or stabilized flow point.) 

Pale The lowest flow rate shall be a rate sufficient to keep 
the well bore clear of liquids. In some instances this may not be 
feasible if a sufficient spread in points is to be obtained, however, 
no appreciable quantity of liquids should be permitted to accumulate 
in the well bore during testing at the lowest flow rate, 

30 Care must be taken to ensure an adequate spread of the 
flow points. (A widely used rule is that the drawdown for the lowest 
flow rate be at least five per cent of the static wellhead pressure, 
and that the drawdown for the highest flow rate be no greater than some 
twenty-five per cent of the static pressure.) Care must be taken at 
all times, to prevent well damage by unnecessarily coning water into the 
well bore at x higher drawdown rates, 

4, An increasing flow rate sequence shall be employed, (If 


the sequence is reversed because it is impossible to obtain the required 


eee 


number of flow rates without accumulating liquids in the well bore at 
the lower rates, then the modified isochronal test method loses accuracy, 
and may not be acceptable. ) 

Do For an isochronal type test, pressure observations shall be 
taken during each of the flow periods at equal time intervals. Each flow 
period must be long enough to ensure steady gtate conditions with respect 
to measurement and other surface facilities, In no case, however, shall 
they be less than one hour, One flow period shall be of sufficient du- 
ration to satisfy the accepted criterion of stabilized flow, or be of 
sufficient length that the appropriate observations may be made to permit 
a calculation to reflect stabilized conditions. By taking appropriate 
observations at time intervals consistent with observations made during 
other flow periods, the extended flow rate may serve as a point in estab- 
lishing both the slope and position of the line. 

For a conventional back pressure test, each flow rate 
shall satisfy the established criterion of adequate stabilization. 

Os For a modified isochronal test, the closed-in time between 
each flow period should be of approximately the same duration as the flow 
period and in any case shall be long enough that the pressure is no longer 
building up at a rate that will appreciably affect subsequent flow tests. 

For an isochronal test the shut-in time shall be sufficient 
to permit the pressure to return to a pressure comparable to that prior to 
the initial flow period. 

For a conventional test, the shut-in time shall be held 
to a minimum, 

fer In cases of deep wells with high flow rates, and particu. 
larly where two phase flow is suspected, wherever possible the flowing 


and closed-in sandface pressures should be directly measured with a bottom 


hole gauge. In other cases, the flowing and closed-in pressures shall be 
measured with a dead-weight gauge and recorded on Figure 3-1, at least 
twice during each flow and closed-in period. One of the recorded pres- 
Sures shall be at the end of each period, 

8. If bottom hole pressures are to be calculated, the gas 
temperature at the wellhead shall be recorded during the closed=in and 
each flow period. To facilitate this a temperature well should be in- 
stalled in the wellhead. 

9. The gas flow rates during each of the flow periods shall 
be measured by an orifice meter, critical flow prover, orifice well tester, 
or other suitable metering device in good Operating condition, The measure- 
ment, and subsequent calculation of flow rates, shall be in accordance with 
the procedures set out in the publications referred to in Part l (72)(99) 
(100) (101) (102)(103)(104)(106), 

10. Where wellhead separation takes place, in addition to 
measurement of separator gas as described in item 9, the separator liquid 
producing rate shall be accurately measured and recorded along with the gas 
flow data on Figure 3-1, The measurement of the produced liquids shall be 
such that it permits the accurate determination of the producing gas-liquid 
ratio for each flow rate, 

ll. Where the bottom hole pressures are to be calculated, the 
specific gravity of the produced gas, as well as liquids where separation 
takes place, shall be accurately determined and recorded, This may be done 
by direct measurement using proper techniques, or alternately by careful 
sampling, analysing and the necessary computations. The determination of 
gravity need be carried out during one flow period only, preferably the 


second or third, 
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The Caleuleting, Plotting end Reporting of Test Rebults | 

This discussion assumes that the Cullender and Smith method 
(see Appendix B) of calculating sandface pressures from wellhead ob= 
servations will be employed. As & result, the suggested form of calcu- 
lation sheet (Figure 3-2) to which the discussion refers, jg designed 
specifically for use with this particular method and is not applicable 
to other methods. Since each of the methods of calculating sandface 
pressures from wellhead measurements bear some similarity to each other, 
the rules set out as standard for use with the Cullender and Smith 
method should be followed where applicable when other methods of calcu- 
lation are employed. 

Figure 3-e serves for converting wellhead observations to 
sandface pressures, SO that it is unnecessary to complete this form 
when bottom hole pressures are measured directly. 

Once again, the specific instructions for calculating the 
‘gandface pressures are presented in point form, 

lig The effective diameter, d, of the flow string is the 
internal diameter of the string through which flow is taking place. In 
the case of annular flow an equivalent diameter is calculated from 
equation (B-31). 

2, The length of the flow string, L, is the distance from 
the wellhead to the mid-point of the producing perforations. 

Bi In cases where the total well effluent is separated into 
a gas and a liquid at the wellhead, but existed as a gas in the well 
bore, the gravity of the recombined effluent is determined. This may 
be done through use of equation (A-4) and Figure A-9, as illustrated in 


Example A-3, all in Appendix A, If an analysis of the recombined well 


effluent is available, the gravity may be determined from the molecular 


ans77_ 


weight of the mixture, as illustrated in Example A-1l. 

4, The pseudo critical properties of the well effluent are 
determined from the recombined analysis by calculation as in Example A-1; 
or from the gravity of the effluent using Figure A-1. 

Do The formation temperature, Toe is based on measured temp- 
eratures in the same reservoir wherever possible. If no such measurements 
exist, the formation temperature may be estimated on the basis of similar 
reservoirs in the immediate proximity, or on the gradients suggested in 
Figure 3-5. 

6. For the static case, the wellhead temperature, T, is 
taken as the mean annual surface temperature if the closed-in period has 
been extensive (greater than 24 hours). If the closed-in period has been 
short, such as will be the case in a modified isochronal test, it is prob- 
ably more accurate to estimate a wellhead gas temperature closer to the 
final flowing gas temperature immediately prior to shut-in. Figure 3-6 
presents mean surface temperatures for the Province of Alberta, 

In the flowing case, the wellhead temperature is measured 
directly. 

ie The gas flow rate calculations are carried out in an 
approved manner, depending upon the method of measurement. Where liquids 
have been produced at a separator, they are converted to an equivalent 
quantity of gas. This may be accomplished by using the vaporizing volume 
ratios as presented in Figure A=-9. The gas-equivalent of the produced 
liquids is then added to the separator gas to get the total gas production. 
It is upon this total that the friction determination is based, 

8. For the static case, and in the absence of a known temp- 
erature profile, the average temperature, tT. is taken as the arithmetic 


average of the mean annual surface temperature and the formation temperature. 


Bic. ee 


In the flowing case, the average temperature is estimated as the average 
of the recorded flowing wellhead temperature and the formation temperature. 

9, The determination of the compressibility factor, Z, is 
carried out by one of the methods 4llustrated in Examples A-1 and A=-2, 

10. For purposes of estimating the viscosity (for friction 
factor determinations), an average well bore pressure and temperature 
(for all flow rates) may be estimated from a knowledge of the wellhead 
pressures and temperatures, the formation temperature, and an estimated 
formation pressure. The viscosity of the gas at these average conditions 
can be determined from Figures A-10 to A-12, as illustrated in Example 
A~4, This average viscosity may be used in the determination of the 
friction factor for all flow rates without introducing appreciable errors. 

11, The relative roughness, 6/d, of the flow pipe is deter- 
mined from Figure B-2. In the absence of special knowledge of the pipe 
condition, the absolute roughness, 6, may be taken as 0.0006. 


12, The Reynolds Number, N__, for each flow rate is calculated 


Re 
from equation (B-18). 

13. The friction factor, f, for each flow rate is determined 
from Figure B-l, 

14, The friction effects as included in the Cullender and 

fe ; 
Smith equation, F , are determined for each flow rate from equation (B-22), 
| 15. The values of Is I. and I, are calculated from equation 

(B-5) in the static case and (B-27) in the flowing case, 

16, The trial and error calculations of P, and Po (from equations 
(B-8) and (B-9) in the static case and (B~29) and (B-30) in the flowing 
case), are repeated until the calculated pressure values are within one 


psi of the assumed values. These calculations are illustrated in Examples 


B-1 and Be2 o 


a0: 


17. In both the static and flowing cases, the final pressures 
are corrected using Simpson's rule, equations (K-10) or (B28), 

The results of the back pressure test shall be recorded on 
Figure 3-3, the summary sheet. The producing rates and the pressure data 
(where sandface pressures are calculated), are recorded under item A, 

If the sandface pressures have been measured directly, they are recorded 
under item B. The shut-in pressures are recorded under item C. -4in the 
case of a conventional back pressure test, only one shut-in pressure will 
be recorded.) The plotting data, item D, are made up of the appropriate 
shut-in and flowing pressures and the total gas flow rates. These re- 
sults are plotted as Q versus chee - Be) as illustrated in Figure 3-}, 
The plot should be made on equal scale log-log paper, and a straight line 
drawn through a minimum of three points, If a straight line is not in- 
dicated by at least three points, consideration should be given to re- 
testing the well, The reciprocal Slope of the line is the exponent 'n', 
If the value of 'n' is greater than 1.0 or less than 0.5, consideration 
should be given to retesting the well, unless experience with wells in 
that pool indicates that a different 'n' value would not be obtained, 

If a well has been retested, and the test is still unsatisfactory, 
the best fit line may be drawn through the points of the test which appears 
to be the most acceptable, If the resulting 'n' is greater than 1.0, a 
Jine reflecting an 'n' of 1.0 shall be drawn through the highest flow rate 
point. If the 'n' is less than 0.5, a line reflecting an 'n' of 0.5 shall 
be drawn through the lowest flow rate point. In any case the relation- 
ship resulting from this second unsatisfactory test is only an estimated 
one, and consideration should be given to a retest within a one year period, 
The retest should involve alterations in test procedure in an attempt to 
obtain a satisfactory relationship. This change may involve direct sand- 


face pressure measurements if two-phase flow appears to be a possibility, 


uO! 


or may lead to another type of flow test. 
In the case of an isochronal test, the line must be positioned 


to reflect stabilized conditions. This is done by plotting the stabi- 


lized value of (as - P.*) versus the appropriate flow rate. A line of 
reciprocal slope 'n' is drawn through the point, as is illustrated in 
Figure 3-4. The data corresponding to this flow point may conveniently 
be designated in the field notes, or on the calculation sheet and the 


summary sheet, as run number 4(a) for a four point test. 


FIGURE 3-1 
OIL AND GAS CONSERVATION BOARD 
BACK PRESSURE TUST — FIELD NOTES ~ PAGE 1 

Well Vea p ta Location _USd.-iwp.-ivre.-''cr. pate fhugust 1, 1003 
Field us “iels Poo] 2718 Pool Tested By 2). b.Omith 

Casing Size 1/) Weight 10-4 Grade _J=2) Producing Interval __0 (32-0) (40 

Tubing! Size.-i ss: Weight — 252 Grade —J=20" Producing Through ubing = 
Well Blown 1 minutes. Description of Blow Water-Condensate Unray, Clear in 5 min. 
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OLL AND GAS CONSERVATION BOARD 
BACK PRESSURE TEST = FIELD NOTES - PAGE 2 


SHUT=-IN PERIOD NO. 3 (Intermediate - for use with isochronal type tests) 
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OIL AND GAS CONSERVATION BOARD 
BACK PRESSURE TEST = FIELD NOTES - PAGE 3 


SHUT-IN PERIOD NO.» (Intermediate - for use with isochronal type tests) 


Pavan, Shut-In Time Tubing Pressure | Casing Pressure Wellhead Temp. 











Date | Time Hrs. bE 
ta emer en emt) ORG T Laagaal | ha Pex | ne 
ec eee ele: Teer Nemo sy Sls me + [avi be | 106 
‘Chapa oe wenn fag en gogo]. oO 105 
SD UL el ee ee eee eee 
indecent alle Melia eclnecnee, Careers lice MA ee 
FLOW PERIOD NO. hi() 
Duration of Run_l!'thrs. Orifice Size 2-000 Meter Run or Prover Size__ te Veb 1.026 
High Pressure Separator - Operating Pressure___90 psig Temp. 90 oF 
Low Pressure Separator — Operating Pressure _________ psig OMY 6. eee 
Date | Time sig psig oF psia oF Prod-sbls_ | Prod-Bbls 
ET IS OY LS) Be 
Augs@o[ueO0pal os | eo | LO | arf ed 
SE EL NN SN CL TT a 
Coes RSL Sat SE SS OT Le OE 
Riera OOne| pee edie leedeene| LOL | 900 | 3] ob 
ree aae0Opa|10nsfetUsi (ee de [e101 [900 [ay [Be 
Aug. 1 [uz00pm] 1h [1029 |__| 100 [900] 37 Total 30.0 
a sserjeshindiec'* 0) "2h eg ag Bd Co 
SHUT-IN PERIOD NO. (Intermediate - for use with isochronal type tests) 
Shut-In Time Casing Pressure Wellhead Temp. 
Date | ‘Time Hrs. ee 
| Le ae a ee ee 
ol PE a ee 
en eee | ee | 
re eee ee ee | 
To"! Sa Te See ee ee 
DUT SPS ks. oo s) il ee) ee ee ee 


FLOW PERIOD NO. 


Duration of Run nrew es Oritpce.ol ze go Meter Runior Prover, size 
High Pressure Separator = Operating Pressure ________ psig Tenn.) 
Low Pressure Separator - Operating Pressure ________ psig Tethoc ee ee 





Wellhead Readings Meter or Prover Data 


Tubing | Casing Temp. Pia hiw Temp. {Condensate Water 
psig psig ok psiu Prod=-bbls_ | Prod-Bbls 





Date Time 
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FIGURE 3=2 


OIL AND GAS CONSERVATION BOARD 


BACK PRESSURE 


TEST = CALCULATION SHEET = PAGE 1, 


Well __ ABC Gas Company No.1 Location _Usd.-Twp»-ive ler, Pleld ing ht 
ee ee ee ese by os us ito Cakleulated by ri r¥Oe 
Cag. SizeS 1/?Wt. 15.5 Grade J-55Producing Throughlubing Length of Flow String_@,/o _ 


Tog. Size_3- Wt._9.3 Grade uwleS5Effective Diameter of Annulus 





Perforati ons ssh seen760. SiGas Uravity 0.691 ~~ Liquid Gravity o.748 


Gas-Liquid Ratio, Rce_154,000. 











Properties of Mixture - GravityQ.Z10 Pc 67) Te 378 


2. ©, li sis-& 2 f 2 
es ay ee 0, Ok u,8 O Ni, 17 fee Analysis-Sample from 2nd flow 
Mean Annual Surface Temp.__37 °F, Reservoir Temp. 103 °F Source _Vemp.Lomb 
STATIC CASE toe Uy P 


0.0375GL Res f 
= w 
paier| “HR re psia pele 


89806 pote tobe beet O.zie} 2h0 
ie Ger cena { Caetene Cumcoa tt be eeemcadt Roem 
poss [531 b.714 b.loed 2193 br7.0 b.786 Dp. 2066} 205 [623 b. 837 Jo.216$ 2hon 
he73 [shy b.7hh b.2064 2187 685.0 b.o7 b.27 of esys {623 637 fo. 2117 239! 
hori [au harrap-nzod iri hoa.s p.toy prolool a3 To D8: .218$ 2382 
8p hoca [505 p.775 b.222u 2172 b94.0 b bor. 2207] 2376 [624 p37 219 2376 
Mittmnmmar fo | 
a I ae ee ee ee ee ee ee 
Ve J he Ee ee ee a ee Ee ee 
i LOWING CASE Viscosity ,y—0:0104 pvsolute Roughness ,6__0.0004 Relative Roughness ,6/d _0,00()2 
vey ae crazy) 
I= prve/ |» + “T000 






Friction Factor 





*Reynolds Number eee 8 Effects 


























Flow = TOCaLI P0011 G Oalef-( Figure B=1) ease, DUO, £ Un 
N SET Te ae Fr = aS : ° 

No. Re odio R 
0.0036 0.00088 
0.003 0.00308 1[2 “310 — 
fetes [ast cua [oust 1 o:oatss Past ; 
0.0101 -2253boue 
eee ee 0.00068 : 





- psia 
ts Pp 


Corrected 


= tt Hw 
n 
22) 





502 _ [0.790 [TS _| eee Oe aro 720 
Parr [1937] 596 [0.819 z Pb | e152 
7.535 O.on5 | ue | 19e3 
Se BS SPE SL (CPL TT 


“For annular flow refer to equations Belba, bella, b-23a & U-2ha. 
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FIGURE 3-2 
OIL AND GAS CONSERVATION BOARD 
BACK PRESSURE TEST = CALCULATION SHEET - PAGE 2, 


CASE 1. FLOW THROUGH PROVER 
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CASE 2. FLOW THROUGH METER RUN 
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FLOW THROUGH METER RUN (Continued) 
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FIGURE 3-3 
OIL AND GAS CONSERVATION BOARD 
BACK PRESSURE TEST - SUMMARY SHEET 


Well __ABC Gas Company Now] Location _Lsd.-'lwp.—Rge.-ler. Fiela_Gas Tield 
POOl tice OC ae eee Nectial by J.L,Smith Calculated by —2.1'.Doe 


Gasing Size - 1.D, Tubing Size 2:20? I,p,_ ss «0D. Produced Thro Tubing 
Perforations _____ 732-0760 Length of Flow Stine Cie eee 
Properties of Well Pri iuenteop,. CVAVity O.710" -°§ Pel 47T Tie oe oe ee 
\. FLOW TEST DATA 
R 


Calculated Gas Flow | Condensate | Total Flow 
Tubing Casing Sandface Press. | 
un No. osia osia t j ] . 











b. BOMB MEASURED PRESSURES Dat in 


Flow Rate Corr. Press. Run | Corr. Press, @ 
No. Date Time Dept Depth = psia Datum = psia Remarks 





C. SHUT-IN DATA (Note: For regular back pressure tests only one shut-in required, ) 


Shut-In Shut-In 
No. TimeeHrs, 
t : 


| Fae a de ad 





Db. PLOTTING DATA 








t Q Total 
psia MMcfd 

Frye eee 

QO) 8.95 

ERE 6.510 a 

536 0.845 5.0h 107? Q 
A.O.F. 2026 ss Med, 1 Oe ee Lh : Note: C= ae SR 
; Oh e - P 

Total volume of gas produced during fest — 202 MMef. ( f Ss ) 
Pemarks: The figures under Gas Flow Rate in Part A,renresent metered gas plus 






‘pus flashed from the stock tank at 0.36!) Mscf/STBb), 
~ | The vapour volume of condensate is 0.420 Mscf/STBb], 
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Q, MILLIONS OF STANDARD CUBIC FEET PER DAY 


FIGURE 3-4 — MODIFIED ISOCHRONAL BACK PRESSURE TEST PLOT 





»@ an 


; a - oP om - - = vad 

; * 

: —* — A— o ; . ‘ 
a a See el i < 


i nr ag - _ . | 
; ae a Pa a OR el rge  e 





- 


| 
e 
cn iene ™ : 
r 
$ 



























































































































































































































































CES (D1 2SR Re GWeY INE ee Ce ee ANSE oa ea ee 
QoS IS (SER Se SER ES AR ERS ERE R REET RRA ERLE AR 





























Suseseedesfoztesesettessasten oafonsestentastan BE RRERRE TS URRRRRRREREREEERERRARERREAAYAO 








1334 30 SONVSNOHL ‘Hid30 


OF 


TEMPERATURE , 





FIGURE 3-5 — TEMPERATURE GRADIENTS - PROVINCE OF ALBERTA. | 
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FIG.3-6— MEAN ANNUAL SURFACE TEMPERATURES 
PROVINCE OF ALBERTA. 
(Redrawn from Beach, 5 ) 
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4, THE INTERPRETATION OF BUILD-UP AND 
DRAWDOWN DATA 


Build-up and drawdown tests at gas wells are particularly use- 
ful in connection with productivity tests. Interpretation of these tests 
permits determination of the effective values of the combined reservoir 
and gas properties which are required for solution of the theoretical 
productivity equations. Moreover, the build-up tests may be used to 
determine the true static formation pressure from pressure-time obser- 
vations taken over time periods shorter than would be required for the 
well to build up to the stabilized pressure. The stabilized pressure 
is necessary in the plotting and interpretation of any back pressure 


test. 


Pressure Build-Up 
The theory of pressure build-up is based on the unsteady-state 


flow theory for the multiple transient case. Equation (C-38), from 


Appendix C, may be written as follows: 


Eat - Pe ro 
tee neo oy) Pee (Ge ed) Po oe 

2 iat De Py ts ne ema 

rE 

a 

ipso ey 
7h aps A (Cc-38) 
where 
ty = total elapsed time since the first flow 


rate began. 
t = total elapsed time since the second flow 


rate began. 


<.:60 2 


t = total elapsed time since the nth flow 


rate began. 


QD, = first dimensionless flow rate. 

does nth dimensionless flow rate. 
n 

P. = dimensionless pressure change number evaluated 
i 

at time tie 

P. = dimensionless pressure change number evaluated 

n 


at time t_. 
n 
The proper manner of conducting a pressure build-up test involves 
prcducing a well for a considerable period of time at an essentially con- 
stant flow rate, The well is then closed-in and pressure-time observations 
are recorded, If the constant producing rate corresponds to a dimension- 
less flow rate of Qn » then dy and following Wnts to qp are equal to zero 
il. n 


2] 
and equation (C-38) when referring to the sandface reduces to 


Gi i Sia (4-1) 





As the pressure in a well builds up toward the shut-in pressure, there 
is a period during which the reservoir boundary effects are not felt and 
the reservoir acts as though infinite. Under these conditions, the di- 


mentionless drawdown, Pio may be expressed in the form of equation (C36) 


P ee l/2Cin t 


2 aut 0, 80907) (c~36) 


where 6, = dimensionless time, 


Combining equations (C-36) and (4-1) 


ees 


Ip 
= 1 
een ln to. dete) (is?) 
2 





Since ty corresponds to the total elapsed time since the first 
1 
flow rate began (which is te + At), and ty corresponds to the total 
2 
elapsed time since the second flow rate began (At), equation (4-2) may 


be written as follows: 








+ 
Bae aah se uo BUC h ela ae | (4-3) 
m ci 5 
At 
where 
ig = measured sandface pressure at any time /P Peapsia: 
m ee 
Py = shut-in formation pressure, psia. 
qd = dimensionless flow rate corresponding to the 


constant flow rate prior to shut-in. 
ty = period of flow prior to shut-in, hours, 
St = shut-in time, hours, 
From the form of equation (4-3), it follows that if the squares 
of the observed pressures, Ba are plotted against the log of (t, + At)/dt, 


a straight line may be expected with the following slope: 





slope = “pf (4-4) 
z 


If equation (C-33) for the dimensionless flow rate, is combined 


with equation (4-4), then the 


Lehe x 10°u a2 TQ 
aaa 
Ek 


slope = 


(4-5) 


The permeability-pay product may be calculated from the slope 
of the straight line portion of the plot of te versus log (t . + At)/At, 


using the following equation: 


Coon. 


5 r 
° fon T 
ye 1 Th re IMO) he a at (4-6) 
slope 
where 
k = £xpermeability, millidarcys. 
h = £4pay thickness, feet. 
Pee average viscosity, centipoise. 
Ze = average gas compressibility factor. 
To = average temperature, aRe 
Q = #£=gas flow rate, millions of cubic feet 
per day at standard conditions of 14.65 
: fe) 
psia and 60 F. 
2 
FP. 
and slope = change in m Per cycle 


2.303 

One problem which arises in processing pressure build-up data, 
is the determination of the proper value of the time of flow, tos and the 
flow rate, Q, that should be used where a well has produced at various 
Reece over a considerable period of time. In this situation, it is recom- 
mended that the stabilized flow rate be taken as the constant rate of flow 
just prior to shut-in. The "pseudo flow time" should then be determined 
by dividing this flow rate into the total production obtained over the 
time interval being investigated, (Where the insitu permeability is being 
calculated and the time of stabilization as calculated from equation (D-16), 
is less than the pseudo flow time, the time of stabilization should be used 
as the length of the flow period, tipo This is not the case though, when build- 
up data is being used to calculate the static reservoir pressure, ) 

A further study of equation (4-3), would indicate that for an 
infinite reservoir, if the straight line were extrapolated to the end of 


an infinitely large shut-in time period, log (t, + dt)/At = 1.0, the static 


=ao 3 - 


formation pressure would be obtained, It is seldom that the pressure 
build-up of a well reflects infinite reservoir behavior during the final 
Stages of build-up, and as a result this approach for determining the true 
static pressure is useful only in the case of a single well located in the 
approximate centre of a large pool, from which little production has been 
taken. Further development of this theory towards the determination of the 
static reservoir pressure in presence of boundary effects is complex, ‘Many 
authors have presented methods applicable under particular circumstances, 


some of which will be discussed later in this section, 


Pressure Drawdown 

A pressure drawdown test consists of measuring and recording the 
declining sandface pressure as a function of time, essentially from stabi- 
lized shut-in conditions, as the well is produced at a constant flow rate, 
The method is particularly useful early in the producing life of the first 
few wells completed in a reservoir, however, it may also be applied to 
wells which have been shut-in for a considerable period of time, 


If equation (C-3)) 


1 rs wept es 
Poop (C=34) 
t 52 
2p f 


is combined with equation (C+36), the following equation results 


2 
Pa 2 rs gor 
P = P= 0.4045 dy Pe SSN age ol th 


Equation (4-7) may be written as 





Pp * = constant = 5 D 
(4.8) 


="Ol = 


If equation (4-8) is combined with equation (C-32) the relationship 


of ty to t 
/ =i 
2.634 x 10 ‘k Pt 
Pabst ty eae cD ee (C-32) 
i 
where 
k = permeability, md. 
Be = average pressure, psia. 
t = time, hours. 
te oe average viscosity, cp. 


> = gas filled porosity, fraction. 


r 


> effective well bore radius, feet. 


the following equation results 


QW Pe 





P e = constant - 
m 


(4-9) 
The form of equation (4-9) leads to the conclusion that a plot of the ob- 
served sandface pressures squared, Past versus the log of the time of the 
observation, t, will result in a straight line relationship with a slope 
equal to dy Date Combining the value of the slope and equation (C-33), 


the following relationship results. 


2 
‘oe LO Weer 8 


kh = 





(4-10) 
slope 


In this case, the flow rate, Q, should be taken as the aver- 


age flow rate during the drawdown test. 


Determination of Static Reservoir 
Pressure from Build-Up Tests 


Many methods are available for using pressure build-up data 


- 65 - 


to estimate the static formation pressure, Some of the better known and 
earlier papers on this subject are those of Muskat (62), Horner (40), 
Arps and Smith (2) and Miller, Dyes and Hutchinson (61). Many other 
papers have been presented which suggest modifications to the various 
methods of interpreting pressure-time data or which investigate special 
problems encountered in conducting build-up tests (26) (44)(52)(58)(59) 
(65)(71)(89). In addition, many authors have discussed the determination 
of well bore effects (commonly referred to as skin effect) from build-up 
and flow data (3)(18)(43)(76)(80). These effects are very important in 
understanding back pressure tests. 

It is not considered appropriate to deal here in full detail 
with the many methods for the conduct and interpretation of build-up tests. 
For illustrative purposes, two of the basic methods are presented, They 
are those proposed by Muskat and by Horner (as modified by Tracy). Applic- 
ations of the various methods to many test data have led to the conclusion 
that each method has advantages under certain circumstances, These circum- 
stances may be related to boundary conditions, formation geometry and homo-=- 
geneity, rock characteristics, and to general conditions prior to and during 
the conduct of the test. Given adequate data, it is believed, thet either 
the Muskat or Horner method should give reasonable results with respect to 


interpreting any pressure-time data taken under normal circumstances, 


The Muskat Method 

The Muskat method has been developed on the basis of the theory 
of the unsteady-state flow of a fluid from a bounded reservoir towards a 
well, starting from the equivalent of equation (C-34), It involves a trial 
and error solution in which the static pressure is guessed, and the log of 


the assumed static pressure minus the measured pressure is plotted for 


bolder © 


each pressure-time observation, If the choice of static pressure is correct, 
the resulting plot will be a straight line. If the estimate of the static 
pressure is too high the plotted points will curve upwards and if it is too 
low the plotted points will curve downwards. (The method was developed for 
use with oil wells and when applied to gas wells should be in the form of 
a semi-log plot of the difference of the pressures squared versus time, ) 

Larson (52) has shown that although the true equation for the 
build-up curve of a bounded well is not the same as the equation of the 
Muskat curve for early times, it does become linear and identical to the 
Muskat curve at dimensionless times greater than 0.08. Accordingly, the 
method should only be employed with data representing dimensionless times 
greater than 0.08, Larson also modified the method somewhat to permit 
the determination of the permeability-pay product of the formation, 

The main advantage of the Muskat method (with respect to deter- 
mination of static pressures), is the relatively good accuracy which can 
be obtained with a simple calculation. Other advantages are that it does 
not rely on an extrapolation of the pressure-time curve, and that the 
time interval for which the method is valid may be easily calculated, 
The fact that the method can be applied only for dimensionless times 
greater than 0.08, results in an answer which is not affected by after- 
production which has been shown to exist only for dimensionless times less 
than approximately 0,02 (61), (After-production is the continued production 
into the well bore due to the compressibility of the fluids in the well bore 
and results from closing a well in at the surface rather than at the sand- 
face. ) 

A major disadvantage associated with the Muskat method is that 
the pressure-time data must be obtained over a longer shut-in period than 
with several other methods. Also, the method does not appear to be practi- 


cal for use in very tight reservoirs. In these situations the very slow 


w OF) = 


nature of the build-up makes it difficult to judge whether the plot corres- 
ponding to a particular assumed pressure is truly a straight line, 

The recommended procedure for conducting a build-up test so that 
the data can be used for determining the static pressure by the Muskat method 
is as follows, 

AA Prior to shut-in, produce the well for a reasonable period 
of time at approximate stabilized conditions. 

20 Calculate the time corresponding to a dimensionless time 
of 0.08, from equation (C-32), The permeability used in this calculation 
should be estimated from core analysis and/or knowledge of the formation. 
Upon completion of the build-up the accuracy of the estimate may be checked 
against the actual effective permeability, 

30 Record pressure observations during the initial stages of 
the build-up. This will enable one to see the upwards curvature which 
always occurs for dimensionless times less than 0.08, and will help with 
the final interpretation of the plot. 

At dimensionless times greater than 0.08, regular pressure 
observations must be recorded, Six to ten data points should be taken 
over a period equivalent to a total dimensionless time of about 0.16, 


4, Estimate an appropriate static pressure, P Thestiyst 


f° 
guess may be based on previous knowledge of the formation pressure and on 


the behavior of the recorded build-up. Calculate values of the estimated 


2 2 
static pressure squared minus the measured pressure squared, Pe - PS . 


for each time of observation, and plot these values on the log axis versus 
the time on the linear axis (see Figure 4-5), If the portion of the line 
beyond ty = 0.08 curves downwards a higher estimate of the static pressure 


is made. If it curves upwards a lower estimate is made. This trial and 


error procedure is repeated until the plotted relationship becomes a 


- 68 - 


straight line. The final assumption will reflect the correct static for- 


mation pressure. 


The Horner Method 

Horner has developed a method for determining the static pressure 
of an oil well by plotting the observed pressure against the log of 
ea At)/Mt. When this method is modified to apply to gas wells, the plot 
is identical to that discussed earlier in connection with determining the 
permeability and pay from 4 build-up. (The modification of the method for 
a gas reservoir is based on work by Aronofsky and Jenkins (1) who compared 
the behavior of the pressure transient for 4 gas well to that of an oil 
well.) 

The basic theory of the Horner puild-up method is similar to 
that developed earlier in this section, except that 4 correction for boun- 
dary effects expected in a finite reservoir has been included, It has 
been mentioned that the infinite reservoir case may be applied to a new 
well in a finite reservoir. This method now deals with wells in reservoirs 
from which production has been taken, In developing the method, it is 
assumed that a previous pressure survey has been conducted on the well 
under consideration. In this case, where boundary effects are felt, the 
equation of flow (as developed by Horner and modified for gas wells) may 


be expressed as 





fe 
a Emenag (Spe ae et AD 
2 Vas eet 
ie ono 5 (in Roe (u)) ei) 
ae = the sandface pressure measured during 


a build-up test, psia. 
P = the pressure at the time of the previous 


survey, PpSire 


ach a 


a variable which depends on reservoir and 
flow characteristics (reservoir permeability, 
porosity, pressure; gas temperature and vis- 
cosity; and on the radius to the external 
boundary ). 


a function of u defined as follows 


y(u) = Ey (-u) + -u 


€ 


cfr 


Horner presented a plot of y(u) for a large 

range of values of u. Tracy (89) has modi- 

fied this plot as presented in Figure 4-1, The 
modification by Tracy is actually a plot of u versus 


y(u)/2.303. 


It follows that if the square of the observed pressures are 


plotted versus log (t + At)/Adt, equation (4-11) may be expressed as 


follows 


p* 


where p* 


- slope of the plot)y(u 
2.303 (ho12)3 


the extrapolation of the straight line portion 
of the plot to log (t = At)/At equal 1.0. In 
the case of a finite reservoir this is a false 
value of the final static pressure (thus desig- 


nated as P*), 


From a knowledge of P*, P? and the slope of the plot, y(u)/2.303 can be 


determined. Values of u corresponding to y(u)/2.303 are shown in Figure 


hol 


For the case of an infinite shut-in time and where the square 


of the observed pressures are plotted versus log (t + At)/At, equation 


(4-11) may be written as 


2) 7O%- 


P Ze) = > 2 slope of the plot (4-13) 
fro Dae Dae 2.303 u 


It follows that the correct static reservoir pressure, Pes may 
be calculated from a knowledge of the pressure at the time of the pre-=- 
vious survey, the slope of the plot of the observed pressures squared 
versus log (t. + At)/At, and the previously determined value of u. 

The recommended procedure for conducting a build-up test and 
calculating the static formation pressure by this method is as follows. 

1. Flow the well at a constant rate for some time prior to 
shutein, 


2. Calculate the producing time, t,, by dividing the total 


f* 
production over the time interval considered, by the constant flow rate 
prior to shut-in, The interval being considered is the time elapsed since 
the previous pressure survey. 
3. Record pressure-time observations during the shut-in 
period. 
4, Plot the squares of the observed pressures against the 
log (t, +o At)7 At 
5. Extrapolate the straight line portion of the plot (from 
which the permeability-pay product may be determined) to log (t, + At)/At = 1.0 
to determine pee 
6. Calculate the slope of the straight line portion of the 


plot ° 


7. Calculate y (u) from equation (4-12), 
26303 


8. Determine u from Figure 4-1. 
9. Calculate P, from equation (413), 
The Muskat and the Horner methods for determining the static 


formation pressure; and the determination of the permeability=pay product 


from build-up or drawdown data are illustrated in the following examples 


a aris 


Example Tee 


A new well in a newly discovered field is produced for 48 hours 
and then closed in. The cumulative gas production is 17,110,000 cubic 
feet, and the rate of flow immediately prior to shut-in is 8.5 million 
cubic feet per day. Drawdown data are recorded during the flow period, 
and build-up data recorded for the first 24 hours after shut-in. (These 
represent actual data taken on a gas well in the Province of Alberta. ) 

If the following reservoir and flow characteristics are known, 


calculate the permeability-pay product from both the drawdown and build- 





up data. 
ee = 0.02 cp 
rd =~ WV ielyate' 
a 
fe) 
Th = 560 R 
Drawdown Data 
Time of Flow Measured Pressure fBele x 107? 
t - hours Raaenets (psia*) 
Gar 983 9.66 
Ose 958 9.18 
eet 944 8.91 
0.6 918 8.43 
1.0 899 8.08 
Sie 885 TaO3 
Brae) 858 7136 
6.0 833 6.94 
10.0 814 6.63 
16.0 798 Gest 
26.0 780 6.08 
36.0 768 5.90 
48.0 758 50 15D 


A plot of mi versus log t is then prepared (see Figure 4-2). 
From equation (4-10) 


? 0,02)(0.88)(560)(8. 


(144,000) /2. 303 
kh = 951 md-ft. 


kh ¢L5x 10 


Pens 


Build-up Data 





Time of Shut-In Measured Pressure (P Va x 107° te + At 
2) Lae OD 
At = hours Po” psia (psia ) At 
e 1.004 25 
3 1009 11018 17 
h TOR; 1.034 13 
a 1026 HF053 10.6 
6 1033 1.067 9.0 
8 1042 1.086 TsO 
Ge 1052 TAO, 54 
13 1061 P2126 Lat 
16 1066 10 WO 
20 1072 1.149 Salt 
2k 1079 ion 3.0 


t. = 8 hours 


haplovsot ae versus log (t. + At)/At is prepared (see Figure 4-3). 


From equation (4-6) 


2 


ol x 10° (0.02)(0.88)(560)(8.5) 


G65 10001) /2.303 


kh = 


kh = 702 md-ft, 
The agreement between the permeability-pay product as calculated from a 


drawdown and a build-up test is considered reasonable, 


Example eae 


A well which produced 1,695 million cubic feet of gas since its 
last pressure survey is shut-in, The static pressure at the time of the 
previous survey was 2250 psia. The well was produced at a rate of 7.5 
million cubic feet per day immediately prior to being closed in, The well 
is located three-quarters of a mile from each of two neighboring wells, 
The following reservoir and gas characteristics are known, 


T = 620°R 
a 


u = 0.018 ep 


bui 


i > 


Z = Oo. 
a 
i) = 0°05 
h = 60 feet 


Determine the static formation pressure given the following 


ld-up data. (These represent actual data taken on a gas well in the 


Province of Alberta, ) 


Build-up Data 


Time of Shut-In Measured Pressure (P ig x10) 
At = hours P47 psia (psiA) 

2 2031 4125 

a 2041 4166 

a 2054 4219 

7 2061 42h8 
10 2071 4289 
15 2082 4335 
20 2089 4364 
24 2092 4376 
30 2098 4ho2 
36 2102 4418 
he ale OM 4439 
48 2110 hse 
60 2114 469 
72 2118 4486 
84 212] 4499 
96 2125 4516 
108 2128 4528 
120 2131 4541 
132 2133 4550 





Horner Method 


ta = 1095 (24) = 5424 hours 
bli 


From equation (4-6) 


y) 
x = ed x10" (0.018)(0.81)(620)(7-5) 2 7 5 na, 


(60) (245,000) /2.303 


The plot (Figure 4-4) is extrapolated to log (t 


Z 


+ At)/At 


Lpe 
130 
114 


=e le 


equals 1,0, where 


pee 4,950,000 


From equation (4-12) 


4,950,000 = (2250)* a 245,000 y(u) 


Pap OG} 
Yiu) = 0,46 
2.303 
Fron Figure a1 
u = Ola sie 
Substituting in equation (4-13) 
Pig = (2250)* - 245,000 
(0.38) (2.303) 
Pi. = 2,187 psia 


The static pressure calculated by the Horner method is 2,187 


psia, 


(bd) Muskat Method 
ASSIS TP ae 
The time corresponding to a dimensionless time of 0.08 is calcu- 
lated from equation (C-32), (In this instance, equation (C-32) has been 


adjusted to reflect the external boundary radius rather than the well bore 


radius. ) 
2 
: 4 £02038 )(9205) (1980) (0,08) Betncuee 
2.634 x 10 (7.5) (2200) 


So the Muskat plot should result in a straight line after approximately 
65 hours, Figure 4-5 is such a plot. The assumed static pressures are 
2,140 psia, 2,160 psia and 2,180 psia respectively, The plot illustrates 
how difficult it is to determine the precise pressure using this method, 
however, it does indicate that the pressure is in the order of 2,180 


psia, as was obtained by the Horner method, 
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FIGURE 4-5—MUSKAT PRESSURE BUILD-UP PLOT. 
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5. THE ESTIMATION OF FLOW BEHAVIOR 
FROM THEORY AND LIMITED DATA 


The previous sections dealt in some detail with the determination 
of the back pressure relationship for a gas well through a series of actual 
flow tests, This section presents methods of estimating the flow capacity 
from a limited amount of data, 

The problem of calculating the anticipated behavior of a gas 
well on the basis of reservoir and flow characteristics has been discussed 
by many authors (14)(19)(31)(35)(39)(46)(70) (80) (84) (85)(86)(88). Only 
the simplest theoretical approach to calculating the back pressure relation- 
ship will be discussed here, although mention will be made of some of the 
shortcomings of the theory and of several of the empirical methods which 
have been developed in an attempt to solve related problems. Included will 
be theoretical methods for estimating: 

Ce) the stabilized position of a back pressure curve from an 
isochronal multi-point flow test, 

(ii) the slope and stabilized position of a back pressure 
curve from a single flow point, and 

(iii) the slope and position of a stabilized back pressure 


curve calculated entirely from reservoir and gas properties. 


Positioning a Back Pressure Curve 
to Stabilized Conditions 


It has been demonstrated how the final flow rate of an iso- 
chronal test may be extended to position the back pressure line to reflect 
stabilized conditions. This positioning may also be accomplished through 
calculations based on the theory of gas flow presented in Appendix C, 


These calculations, however, are based on the knowledge of the interwell 


= 86:2 


characteristics of the formation and ideally should be employed only when 
guch information has been either obtained through a properly conducted 
build-up or drawdown test or can be accurately estimated. In addition, 

a knowledge of the true static formation pressure and of the time required 
to attain adequate stabilization is necessary. The detail of estimating 
this time is included in Appendix D, but in general may be taken 4s the 
lesser of 15 days or the theoretical time to stabilization as calculated 
by equation (D-16). (It is recognized that the true static formation 
pressure can be calculated from data taken over 4 time period of less than 
15 days. In fact, a maximum of only 6 days of build-up data is normally 
sufficient to calculate the reservoir pressure. ) 

The following is a step by step trial and error procedure which 
results in the positioning of a back pressure line of known slope to re- 
flect stabilized flow conditions. The method assumes that data of an 
isochronal test is available, however, it is based on only one of the 
series of flow tests. The particular flow rate is selected for its re- 
liability, and should be within the limits of laminar flow. 

1. Calculate from equation (C-32) the dimensionless time for 


the sandface, corresponding to the short term flow test. 


OEGa ex 1074 k Pit 


D a SS (C232) 


where 


k = permeability, md, 


Ee = average pressure, psia, 
t = time, hours. 
oe average viscosity, cp. 


> S gas filled porosity, fraction, 


. By 


ry = effective well bore radius, ft. 


To accomplish this, the average pressure during the short term flow 
period is determined and the porosity, average viscosity of the flowing 
gas and the effective well bore radius are estimated on the basis of the 
best information available. The well radius may simply represent the size 
of the casing, or may incorporate adjustment for well bore damage based on 
either build-up, drawdown, or drill stem test data (76)(85)(90)(95). ‘The 
effective permeability should be based on a build-up or drawdown test, or 
lacking one of these may be estimated from core analysis. 


Zo Determine the dimensionless drawdown, P,, in accordance 


+? 
with the appropriate boundary conditions, from Figure C-2, 
30 Calculate the dimensionless flow rate, Qs for the short 


term flow period, from equation (C-34), 





1 ar z ay 

p a (C=34) 

t qp p 2 

f 
where 

P, S dimensionless drawdown, 

or = dimensionless flow rate. 
P ~ formation pressure, psia, 
V4 = flowing sandface pressure, psia,. 


In the definition of the dimensionless flow rate by equation (C-33) 


6 
Lote se 0 wie, ,8 


Q = 
D on Pp.” (c-33) 


where 


average viscosity, cp. 


5 = 
il 


OO. te 


Ze = average compressibility factor. 
qT. = average temperature, oe 
Q = gas flow rate, millions of cubic feet 


per day at 14.65 psia and 60°F. 
k = permeability, md. 
h Spay tts 


Pe = formation pressure, psia. 


it is obvious that only the viscosity and compressibility factor will change 
with an increase in the duration of the flow test. Moreover, the variation 
in these gas properties will be small if the change in the pressure is not 


too great, so the value of, dy)» may be considered constant. 


4, Calculate the theoretical time to attain an acceptable de- 
gree of stabilization corresponding to a radius of drainage of one-half 


mile, from equation (D-16). 


KP, 
where 
te = time to stabilization, hours. 
tes exterior boundary radius, ft. 
LB arg average viscosity, cp. 
> = gas filled porosity, fraction, 


k = permeability, md. 


70) 
ut 


formation pressure, psia, 


3 Assume a stabilized pressure, Po? calculate Po and from 
equation (C-32), calculate the dimensionless time corresponding to the 


lesser of the time clculated on item 4 or 15 days, 


460... 


6. From Figure C-2, determine a P, corresponding to the time 


t 
of stabilization. 

7. From equation (C-34) and the knowledge that the dimension- 
less flow rate is constant, determine the sandface pressure. Check this 
against the value assumed in step 5, and if the difference is significant, 
repeat steps 5 to 7. 

The stabilized value of Pe - ES is then plotted against the 
constant flow rate, and a line having the slope of the isochronal line 
drawn through the point. The stabilized coefficient 'C' may readily be 
determined from the line. 

The major drawback to this approach, is the difficulty experienced 
in accurately estimating certain reservoir properties, such as the effec- 
tive well bore radius and the effective permeability. Another problem 
is related to turbulence, If the slope of the isochronal line is other 
than unity and thus reflects turbulence, it is logical to expect that an 
extension of the flow period to the time of stabilization would reduce 
the turbulence and therefore the slope of the line. This problem may be 
accounted for theoretically, but, is not serious under normal conditions 
and the method outlined should give reasonable results, particularly if 
applied to a short test where the flow in the reservoir is mainly laminar. 

Carter, Miller and Riley (14), Cornell (19), Haymaker, Binckley 
and Burgess (39), Poettmann and Schilson (70), Swift and Kiel (84), Smith 
(80), Tek, Grove and Poettmann (85), and Winestock and Colpitts (96), have 
published methods for determining the stabilized coefficient on the basis 
of a short time flow test. Each of these methods is in general, based on 
the flow theory previously presented, Most methods also incorporate em- 


pirical relationships which have been developed to overcome the complexity 


associated with a purely theoretical approach. 


BO Gas 


Establishing a Back Pressure Curve 


From @ Single Flow Point 


Theoretically, the slope of a stabilized back pressure line 
may be calculated from the combined steady-state laminar-turbulent flow 
equation. This equation relates the total pressure drop associated with 
flow, to the flow rate and to certain properties of the reservoir and the 
flowing gas. According to this theory, which assumes complete homogeneity 
throughout the reservoir, the slope of 4 stabilized back pressure curve 
should be unity under laminar flow conditions and any departure from this 
slope should be ascribed to turbulence. (Turbulence as here used includes 
all deviations from strict laminar flow and specifically those due to eddy 
currents, kinetic energy degradation, as discussed by Houpeurt (41), 
expansion and contraction losses, and the like. The turbulence is con- 
centrated in a zone near the well bore where flow velocities are greatest.) 
Experience, however, indicates that most gas wells, even when tested under 
completely atabilized conditions, yield a back pressure line of steeper 
‘slope than may be explained by using the laminar-turbulent flow equation 
which includes the effects of turbulence. It seems clear that this 
additional deviation from steady-state theory may only be explained through 
serious deviations from the concept of complete uniformity of reservoir 
properties caused by variations in permeability, presence of channels and 
fractures, partial impregnation of the formation with drilling fluids, and 
variable saturation with respect to liquids. 

Notwithstanding the above mentioned problems, it may at times be 
necessary to estimate, by calculation, the slope of the back pressure re- 
lationship. In this situation, if it igs assumed that the formation is 
reasonably uniform throughout, and that the effective well bore is approxi- 
mately circular, equation (C-15) may be ae to calculate the slope of the 


back pressure curve. 


se 


6 r 
ip) 
Pp? =a 1.418 x 10 02,8 (=) . 
kh s 
(C#15) 


-6 
rae cere ae c = 
aa oe a 


h 
where 
Pe = formation pressure, psia. 
P = flowing sandface pressure, psia. 


k = permeability, md, 


h = pay, ft. 
oT On average viscosity, cp. 
da = average temperature, aa 


Z = average compressibility factor. 


8 

ae exterior boundary radius, ft. 
ry = effective well bore radius, ft. 
B = turbulence factor, 


Q = gas flow rate, millions of cubic feet 
per day at 14.65 psi and 60°F, 
Govier (35) has presented a convenient means of using this equation 
through the introduction of the turbulence contribution factor. As dis- 
cussed in Appendix C, the laminar-turbulent radial flow equation may be ex- 


pressed in the following form: 


y - P . ea Q Fy 
a (c-22) 


the turbulence contribution factor, is described by the following 


F, haope —= GQ i (c-23) 
y_k Pv oll 
a 8 


ln (r,/r,) 


where Fis 


equation. 


eee 


It should be noted at this time, that the above equations have been 
developed on the basis of steady-state theory and that flow in a reservoir 
is an unsteady-state phenomenon. Nonetheless, use of steady-state equations 
4s satisfactory provided it is assumed that flow originates at the apparent 
or steady-state radius of drainage. Janicek and Katz (45) have suggested 
a definite ratio of the steady-state radius of drainage and the effective 
radius of drainage, (0.749 when the effective radius js moving and 0.606 
when it is stationary). 

The following is the recommended procedure for calculating a 
stabilized back pressure curve from a single flow test. 

at Plot the difference of the static and the flowing pressures 
squared, (P, - Po), versus the flow rate, Q, on log coordinates. 

ae Position the point to reflect stabilized flow conditions in 
the reservoir by the method presented earlier in this section, or by some 
other acceptable approach, 

sh Calculate the turbulence contribution factor corresponding 
to the flow rate, Q, from equation (c-23), and correct the stabilized point 
as determined under item 2, to reflect laminar flow conditions. The perme- 
ability used in this calculation should be the effective permeability as 
determined from a build-up or drawdown test, and the effective well bore 
radius may be adjusted to incorporate any damage or stimulation that may have 
occurred during drilling or subsequent operations at the well. 

h, Draw a line representing laminar flow conditions (of slope 
equal to unity) through the stabilized laminar flow point. 

Bh Calculate the turbulence contribution factor as a function 
of Q, from equation (C-23). 

6. Select a series of flow rates representing the range over 


which the well will normally produce and determine the value of (P . «- f a) 
f 8 


for laminar flow conditions from the line of slope equal to unity for each 
of the selected flow rates. 

te Calculate the turbulence contribution factor for each of 
the flow rates, 

8. Correct the previously determined values of 


FE) 


2 
(Ba 4 ms 


ramivac for turbulence by multiplication with the appropriate 
turbulence contribution factors calculated in step 7. 

9. Plot the adjusted values of [pe -~ ey on log coordinates 
against the selected values of Q. The best fit straight line is drawn 
through these points. This line should represent a stabilized back pressure 
curve, 

Several authorities (14)(19)(31)(80) have presented methods for 
calculating the slope of a back pressure curve from a single flow test, 

Most of these methods incorporate some form of empirical relationship 

which attempts to account for the extreme effect that the immediate area 

of the well bore has on the slope of the curve. One example of such a 
method is that developed by Smith (80). The method accounts for the devi- 
ation of gas well behavior from laminar flow, by combining the theory of 
flow with a procedure for measuring and calculating well bore effects for 
individual wells. Smith concluded that the Slope of a back pressure curve 
can be described by the equation for the afiie? transient flow of gas in 

a radial system, modified to include a completion or "skin" factor and a 
"rate of flow" factor, The rate of flow factor accounts for the energy loss 
associated with the turbulence in the region of the well bore, and along with 
the skin factor is determined from a pressure build-up test. In cases 

where the slope of a back pressure curve cannot be directly measured, and 
where it is impossible to accurately estimate the effective permeability 

or well radius, and where serious deviations from homogeneity are expected, 


a method such as that proposed by Smith, is recommended for estimating the 


= Ole e 


slope» 


Establishing a Back Pressure Curve from 


Gas Properties and Reservoir Characteristics 


Study of equation (C-15) results in the conclusion that theo- 
retically the pressure drop associated with any chosen flow rate can be 
calculated for a given reservoir, provided certain properties of the flow- 
ing gas and the reservoir can be estimated. Tek, Coats, and Katz (88) have 
applied this approach to develop a graphical representation of the theo- 
retical back pressure relationship for a gas of certain gravity, 98 a 
function of the flow rate and the pay, as well as the permeability and the 
discovery pressure of the reservoir. This graphical repre epterionnle pre= 
sented as Figure 5-2 for 4 0.6 gravity gas and is useful for estimating the 
potential of a well to produce when no flow test data is available. The 
problems and limitations associated with this theoretical calculation are 
essentially the same as those encountered in calculating the slope from a 
single point test. It should be emphasized, that this theoretical compu- 
tation assumes reservoir homogeneity and a circular well bore. 

The following is a recommended procedure for calculating a back 
pressure curve by this method, 

bs Estimate the average viscosity and compressibility factor 
on the basis of the static reservoir pressure, the approximate drawdown 
at which the well is normally produced, the gas gravity and the bottom 
hole temperature. 

Co Estimate the effective permeability, pay thickness, well 
bore radius and radius of drainage from core analyses, logs, build-up or 
drawdown data if available, and from the spacing of wells in the area. 


Bh Calculate the turbulence factor, 8, from equation (C=16iie 


B = an a aka (C#16) 
free 
where 
k = permeability, md. 
4, Assume a reasonable flowing sandface pressure, P.. 
do From equation (D-16) and an average pressure, calculate 


the theoretical time to adequate stabilization, The average pressure 


may be taken as 


P +P 
peek 8. 
2. 





in other cases, 


6. Calculate the dimensionless time corresponding to the lesser 
of the time to stabilization or 15 days, from equation (C-32) and the 
assumed sandface pressure, 

(e Determine from Figure C-2, a dimensionless drawdown, Pie 
corresponding to the time of stabilization, 

8. Calculate the dimensionless flow rate, Qn» from equation 
(C-34). 

96 Calculate the flow rate, Q, corresponding to 'an' from 
equation (C-33), 

10. Calculate the sandface pressure, Pos from equation (C-15), 
If this pressure differs significantly from the assumed pressure, steps 
4 to 10 are repeated, using the calculated P in step 10 as the assumed 
He in step 4, (Steps 4 to 9 inclusive, calculate a stabilized flow rate 


corresponding to a reasonable sandface pressure, These steps may be omit- 


ted, a stabilized flow rate assumed, and a sandface pressure calculated as 


OO ae 


in step 10, If the calculated pressure is not a reasonable one, another 
flow rate is assumed and step 10 is repeated until the resulting sandface 
pressure is reasonable. ) 

11, Plot on log coordinates, the resulting (P, - Po”) against 
the 'Q' calculated in step 9. Determine the turbulence contribution 
factor, Fy at this flow rate, and adjust the point to reflect only laminar 
flow. A line representing laminar flow conditions (slope equal to unity ) 
is then drawn through the point. 

IGE, Adjust the line for turbulence by the method previously pre=- 
sented. The resulting curve should represent the capability of the well to 
produce under stabilized conditions. 

The following example illustrates the recommended procedures for 
determining the proper slope of a back pressure curve and for positioning 
the stabilized curve. It also illustrates the mechanics of the calculations 
involved in constructing a stabilized back pressure curve from gas and 


reservoir properties. 


Example oat 


A well is flowed at a rate of 5.0 MMcfd for three hours. The 
sandface pressure after this time is 2890 psia. The well is drilled on 
one section spacing and the following reservoir and flow characteristics 


are known, 


ad 
it 


3030 psia (from previous pressure build-up). 
G mt dO rad 
u = 0,02 (at estimated average reservoir 
conditions). 


- = 14 md (from previous pressure build-up). 


OT = 


h = Ws ft. 

7 fe) 

Pe = 630 -R, 

Zo = 0,81 (at estimated average reservoir 


conditions), 
Construct a stabilized back pressure plot of correct slope. 
Solution 


At Q= 5.0 MMcfa 
Pe - Po“ = 829,000 


This point is plotted in the conventional manner (see Figure 5-1), 
The time to which the line is positioned to reflect adequate 


stabilization may be calculated from equation (D-16). 


2 2 
te =eT00” f Hee = 1009 (2640) °(0.02)(0.07) | n54 youre 
k Fe (14) (3030) 


The dimensionless time corresponding to three hours is determined from 
equation (C-32). 
al 
tas 2.634 x 10 14) (2960) (3 = 260,000 
2 
(0,02)(0.07)(0.3) 
From Figure C-2 


Fe = 6,60 


and from equation (C-3}) 


2 2 
a, (2 a Sama a, 


3030° 6.60 


Assume a sandface pressure of 2800 psia after 230 hours. 


P = 3030_+ 2800 = 2915 psia 


a ve! 


oo leh 


From equation (C-32) 


ty 2.634 x 1074 14) (2915) (230) _ 1.96 x 10" 
(0.02)(0.07)(0.3)* 


From Figure Ce2 


igh pang Bono 


Substituting in equation (C-34) 


BOSC eee 


3030° 





= 207013700075) 


Pe = 2843 psia 


So after 230 hours, and at a flow rate of 5 MMcfd, PA - Bag = 1,100,000. 


The apparent radius of drainage may be determined from the 
relationship developed by Janicek and Katz (45). 


ees OE TUO(PONO le =elO7 Trt. 


me is then evaluated from equation (C-23). 


F, a = 1+ 0.0124 Q 
CO OS) (LN) mre (O93 G45) 


1n(1977/0.3) 
At 5 MMcfd, the turbulence contribution factor is 


re = 1 + 90,0124 (5.0) = 1.062 


2 a : 
So the value of Pe - Ps » under stabilized laminar flow conditions would 


be 


p 2 _ p 2 = 242004000. 1 436 009 


i iS 1.062 


This point is plotted on Figure 5-1, and a line of slope unity is drawn 


through it. 


2 2 : 
The value of, Pe - Ps » for various flow rates under laminar 


20 


conditions ig then multiplied by the turbulence contribution factor to 
correct for turbulence. If we assume that the well in question will 
normally be produced over the range 2 to 20 MMcfd, the turbulence core 


rection is calculated accordingly. 


iy : P.*) Share (p* = By} “aioe g 
Q-MMcfd fs (fron laminar plot) (adjusted) 
2 1.025 O41 0.42 
5 1.062 1.04 1,10 
10 EA. Zul 2.37 
15 1.186 3.20 3.80 
20 1.2h8 4.25 5.30 


2 
f 


best fit straight line is drawn through the points. The resulting value 


The adjusted values of P,” - Bf are then plotted on Figure 5-1] and the 
of 'n' is 0.93. This line represents the stabilized back pressure line, 
Had the flow point of 5 MMcfd at a pressure of 2890 psia not been 
known, it could have been estimated, or the stabilized point could be esti- 
mated as follows, 
The turbulence factor, 8, is calculated from equation (C-16), 


10 

= Kellie 10- 12Ps29xq 10° 
(14) 4/3 

=) 


The value of ina ~ F. is then calculated for the condition 


of stabilization, from equation (C-15), 


6 
P 2 LP est 1.418 x 10 (0.02) (630)(0.81) (5.0) 1n(1977/0. 3) ; 


; (14)(45) 
3.14 x is (12.2 x 10°) (0.78) (0.81) (630) (52) ely 1 
(45°) 0.3 1977 


= 1,073,000 


- 100 - 


Knowing the formation pressure is 3030 psia, 


eared = 2847 psia 


This compares with the stabilized sandface pressure of 2843 psia which was 
calculated from the pressure measured at three hours. 

All of the methods presented or referred to in this section are 
approximations, consequently the resulting estimation of flow behavior is 
of limited use. A reasonable application of these methods, however, could 
be the interpretation of single point tests conducted periodically through= 
out the life of a well, provided, they are spaced appropriately with regard 
to multi-point tests. Another important use might be interpreting tests 
where the results of 4 multi-point test appear to be erroneous. In these 
situations, calculations may lead to a more representative back pressure 
curve than that resulting from a series of flow rates. In any case, due 
largely to lack of reservoir homogeneity, theory of flow in a reservoir has 
not yet advanced to the point where methods of calculation can replace 


actual flow tests. 
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Q, millions of standard cubic feet per day 


FIGURE 5-|— CONSTRUCTION OF A STABILIZED BACK PRESSURE CURVE 
FROM A SINGLE FLOW POINT. 
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FIGURE 5-2-PERFORMANCE CURVES FROM CORE DATA FOR 0.6 GRAVITY GAS 
(Reproduced from Tek Coats and Katz (88).) 
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APPENDIX A 


NATURAL GAS PROPERTIES 


This appendix is intended to familiarize the reader with those 
physical properties of natural gas that are of importance in the evaluation 
of the performance of gas wells, It also presents a summary of the formulae 
and commonly used graphs and tables related to these properties, 

The properties of a natural gas may be determined either directly 
from laboratory tests or by calculation from the known chemical composition 
of the gas. In the latter case the calculations are based on the proper- 
ties of the individual component of the gas and upon physical laws relating 


the properties of the components to those of the gas mixture, 


Properties of Constituent Components 


The relevant physical properties of the important constituent 
components of natural gas are listed in Table A-l, This table presents 
the recognized values of the critical temperature, critical pressure, 
molecular weight, specific gravity, vaporizing volume ratio and heating 
value of these components, 

The critical pressure and the critical temperature are the 
pressure and temperature under which the gaseous and liquid states be- 
come indistinguishable. For a single component system, the critical 
temperature is the temperature above which a gas cannot be liquified by 
the application of pressure alone. The critical pressure is the vapor 
pressure at the critical temperature, 

The molecular weight is the sum of the atomic weights of the 
atoms in a molecule of the substance, the weight of an atom of oxygen 


being taken as 16. 


= ele. 


The specific gravity of a gaseous substance is the ratio ofthe 
weight of a volume of the gas to that of the same volume of pure ary aiz., 
both at a standard pressure and temperature. The specific gravity may be 
determined by direct measurement, but, from the above definition, it may 
readily be demonstrated that the specific gravity of an ideal gas is the 
ratio of the molecular weight of the gas to the molecular weight of pure 
dry air which is 28.966. For practical purposes, it may be considered 
that this relationship also applies to real gases. 

The vaporizing volume ratio is the ratio of a volume of gas at 
@ standard pressure and temperature to the volume that it would occupy in 
the liquid state at the same pressure and temperature, 

The gross (or higher) heating value of a gas is the total heat 
liberated when a unit volume of the gas is burned under specified con- 
ditions. The net (or lower) heating value is the total heat liberated, 
less the latent heat of vaporization of the water formed in the combustion, 
This property is usually expressed in British Thermal Units of heat per 


cubic foot of gas at standard conditions. 


The Gas Law 


In the case of gases a relationship between pressure, tempera— 
ture and specific volume exists. This relationship is commonly known as 


the Gas Law and may be written as follows. 


ZRT 


ae Gah oak art eae 
where 
P = pressure under which the gas exists, 
psia, 
Vv = specific volume of the @a8s, cubic 


feet per pound, 


Pty ae 


compressibility factor for the gas. 
absolute temperature, degrees Rankine 
(degrees Rankine = degrees Fahrenheit 
+ 460), 

Gas Law constant, 1545 foot pounds 
per pound mole degrees Rankine, 
molecular weight of the gas, pounds 


per pound mole, 


Equation (A-1) can be rearranged in a more convenient form by 


introducing the relationship 


where 


then 


where 


< [he 


28,966 G pees 


specific gravity of the gas. 


144 GP 
53.34 TZ (A=3) 


density of the gas, lbs. - mass per 


cubic foot, 


Compressibility of Natural Gases 


The compressibility factor, Z, is defined by equation (A-1) 


and in fact, is the factor which validates the equation, It is best 


determined experimentally for any gas and is found to be a function of 


temperature, pressure and composition. 


According to the theorem of corresponding states, two sub- 


stances should have similar properties when at "corresponding states", 


The corresponding state is taken as the state existing at corresponding 
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fractions of the critical temperature and the critical pressure. ‘The ratio 
of the actual temperature to the critical temperature is called the reduced 
temperature and the ratio of the actual pressure to the critical pressure 
is called the reduced pressure. Thus, one would expect that at the same 
reduced temperature and pressure, the compressibility factors of two hydro-= 
carbons would be the same, 

For mixtures such as natural gas the molal average critical tempera- 
ture and pressure may be used in place of the true critical temperature and 
critical pressure for the purpose of establishing corresponding states. 

These molal average properties are called the pseudo critical temperature and 
the pseudo critical pressure and may be calculated from the composition of 
the gas as illustrated in Examples A-1 and A-2, 

Matthews, et al (57) found a reasonably systematic relationship 
between the pseudo critical properties and the gas gravity for a wide range 
of gases containing less than about three per cent non=-hydrocarbons, Further 
work (10) has led to the development of the correlation presented in Figure 
AY. The commonly used upper curve of this figure is reproduced in tabular 
form in the Interstate Oil Compact Commission (103) and the Kansas State 
Corporation Commission (104) test manuals. Figure A-1 is useful for esti- 
mating the pseudo critical properties in cases where the gas analysis is not 
available. 

Standing and Katz (82) have developed the most. widely accepted 
correlation of the compressibility factor in terms of pseudo reduced temp= 
erature and pseudo reduced pressure for natural gases. This correlation is 
presented in graphical form in Figure A-2 and appears in tabular form in the 
previously mentioned back pressure testing manuals, The Standing and Katz 
correlation is restricted to natural gases containing small amounts (say 


less than two per cent in total) of carbon dioxide and hydrogen sulphide, 


BiOo: - 


Recent work at the University of Alberta (73) sponsored by the 
Alberta Oil and Gas Conservation Board, has resulted in a method for 
estimating the effect of the acid gas constituents on the compressibility 
factor, ‘Ihe method involves the use of the compressibility factor cor- 
rection ratios presented in Figures A-3 to A=8, inclusive, The procedure 
is as follows; 

15 Calculate the pseudo reduced properties of the whole gas, 

an Determine the approximate compressibility factor from 
Figure A-2, 

Sis Determine the compressibility factor correction ratio from 
Figures A-3 to A-8, 

h, Calculate the final compressibility factor by dividing the 
quantity determined in item 2 by that determined in item Bz 

The compressibility factors resulting from the correlations in 
Figures A-3 to A-8 have been tested against actual experimental measure- 
ments for a number of sour natural gases, and are believed reliable to 
within one per cent for sour gases containing less than five mol per cent 
pentanes plus, 

Examples A-1 and A-2 illustrate the use of the above mentioned 
methods for determining the compressibility factor for a sweet and a sour 


natural gas, 


Example Aew-1 


(a) Compute the compressibility factor at a pressure of 1500 
psia and a temperature of 100°F of a sweet natural gas having the com- 
position shown below. 

For convenience, the following symbols will be used in the 


solution of this problem and, in fact, throughout the manual, 


= LOY = 


Po = critical pressure, psia. 
fe) 
a = critical temperature, RK. 
a = reduced pressure or the ratio of the 


actual pressure to critical pressure. 
iD = reduced temperature or the ratio of 
the actual temperature to the critical 
temperature. 
The above terminology will apply whether the property referred 
to is a true property for 4 pure component of a natural gas or a pseudo 


property for a mixture of these components. 


Solution 
The pseudo critical properties of the natural gas are calcu- 
lated by summing the products of the mo] *yaction times the critical 


properties for each component. 


























MW. CRITICAL CRITICAL MOL MOL 
MOL. x MOL | PRESSURE TEMPERATURE FRACTION FRACTION 
COMPONENT WI. MOL % FRACTION P. I. Kor valli 
N,, 28.016 2425 0.6304 492.3 226.9 11.08 Sys 
CH, 16.042 86.89 13.9389 673.1 343.2 584.86 298.21 
CoH 30.068 6.20 1,862 708.3 549.8 43.91 34.09 
Caile ily Gis (eee.5on 9 251288 617 22 666.0 Wetete¢) Ui. 05 
ACH 58.120 O.49 0.2848 529.0 Tatet 2.59 3.60 
nC, Hi SieloOmOcy 4 0.4301 55 aw 765.3 4.08 5.66 
iC Mio Toeow  aOre4 es 0.1132 483.5 829.7 Tee 1.99 
nC), TovinG a0. 1un | 0.1299 489.4 845.6 0.88 le 
*C (Ht 100.198 0.45 0.4509 396.8 972.3 ee pacts) 
19.0312 Peeudo P_ Pseudo I, 





= 666 = 372 
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* There is a problem in knowing what critical properties should be used 
for the residual mixture of components. Kay (50) devised a method for 
obtaining the pseudo critical properties of liquids through a relation- 
ship to the molecular weight and specific gravity of the mixture. Howe 
ever, it is convenient (and normally quite accurate) to use the critical 
properties of the next higher component, In this example the properties 
of heptanes have been used for the hexanes plus fraction. 


R 3 156059 
_ 372 ube 


From Figure A-l, at Be = 2,25 and T.. = 1.51, the compressi- 
bility factor, Z = 0.812. 

(ob) Assuming only the specific gravity of the gas were known, 
calculate the compressibility using the Brown et al correlation of pseudo 


critical properties and specific gravity. 


Solution 


The specific gravity of the gas = 5 ae = 0,657 


From Figure A-1 for 


QQ 
} 


= 0.657 
P = 669 psia 


376°R 


woe on = 2,2) 
- aes = 1,49 


= 0.805 (Figure A=2) 


= 
it 
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The excellent agreement indicates the relative accuracy frequently 
obtained through the use of Figure A-1 for a normal sweet gas, but should 


not be taken as indicative of the accuracy that will always be attained. 
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Example Aw2 
Compute the compressibility factor at a pressure of 1500 psia 


and a temperature of 100°F of a sour natural gas with the composition 


shown below, 


Solution 
Since the gas contains greater than three per cent non-hydro- 
carbons, Figure A-1 is not suitable to estimate the pseudo critical proper- 


ties. Consequently, the molal average of these properties must be calcu- 























lated, 
M.W. CRITICAL CRITICAL MOL MOL 
MOL. x MOL PRESSURE TEMPERATURE FRACTION FRACTION 

COMPONENT WT. MOL % FRACTION EA TA XEEs tN 
HS 34.076 1630) 7h 9001 1306.5 672.4 187.87 96.69 
CO, 44,010 0.30 0.1320 102.8 SLi T cre 1.64 
N, 28.016 0.460 (0.12989 492.3 226.9 2826 1.04 
CH) 16.042 84.14 13.4977 673.1 343.2 566.35 288.77 
ChHe B0.0GUm 80059 0.1774 708.3 549.8 4,18 3.2h 
CH, 44,094 0.08 0.0353 61722 666.0 0.49 0.53 
iC) H, 4 50.120, 0.03 0.017) 529.0 Tauet 0.16 0,22 
nC,H 

410 50.1208 0.02 {050116 551.1 765.3 Ov 0.15 
Cee Trace 
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18.9004 Pseudo A Pseudo T 
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Ol = 765 = 392 


Sa SOCEE 
P*Fbp = 1496 


pyplel 
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From Figure A-2, the approximate compressibility factor is 
determined, 

Z = 0,790 

From Figures A-3 and A-h the compressibility correction ratio 


= 0.960 


Corrected Z = a roximate) = at 5 = 0,823 


Correction Factor 
It can be seen from the preceeding example that for sour natural 
gases appreciable errors are involved if only the Standing and Katz cor~ 


relation (Figure A-2) is used for calculating the compressibility factor, 


Specific Gravity of Condensate Well Effluents 


Many gas wells produce an effluent which Separates into a gas 
and a condensate at the wellhead operating pressure and temperature, 
Frequently the separated liquid results from retrograde condensation and 
the fluid in the well bore itself is completely in the gaseous phase, 
The specific gravity of such a well bore gas may be calculated, assuming 
additive volumes, from the specific gravity of the gas separated at the 
wellhead, the gas condensate ratio, and the specific gravity and the 
vaporizing volume ratio of the condensate at separator conditions. The 


equation relating these factors is 


co + 4608 G. 
8 R 
Gee Sy ee (Anh ) 


1+ VOPR, 


~ 11% - 


where 

G = specific gravity of the single phase 
gas in the well bore (Air = 1). 

G = specific gravity of the separator 
gas (Air = 1). 

G = specific gravity of the separator condensate 
(Water = 1). 

R = gas condensate ratio, cubic feet per 
barrel (14.65 psi and 60°F), 

V = condensate vaporizing volume ratio, 
cubic feet per barrel (14.65 psi and 


60°F), 


Condensate Vaporizing Volume Ratio 


The condensate vaporizing volume ratio varies with the com- 
position of the condensate which in turn is reflected by its specific 
gravity. Figure A-9 presents the relationship between the vaporizing 
volume and the API gravity for pure pariffinic hydrocarbons. (For con-= 
venience in converting API gravities to liquid specific gravities, an 
abridged conversion table is presented as Table A-2),. This relation- 
ship may be used for condensate of normal composition without appreci- 
able error, 

Equation (A-4) makes no allowance in the calculation of 
specific gravity for water vapor that may be produced with the gas, 
Experiments and calculations both indicate that in the absence of 
liquid water in the flow stream, even under conditions of maximum water 


vapor content, the error resulting from this assumption is negligible 


~ es 15"= 


The following example demonstrates the use of the equation, 


Example A=3 


Given the following information with respect to separator pro- 


ducts, calculate the specific gravity of the single phase gas in the well 


bore, 
G. = 0,70 
G. 56° API 
R, = 15,000 cu.ft./bbl. at 14.65 psi and 
60°F, 
Solution 
G = 0.755 (Table A~2) 


V = 575 cu.ft./bbl. at 14,65 psi and 
60°F. (Figure A-9) 


Substituting in equation (A-h) 


0.70 + 4608(0.755) 


s m 15,000 
m thot Bee 
15,000 


G_. = 0,898 


Viscosity 


The viscosity of a pure gas depends upon both pressure and 
temperature, and in the case of a natural gas also depends upon the 
composition of the mixture, Carr, Kobayashi, and Burrows (13) have 
developed correlations for determining the viscosity of a natural gas, 


These correlations appear as Figures A=-10, A=-11 and A-12, Figure A=10 


ze) Ney 


is a plot of the viscosity at one perrrr rene pressure as a function of the 
molecular weight and the temperature of the mixture. Figures A-1l and 
Aw=1l2 give the viscosity ratio as e function of the pseudo reduced tempera~ 
ture and the pseudo reduced pressure. The viscosity ratio is the ratio of 
the viscosity at a given pressure to the viscosity at one atmosphere 


pressure. 


SS 


Use of the above mentioned correlations give reliable results 
for gases that contain no nitrogen, carbon dioxide or hydrogen sulphide, 
Corrections for the presence of these constituents may be made through 
the use of the insert plots of Figure A~10. 

The following example illustrates the use of the above mentioned 


correlations. 


Example A-4 
Calculate the viscosity of the natural gas of Example A-2 at 


i40°F and 2000 psia. 
The following symbols and definitions are for convenience and 


will be used throughout the manual, 


u = viscosity, centipoise (at any given 
pressure and temperature condition). 

ul = viscosity, centipoise (at a pressure 
of one atmosphere and any temperature 
condition), 

u/ul = viscosity ratio as determined from 


Figures A-11 and A-12, 


- 117 - 


Solution 
From Figure A-10 for a gas with a molecular weight of 18.9 
and at a temperature of 140°F 
vl = 0.0116 
Correction added to viscosity for 0.652 gravity gas with 
14, 38% HS 
= 0,0003 (Figure A~10) 
Corrected u1 = 0.0116 + 0.0003 = 0.0119 ep 
P = 765 psia 


Cc 


4 fe) 
T, 7 393 R 


a eC.’ 
P. = ha AM Sal 


600 


TY. 392 


= 1.53 


From Figure A-il for P| = 2.61 and Tem 1.53 
w/ul = 1.37 


u = 1.37 (0.0119) = 0.0163 cp. 
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TABLE A-1 


PHYSICAL PROPERTIES OF IMPORTANT COMPONENTS OF NATURAL GAS 


S| 


Ft Gas /Imp 
p rr Gal @ 60°F Gross ligating Value, 
c c Molecular and 14,65 Btu/ft> @ 60°F ana 

Component (psia) ORY Weight psia 14.65 psia (dry) 
HS 1306.5 672.4 34.076 88.22 678 
CO,, 10fene SAT. 7) @hl.010 70.18 = 
N,, 492.3 226.9 28,016 - - 
He 33.1 9.2 4,003 - - 
CH), Ossie Gage. 6, 0he ileGe 1007 
CoH 70633 549.8 30,068 47.25 1763 
CH, 617.2 666.0 hh,ook 43.76 2510 
iC) H, 72950 T3457  §8.120 36.83 3246 
nC) H, 4 Soiee «=~ (65.3 BO.120 38.22 3254 
Butanes * oueeg 755.1) 58.120 Sf. 70 3251 
iC.H4 5 483.5 829.7 72.146 32.98 3988 © 
nC Hy 5 489.4 845.6 72.146 83.32 3996 
Pentanes** 486.2 836.9 72.146 83.12 3992 
nC,H ), 440.0 914.2 86.172 29.35 47h 
nCoH ¢ 396.8 972.3 100.198 26.16 5486 
nCgH, g 361.5 1024.9 114,22 23.57 6230 
——<——— 
CoH 345.0 1073.0 128,250 21.434 6975 
Ci oflos 320.0 1115.0 142.276 19.652 T7720 


NN LL SSA ER SSS ie SSS si se tsisssitemenestenstunctnttiestaenesstoeetatonase 


* Butanes = 1/3 iC), and 2/3 nC), 


** Pentanes = 55% ic, and 45% nC, 


- Mid a 





BAO SAmUTAR 40 QEWAMOSMOD TRATHOUMT tC BATINTIORT ART Oras 


ativan 2% 


ona 0 £ eal pets hg 8 a oS te 

aya ga.88 le asta’. 7 .90e sr 

- BL.oY Gro.ui  7.Tee 8.870! Ro) 

- 01040 ,o8S5 ,ooe gf 

iil eee ee 
* - FOOsa »t ‘ $i 

yoo. $0.47 oso,af- 5.cee 2-679 jhe 

EATS en.Te 30, OF ph Ff .40) = H.9 

Giz sted do, ud ; ate 

| catenin TO a EE 
dace ‘$B. ee GSi tM Patt o.0S¢ Ay tye 

deS£ 65 8 ORE 2 aay £, [ke oho 

tes PVE eu.be- 1 TT.) “roasted 








, BORE 8Q.SE oifeer y.es6- .¢,£0e gpa 7 
3QCE : SE EE sar eSy ae », Cw eT 
SROE SE,££ @L.SY io oe SOs at **sonatast 


ir 


Syteod die Oi odd 

“Berson ESTO” SRO 

ee eso, 219 
wa be i 


a 
on 





- lel - 


TABLE A=2 


CONVERSION OF API GRAVITY TO LIQUID SPECIFIC GRAVITY AT 60°F, 


WOON AU FWNF- CO 


Specific Gravity = pa) vhs 


Specific 


Gravity 


1.076 
1.068 
1.060 
1.052 
1.044 
1.037 
1,029 
1,022 
1,014 
1.007 
1,000 
09930 
9861 
09792 
9725 
9659 
09593 
09529 
9465 
9402 
09340 
9279 
09218 
09159 
+9100 
902 
8984 
8927 
» 8871 
» 8816 
- 8762 
» 8708 
» 8654 
» 8602 


Oo 


ASiehtke APT 


Specific 
Gravity 


» 8550 
8498 
8448 
» 8398 
° 8348 
«8299 
28251 
- 8203 
8155 
. 8109 
8063 
» 8017 
7972 
01927 
- 7883 
«71839 
cules 
wage 
- 1628 
ey ble. i 
© TOUT 
«7507 
o TH67 
7428 
«7389 
Piva jep 
<(oL5 
f2q5 
$7201 
$TLbS 
~7128 


Specific 
Gravity 


» 1093 
«1057 
at0ee 
- 6988 
~6953 
°6919 
> 6886 
26852 
© 6819 
06787 
6754 
6722 
>6690 
- 6659 
26628 
6597 
- 6566 
20536 
26506 
6476 
o OLLE 
Vetta ley 4 
» 6388 
» 6360 
06331 
- 6303 
06275 
6247 
6220 
26193 
26166 
~6139 
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FIGURE A-l|— PSEUDO CRITICAL PROPERTIES OF MISCELLANEOUS NATURAL GASES. 
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(Redrawn from Brown etal, 10) 
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PSEUDO REDUCED PRESSURE, P, 
FIGURE A-2 - THE COMPRESSIBILITY FACTOR FOR NATURAL GASES. 


(Redrawn from Standing and Katz,82. Courtesy AIME.) 
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REPRODUCED FROM ROBINSON et al,(73) 
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CONDENSATE VAPORIZATION VOLUME,CUBIC FEET PER BARREL 
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FIGURE A-9 — RELATIONSHIP 
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FIGURE A-I1 - VISCOSITY RATIO VERSUS PSEUDO REDUCED TEMPERATURE. 


(Reproduced from Carr, Kobayashi and Burrows, !3. Courtesy AIME.) 
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PSEUDO REDUCED PRESSURE , P, 


FIGURE A-I2 —VISCOSITY RATIO VERSUS PSEUDO REDUCED PRESSURE. 
(Reproduced from Carr, Kobayashi and Burrows, !3. Courtesy AIME.) 
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APPENDIX B 


CALCULATION OF BOTTOM HOLE PRESSURES IN GAS WELLS 
LS LLY I LE AS AS SSE Recta ee atnesanaletioatemase ee aa 


In order to evaluate the performance of a gas well it is neces- 
Bary to know the pressure at the sandfacc of the formation under both flow- 
ing and atatic conditions. Recognizing that it is not always practical to 
obtain these pressures by direct measurement, it becomes necessary to cal- 
culate bottom hole pressures from wellhead observations, 

This appendix introduces the basic theory supporting the calcue 
lation of bottom hole pressures from wellhead pressures in the case of a 
ges phase only in the well bore, and also discusses briefly different methods 
of calculating these pressures, It discusses the main advantages and dis- 
advantages of each of the methods, and sets out a recommended procedure of 
calculation illustrated by appropriate examples. 

Also included in the appendix is a brief discussion of some of the 
methods for calculating flowing bottom hole pressures from top hole measure- 


ments for wells with both gas and liquid phases in the well bore, 


A. ONE PHASE (GAS ) PRESENT IN THE WELL BORE 


Calculation of Static Bottom 
Hole Pressures 


The static bottom hole pressure in a gas well is simply the well- 
head pressure plus the pressure due to the weight of the column of gas in 
the well bore, The starting point for the calculation of this pressure is 


the familiar mechanical energy balance equation, 





Ah aP ee ey (B-1) 


p & 


Aly Fe RS 


where 
L =» height above an arbitrary datum 
elevation, ft. 
p = density of gas, lbs - mass per 
cubic feet, 
V o velocity of flow, feet per second. 


8 = dimension conversion factor 


a 32, 17idoamass ) ft. 


(1b. force) ace 
ar = energy losses due to irreversi- 
bilities (friction). 

P = pressure, pela. 
The terms in equation (B-1) represent respectively the pressure-volume 
energy, the potential energy of position, the kinetic energy and the 
energy losses due to irreversibilities, 

For static conditions the kinetic energy and energy losses due 


to irreversibilities are zero and equation (B-1) may be written as 


dP = - ofr ab aan 


The density of the gas may be expressed in terms of tempera- 
ture, pressure, specific gravity and compressibility factor through the 
Ges Law, as expressed in equation (A-3), and when substituted into equation 
(B-2) yields 


dP G dL (B-3) 


ey >) SRR 

where 
G = specific gravity of the gas. 
aK = temperature, °R. 


Z 2 compressibility factor of the gas, 
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The only assumption involved in this equation is a single 
phase gas described by the Gas Law. Equation (B-3) is the basis of all 
methods for calculating static bottom hole pressures from top hole values. 
The equation is in differential form and must be integrated, The tempera- 
ture and compressibility of the gas in o well bore varies with the depth, 
which means that to integrate equation (3-3) a knowledge of the nature 
of these variations or some assumption regarding them must be made. In 
addition to this, the gravity of the gee may very slightly with depth, 
Since knowledge regarding this variation is seldom available, in the re- 
mainder of this appendix it will be assumed the Gravity is constant. A 
number of the methods for integreting this equation, slong with the 
necessary assumptions (assuming only a @&8 phase in the well bore) are 


discussed in the following. 


The Average Temperature and 
Gompressibility Method 


This method includes the well-known Monograph 7 method (72), The 
Railroad Commission of Texas method (106), and the Vitter method (93). The 
approach ig based on the assumption that both tne temperature and compres~ 
sibility factor are constant at some average values, a and Za? over the 
depth of the well. 

The average temperature, The ig usually taken as the arithmetic 
average of the reservoir and the wellhead temperature. A refinement which 
has been proposed (49) is the use of a logarithmetic mean temperature in- 
stead of the arithmetic mean, If the compressibility factor were truly 
constant, and if the temperature profile were linear with depth, the log- 
arithmetic mean temperature would be correct and its use would yield an 


exact answer. On the other hand, since the compressibility factor does 


oi 


vary and depends on temperature, there is no assurance that the use of a log 
mean temperature rather than an arithmetic mean would give a better answer. 

There are several ways of calculating an average compressibility 
factor, but the most common is to use the compressibility factor at the 
arithmetic average pressure and temperature. This is not necessarily better 
or worse than using the arithmetic average of the compressibility factors 
at bottom hole and at top hole conditions. 

The average temperature and compressibility method involves a 
trial and error solution, and may be applied in the form of a one-step 
calculation from the wellhead to the sandface or a multi-step (usually two) 
modification, This two step modification involves dividing the well into 
two equal sections and calculating the pressure at the mid-point of the 
flow string. The pressure at the sandface is then calculated from the mid-= 
point preesure. 

Some of the major advantages of this method are that it is well- 
known and widely used, and that it may be employed with a sour gas by ine 
corporating the compressibility correction ratios into the calculation, 
Disadvantages include the necessity of a trial and error calculation and 
the fact that the method is not accurate where there is a large variation 
in either the temperature or the compressibility factor. For the latter 
reason the method is not suitable for deep wells, and is therefore not 
recommended as a standard procedure for converting top ote pressures to 


bottom hole conditions for the static case, 


The Average Density Method 


This method proposed by a University of Alberta graduate class 
(see Preface), is based on the assumption that the density of the gas in 


the well bore mey be assumed constant at the arithmetic average of the 


ee 


bottom hole and top hole values, The method is not commonly used, although 
it has the advantage over the average temperature ond compressibility 
method of introducing only the one overall assumption concerning the key 
variable (density), as opposed to the two independent assumptions with re- 
8pect to temperature and compressibility, It may be applied in a one-step 
or multi-step calculation. 

The method involves a trial and error procedure. In cases where 
the fractional variation of density with depth is appreciably less than the 
fractional variation of temperature and/or compressibility factor, this 
method is probably superior to the previously discussed method, However, 
Since the density depends upon the temperature and compressibility factor, 
the method is not accurate where there ie a large variation in either 
temperature or compressibility; therefore, it is not suitable for deep wells 


or for adoption as a standard, 


The Poettmann Integral Method (68 


This method is based on the assumption that the temperature in 
the well bore may be assumed constant at some appropriate value with full 
recognition being given to the variation of the compressibility factor with 
pressure. Either the arithmetic or the logarithmetic average temperature 
may be used. For reasons previously discussed there appears to be no 
assurance that the latter would give a better answer than the former, 

The necessary integration of the basic equation has been carried 
out on a computer, The results are presented in tabular form (49)(68) and 
may be used directly to solve for a sandface pressure, However, the inte~- 
gration employs compressibility factors which have not been corrected for 
acid gas constituents, and as a result the method may be used only for 


sweet gases, Other disadvantages are that the solution requires double 


er rHom- 


interpolation of integrals within the Poettmann table and that inaccuracies 
are introduced where there is a large variation in the temperature, The 
method is not a suitable standard to be applied to all calculations of static 


bottom hole pressures. 


The Cullender and Smith Method (24) 


This method enables account to be taken of the variation of both 
temperature and the compressibility factor with depth. It involves numeri- 
cal integration for any number of increments down the flow string. The 
method of solution is a trial and error procedure, and the calculation may 
be made for a sour gas through incorporation of the appropriate correction 
ratios. Since each of the previously discussed methods must assume that 
some variable or combination of variables such as temperature, compressi- 
bility factor or density are a constant, it is concluded that this method 
is theoretically the best available for calculating static bottom hole 
pressures, 

It is recognized that for many wells (particularly those less 
than four or five thousand feet deep), any of the methods mentioned should 
result in static bottom hole pressures which are satisfactory. However, 
since the Cullender and Smith approach is the most suitable for deep wells, 
it is therefore adopted as the standard method and is discussed in greater 
detail in this manual. 


Cullender and Smith have rearranged equation (B-3) to 


eo TZ dP 


~ ALT) 


where 
L © length of flow string, ft. 
= the sandface pressure, psia, 
Ee = the wellhead pressure, psia. 


The right-hand side of the equation may be integrated numerically 
by determining the average value of (TZ2/P) for each of any number of in- 
crements in P between ar and Poe In general, for limits in pressure of 


me and re and designating 


P (B-5) 
the right-hand side of equation (B-4) may be expressed as 


e 


n 

22 ap = 1/2 [((Porn(PL)(T at IP}E+ (RJ seuRs) (Dnesle)o< 

P 1 © eel re) 2 Saeete icatinite? 
ie * ire - Rey ee + seen a) 


In the case of a two step calculation where only one intermediate 
value of pressure, that at the mid-depth, P, is considered, equation (Bed) 


May be expressed as 








P 
Ss 
GL Lys (Pitan tet Tey eat Ps uP a( rere Tea) 
33, 3h oz p dP = a Ww a W 8 a 8 a 
P 
or Ww 
0.0375 GL = sig Ae A ah + 7) + (P., - ake, + I.) (B-7) 


Equation (B-7) may be separated into two expressions, one for the upper 
half of the flow string and the other for the lower half. For the upper 
half 


Pardes 9 a 
0.0375 G 5 (Po - eR + La, (B-8) 
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and for the lower half 


g L Pry ; Bee 
0.0375 G 5 (p. 4 Podlae + re) (B-9) 


This relationship is equivalent to assuming the specific volume 
of gas constant at its arithmetic average value and is the specific volume 
counterpart to the average density method applied in two steps. While the 
method may be used with any number of steps, Cullender and Smith have demon- 
Btrated that the equivalent of four step accuracy may be obtained with a 
two step calculation and parabolic interpolation (or the use of Simpson's 


rule). This results in 


(P. - Ee 
0.0375 GL = [Sage aie - a + I,) (B~10) 


The solution of equation (B-7) or (B10) requires knowledge of 
the temperature at points in the well bore where the quantity, I, is to be 
evaluated, Except in the rare case where the temperature profile is known, 
it is necessary to either assume or to calculate one. Lesem et al (53) 
have presented a method of calculating the temperature-depth relationship, 
but, sufficient accuracy will normally be attained by assuming the relation- 
ship to be linear. 

The following procedure is recommended for the solution of 
equation (B-10), 

1. Calculate the left-hand side of equation (B~8) for the upper 
half of the flow string. 

| 2. Calculate I, for wellhead conditions from equation (B-5), 

So As6ume I. a i for the conditions at cre average well depth 
or at the mid-point of the flow string. 

4, Calculate ie from equation (B-8), 


Do Using the value of Po calculated in step 4 and the 


Pe oie 


arithmetic average temperature, determine the value of I. from equation 


6, Recalculate P, from equation (B-8) and if this recalculated 
value does not ayree within one psi with the De calculated in step 4, re- 
peat the procedure with the new calculated value of i until agreement is 
reached with the previously calculated value. 

Hie Repeat the entire procedure, using equation (B-9), for the 
lower half of the well and obtain a value of the sandface pressure, 

8, Apply Simpson's rule as expressed in equation (B-10), to 
obtain a more accurate value of the sandface pressure. 

The following example illustrates the use of the Cullender and 


Smith method. 


Example B-1l (A Deep Sweet Gas Well 

A well produces 0.75 gravity gas of a known composition from an 
average formation depth of 10,000 feet. The shut-in wellhead pressure is 
2500 psia. Wellhead and formation temperatures are 35°F and 265 °F re- 
spectively. Calculate the static bottom hole pressure by the Cullender 
and Smith method, using a two step calculation and a parabolic interpo- 


lation with Simpson's rule, 


Solution 


P = 667 psia . 
(From analysis) 
408°R 


La | 
i] 


y. T 495 + 705 
= si— = 600°R 








O 
ee [2 = 1.728 


i oper 


_ 600 
oy ~ hot 1.571 

4 rc 
to Bho 123 

2500. 
aries “fay = 3. THB 
Z =# 0,593 (Figure A-2) 
WwW 


Substituting in equation (B-5) 


i a (495)(0.593) = 0.117% 


Solving the left-hand side of equations (B-8) and (B-9) 


0.0375 G Ze 0.0375 (0.75) £20,000) 140.63 


Step al (The upper half af the flow string) 


First trial 
Assume I = I « 0.1174 
a w 
Solving equation (B~8) for Po 
. UOLO3IO wt (P_ = 2500)(0.1174 + 0.1174) 


Li = 3099 psia 


Second trial 
ta eee ee : 
Pe oa 4.646 


Zo = 0.780 (Figure A=2) 


600)(0, 80 


Once again SOlving equation (B-8) for re 
140.63 = (P. = 2500)(0.1510 + 0.1174) 


ie = 3024 psia 
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Third trial] 


Pp = gree = 53h 
Be 667 
Ze = 0.775 
aR Pe (600)(0.775) = 0.1538 
(302k) 


From equation (B-8) 


140.63 = bPe - 2500)(0.1538 + 0.1174) 


un = 3019 psia 


Fourth trial 


ph td Bs 
ot = B67 4.526 
Z = 0.775 


a 


600) (0, 
Ios sore} Tie) 0.1540 


From equation (B-8) 


140.63 = ae ~ 2500)(0.1540 + 0.1174) 


Po = 3018 psia 


Step 2 (The lower half of the flow strine) 


First trial 
Assume I =I = 0,15h0 
8 a 
From equation (B-9) 


140.63 = (Ee - 3018)(0.1540 + 0.1540) 


i = 3475 psia 


(Figure A-2) 


(Figure A-2) 


- 152 - 


Second trial 


Zee 0.894 


0 0,894 
a ca 2 e 814 


From equation (B-9) 
140.63 = (P. - 3018)(0.1540 + 0.1814) 


Be = 3437 psia 


Third trial 
Z = 0,892 (Figure A-2) 


From equation (B-9) 


140.63 = kre = 3018)(0.1540 + 0.1830) 
a = 3435 peta 


Simpson's rule is now applied through use of equation (B-10) 


to improve the accuracy of the calculated sandface pressure. 


- P 
8 Ww 


3 


281.26 = 





(0.1174 + 4(0.1540) + 0.1830) 


P,- P, = 921 peia 
ce = 2500 + 921 psia = 3421 peia 
For comparative purposes, the same method has been used on the 


same well with only a one step calculation. This approach involves a 


direct assumption of bottom hole conditions rather than those at a mid- 


etHiss. 


point, and trial calculations until the required agreement is reached. The 
sandface pressure calculated this way is 3436 psia. This and other similar 
examples have indicated that a one-step calculation would not yield satisfactory 
results for deep, relatively high pressure wells, On the other hand, a series 
of calculations for shallow sweet gas wells (less than 5000 feet deep), has 
shown that a one step calculation can be used without significant loss of ac- 
curacy, 

The method may also be used to calculate the static bottom hole 
pressure for a well which produces sour gas. In this case, a correction is 
applied to the compressibility factors prior to their use in equation (B-5). 

To support the theoretically based conclusion that the Cullender 
and Smith method is the best available, the Board has carried out a series 
of calculations of bottom hole pressures for a shallow sweet gas well, an 
intermediate depth sweet gas well, a deep sweet gas well and a deep sour 
gas well. Although the results of these calculations are not completely 
conclusive, they support the use of the Cullender and Smith method as a 


standard for converting static top hole pressures to bottom hole conditions, 


Calculation of Flowing Bottom 
Hole Pressures 


The starting point for flowing bottom hole pressure calculations 


is the previously mentioned mechanical energy balance equation. 


uy S + ab + 2+ aF = 0 (B-1) 


c 
The kinetic energy term, VaV/e has been shown to be negligible 
compared to the other terms in the equation and is disregarded in further 
derivations. The energy loss term, dF, represents the mechanical energy 
degraded to thermal energy because of irreversibilities attending the flow 
and is commonly called the friction loss term, In fact, it is the term 


required to "balance" equation (B-1) and is defined by the equation, It 


2154. = 


ie convenient to express the friction loss in terms of another quantity, 


the friction factor®, f, by the well-known Fanning equation 


: | 
Fs eo (B11) 
6. 


where 
f a friction factor. 


& = dimension conversion factor 


( = ) ft, 
- 32.17 lb-mass A 


(1b. force) sec. 
D = inside diameter of pipe, ft. 
V = velocity of flow, ft. per sec. 
The actual lineal velocity of the gas at any point in the well 
bore, in feet per second, may be related to the production, Q in MMcfd at 


14.65 peia and 60°F through the equation 
a8 p”(3600)(24)(520) PV 
(4) 10° (14.65) TZ (Be12) 


which reduces to 


yw e452 T2Q 
2 
PD (Be13) 


* A similar friction factor known as the Blasius friction factor is defined 
by the equation 


26D 


This is commonly used in civil engineering references. It is to be noted 
that f, = hf, The Fanning factor is used throughout this manual. 
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Combining equations (B11) and (U-13) and the density relationsliip given 
in equation (A-3) with equation (H-1) the following is obtained, 


Ag pene 2 
2a TZ (2): wee ar (=) (o-1328 =O (pe 1b) 
is D 


which reduces to 


53.34 TZ “ai ae 2234s van (m2? o? = 0 (8-15) 
G \P gD P 


Equation (B-15) may be considered as the basic differential 


equation relating pressure to flow rate, pipe diameter, vertical position, 


gas gravity, density and compressibility. The only assumptions in this 
equation are: 

up Steady state flow. 

en Single phase gas ‘described by the Gas Law. 

3. Negligible kinetic energy cffects. 

h, No work done on or by the gas over the section under 
consideration. 

If 'T' and 'Z' are taken as constant at average valves, and in 
the normal case where 'G' is constant and where there is no choice but to 


assume 'f' to be constant, the equation may readily be integrated to give 


Q = 0.1000 (B16) 





where 


20GL 
53.34 ie 


S z 


L #2 length of the gas column, feet. 


d= inside pipe diameter, inches. 


oP boom. 


Aa will be discussed later the assumption of soa tarnichl and 
'7Z* 4g not a necessary one, but it is inherent in several of the methods 
to be discussed. For this reason equation (B-16) is sometimes considered 
the "basic" equation. As compared with equation (Be15) it incorporates 
the following further assumptions, 

Se The gas temperature is constant at some average value, 
ustially taken as the arithmetic average of the terminal values, 

6, The compressibility of the gas is constant at some average 
value, usually taken as that corresponding with the arithmetic average of 
the terminal pressures and temperatures. 

7, The gas gravity is constant (not usually a serious 
assumption). 


83 The friction factor 'f' ia constant. 


Regardless of whether equation (B-15) or equation (B-16) is 
employed, one of the critical quantities is the friction factor 'f’, 
This factor is defined by equations (B-1) and (B-11) and is the factor 
which validates these enisci sans 

Much experimental work has been carried out, particularly in 
horizontal pipes to determine the variables which influence ‘f' and in 
an effort to develop methods of predicting it. In recent years this 
work has been augmented by theoretical studies, It is now well estab- 
lished that in the region of turbulent flow, as encountered in gas wells, 
the friction factor depends upon two quantities, the relative roughness 
and the Reynolds number. 

The relative roughness, 45/d, is the ratio of the absolute 


roughness, § (the distance from the peaks to the valleys in pipe wall 


a hy ae 


irregularities), to the internal pipe diameter'd', 
The Reynolds number is @ grouping of those variables that deter- 
mine the flow pattern in geometrically similar Systems, and may be expressed 


as 


BEN Es (B+17) 


where 
os viscosity of the gas, lb-mass per ft. 
Bec. 


For steady state flow this may be expressed in terms of 'Q' to give 


yw 20011 69 (B-18) 


Re wa 


where 

G ba gas gravity. 

Q Sal flow =~ millions of cubic feet per 

day at 14.65 peia and 60°F. 

u o viscosity of the gas, centipoise, 

a « pipe diameter, inches. 

Weymouth (94) developed one of the earliest correlations for 
the friction factor 'f', He found that for high rates of flow (high 
Reynolds numbers) the friction factor depended largely on the pipe diameter 
and was inversely proportional to the 1/3 power of the diameter, 


Weymouth's equation for the friction factor is 


pw 22908 
aor 
This equation does not take account of variations in absolute roughness or 
of the effect of Reynolds number, At the high Reynolds numbers normally 


encountered in gas wells, the equations give a friction factor corresponding 


56s 


to very dirty pipe (see Figure Be2). The Weymouth friction factor is 
therefore, not suitable for general application. 

Probably the best equations now available for the prediction of 
the friction factor are those based on work by Von Karman (47) and 
Nikuradse (64), or ones closely related to them. The Colebrook (16) 


equation 








d 

4.67 = 
fam = 4,0 log $+ 2.28 - 4.0 log 1+ (B-19) 
No if. 
Re 
a 
4s such a relationship. For values of less than 0.01 the Reynolds 
pet. 


number is found to have no further influence on the friction factor and 
the equation reduces to 


fa = 4.0 log $+ 2,28 (B-20) 


These equations are awkward to use as such, but, may be represented 
graphically as in Figure B-l taken from Knudsen and Katz (51). This 
figure may be used with confidence to predict the friction factor in 
cases where the Reynolds number and the relative roughness are known, 

The Reynolds number presents little problem being readily calculated from 
equation (B-18), however, the relative roughness is not so easily deter- 


mined, 


Relative Roughness 

| For clean new pipe the relative roughness is determined by the 
method of manufacture and usually reflects on absolute roughness, 6, of 
0.00055 to 0.0019 (22) (78) (79). For new pipe or tubing used in gas wells 


the absolute roughness has been found to be in the order of 0.0006 or 
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0.00065. 

A convenient way to present the relative roughness is in the form 
of a plot of the absolute roughness versus the pipe diameter. Figure B-2 
is such a plot and includes curves for well tubing, commercial steel and 
very dirty pipe. Also included are actual data points reported by Smith 
and co-workers (78) (79), 

Figure B-2 also shows the relative roughness which corresponds 
to the use of the Weymouth friction factor at four levels of Reynolds 
numbers, 10", 10°, 10° and 10!” It is to be noted that the Weymouth 
friction factor corresponds to an unreasonably high relative roughness in 
small pipes and a high or low relative roughness in intermediate diameter 
pipes, depending upon whether the Reynolds number is high or low. 

The Figure is recommended as a guide in the estimation of the 
relative roughness. In the absence of special knowledge of the condition of 


the pipe, the line corresponding with an absolute roughness of 0.00060 


inches is recommended, 


The Average Temperature and 
Compressibility Method 


This method is similar to the Monograph 7 method, the Railroad 
Commission of Texas method or the Modified Vitter method, the only real 
difference being in the evaluation of the friction loss, The method is 
capable of incorporating any friction factor relationship. However, as 
employed in the Monograph 7 and Texas Railroad Commission methods, it in- 
corporates the unrealistic Weymouth relationship. The Modified Vitter 
method also uses the Weymouth factor, but, introduces a compensating 
correction factor. 


Each of the methods incorporate the basic assumption of constant 


temperature and compressibility, and ag a result is not accurate forrjdeep 
wells where a large variation in the temperature and compressibility factor 
may be expected. In addition, where the eavaKcr factor has been assumed 

a function of the pipe diameter only (through use of the Weymouth equation), 
this will usually introduce additional error into the calculation. If the 
friction factor used is an appropriate one (such as the Miller (60) cor- 
rected Weymouth factor used in the Modified Vitter method, or preferably 

4 friction factor calculated from the Colebrook equation), this method of 
calculation gives reasonably accurate results for relatively shallow wells. 
However, since it is not an acceptable approach for deep wells, it cannot 


be used as o standard method for calculating flowing bottom hole pressures. 


The Average Density Method 
As in the static pressure case, the flowing calculation by this 


method is based on the assumption that the density of the gas in the well 
bore may be assumed constant at the arithmetic average of the bottom hole 
and top hole values. The friction factor is based on the appropriate 
relative roughness and the Reynolds number, and is determined from Figure 
bel, 

As in the static case, this method gives accurate results for 
relatively shallow wells, but, due to the dependence of the density upon 
the temperature and compressibility, is not acceptable for deep wells or 
as G standard procedure, (It may be more accurate than the previously 
discussed method, if the fractional variation of density with depth is 


appreciably less than the fractional variation of temperature and/or com- 


pressibility), 


all 


the Sukkar and Cornell Method (63) 

This method is similar to the Poettmann Integral method as pre- 
sented for the static bottom hole pressure calculation, and is based on the 
assumption that the temperature in the well bore may be assumed constant at 
some average value, 

The basic equation (B-15) has been solved by numerical inte- 
gration and the results are presented in tabular form (49) (83). The 
friction factor as used in this method may be determined from Figure Bl, 
on the basis of an appropriate relative roughness and the Reynolds number. 
Since the evaluation of the integral has been based on the Standing-Katz 
correlation, the m-thod is not applicable for sour gases. Also, the 
method is not accurate where there is a large variation in temperature 
(deep wells), so is not acceptable as a standard procedure for calculating 


flowing sandface pressures, 


The Cullender and 
Smith Method 


As in the case of a static well, Cullender and Smith have de- 
veloped a numerical integration method for calculating flowing bottom 
hole pressures from wellhead measurements. Again, the method enables 
account to be taken of the variation in temperature and compressibility 
factor with depth. It permits a choice of friction factor to suit the 
relative roughness and the Reynolds number, and also may be used for a 
sour gas by incorporating compressibility correction ratios into the 
calculation, 

Since this method involves fewer assumptions than any of the 
other methods, and because it is applicable to a wider range of well 


conditions, it is adopted as the standard method for calculating bottom 


Se oS 


hole pressures and is considered in some detail in the following. 
By substituting for the density terms and the friction factor 
term in the basic differential equation (B-15), Cullender and Smith have 


developed the following general formula for the calculation of the flowing 


pandface pressure. 


P 
8 
P 
1000 GL Ts ae 
3h * rae EMO EAN OR. 
ye p 
2 TZ | 
eat (Be21) 
P 
WwW 
where 
2 
re g 2.6665 f (p-22) 
a? 


and all other terms are as previously defined. 

Equation (B-21) contains no assumptions or approximations other 
than those of the basic differential equation. As in the case of equation 
(Bel), the right-hand side of equation (B-21) may be integrated numerically 
for any number of increments down the flow string. Precise solution re- 
quires a detailed knowledge of the temperature |}: ofile. ae in the 
rare case where it is known, the profile is normally assumed to be linear. 

Equation (B-22) may be solved for 'F' directly by determining 
the average viscosity (at the average temperature and pressure), the 
Reynolds number and the friction factor. However, a simplification based 
on an absolute roughness of 0.0006 inches, as calculated from the 
Nikuradse friction factor equation for fully developed turbulent flow 


gives 


; 0,10 
Be Ay OnE: (B23) 


ma to a has 


for internal pipe diameters less than 4,277 inches and 


fee Ono 
OEE} en (Be2h) 


for internal pipe diameters greater than 4.277 inches. 

Values of F/Q (called F by Cullender and Smith nd calculated 
from equations (B~23) and (B-24)) for various tubing and casing sizes are 
presented in Table Bel, 

Equation (B-21) may be solved approximately by numerical inte- 
gration between the limits P, and Lm by employing a two step procedure 
using an average well pressure, Pa at half the well depth. 


The equation may then be written as 











Bs 
P 
1000 GL Toa te : (Pla Py) (L + 1) +e lieear id Cede 1) 
53. 3u 2 p 2 2 2 
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1000 (B25) 
p 
W 
or..37.5GL = (P= P.) (1 +3) + (P, =P) (tod) (B26) 
where 
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mupe geocele (B-27) 
2 
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Equation (B-26) may be compared to equation (B-7) which was 
developed for static bottom hole pressure calculations. 

A two step procedure and the use of Simpson's rule will also give 
the equivalent of four step accuracy for the calculation of the flowing 


bottom hole pressures. For the following case, the use of Simpson's rule 


~\bA- 
results in 
p seuciop 
37.5GL «# ————— Cae LA bag (B~28) 
3 


Since a two step procedure will normally be used, it is convenient to divide 
equation (B-26) into two expressions, one for the upper half and the other for 


the lower half of the flow string. For the upper half 
Lee Be 
OTs Suche (ihe maple (Levit } (Be29) 
and for the lower half 


37.5 GE (P, =P.) (I, + I,) (B-30) 


As in the static calculation the assumption that the tempera- 
ture gradient is linear in the well bore is also used to solve equation 
(B-26) or (B-28) for the flowing bottom hole pressure. 

The following general procedure is recommended for solving 
these equations. 

1. Calculate 37.5G 2, (for the upper half of the flow string). 

De Calculatelr. efron equation (Be23) or (B-2h). 

3. Determine I for wellhead conditions, from equation (B=27). 

4, Assume I. = I, for the conditions at the average well depth. 

5. Calculate P_ from equation (Be29). 

6, Using the value of re calculated in step 5 and the arith- 
metic average temperature, determine the value of I, from equation (B-27). 

1. Recalculate P| from equation (B29) and if this recalcu- 
lated value does not agree within one psi with the a calculated in step 5, 
repeat the procedure with the new calculated value of Pe until agreement 
is reached with the previously calculated value. 

8. Repeat the entire procedure using equation (B-30) for the 
lower half of the well and obtain a value of the sandface pressure, 


9. Apply Simpson's rule as expressed in equation (B28) to 


16st 


obtain a sandfnce pressure equivalent to a four step calculation, (It 
should once again be noted that for wells less than 5000 feet deep, @ one 
step calculation can be used without significant loss of accuracy, ) 


The following example illustrates the use of this method, 


Example Be? (A Deep Sweet Gas Well) 


A gas well produces 0.75 gravity gas of a known composition 
from an average formation depth of 10,000 feet through 2.441 I.D. tubing. 
The wellhead pressure at a stabilized flow rate at 4.915 MMcfd is 2000 
psia. The flowing wellhead temperature is 110°F and the formation tempera- 
ture is 2h5°F, Calculate the flowing sandface pressure by the Cullender 
and Smith method using two step calculation and adjustment with Simpson's 


rule. 


Solution 


Rea {3 667 psia 


(From analysis) 
Pes 4O8°R 


37.5 G “ = 140,625 


From equation (B-23) 


F od Leh hee ] 0.05158 
(2 .4h1)"" 


F° = 0.002660 


tT. = tip = 1.397 


W 
. 2000 
r 667 


Ww 


ie 0.705 (Figure A-2) 
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Sore 0 


7 Bs ~ = 638°R 
a Gc 
638 
= “ = O h 
Los 1.56 
a 
0 
1 = Tos = 1.728 


p 2000 
aac meals 
(1z), (5T70)(0. 705) gt 


From equation (B-27) 


Step 1 (The upper half of the flow string) 


First trial 


Assume 1 # I 
a Ww 


Solving equation (-29) for Pe 
140,625 = (P. - 2000)(181.4 + 181.4) 


ie 2388 psia 


Second trial 


2388 
Pr oeT 3.580 
G = 0.800 


(t\ - <8 = u.675 
(638)(0.U00) 


(Figure A-2) 


= JOT 


From equation (bB-27) 


I ey ae ] 190.6 
(0.002660) + 4:02)" 


1000 
Solving equation (B-29) 


10,625 = (P_ = 2000)(190.6 + 181.4) 


= 2378 psia 
Third trial 
Pe s 2378 . 3.565 
a 667 
yA = 0,800 (Figure A-2) 


From equation (B-27) 


ui = ie adopt kh) 650 & shanusr 191.2 


2 
(0.002660) + Stee). 


Solving equation (B-29) 
140,625 = (P_ - 2000)(191.2 + 181.4) 


Pa = 2377 peia 


Step 2 e lower half of the ow strin 
First trial 
Assume a ~ oa 


Solving equation (B-30) for P 


140,625 = ie ~- 2377)(191.2 + 191.2) 


P, = 2745 peia 
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Second trial 


IThS 
> = oe = 4vll 
he o6T ? 


= 0,869 


fr) = raaetrecugay 7 +462 
Tz : 7105)(0. 9 


From equation (B-27) 


WS 


| c —__+__.. - 197.1 


8 
(0.002660) + (4.uB1)" 


1000 


Solving equation (B-30) 
ON eeey Ge Bs ere 2377)(197.1 + 191.2) 
P = 2739 psia 


Third trial 


ac) 
fl 


4 ate = 4,106 


0, 869 


p On 
iS = (705 0. 9 4.471 


From equation (B-27) 


~ 
| 


1 ed pee uh) Nad ee ee 197.4 


Solving equation (B-30) 
1h0,625 = (P_ - 2377)(197.4 + 191.2) 


P| = 2139 pala 


(Figure A-2) 


(Figure A-2) 
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Simpson's rule ig now applied through use of equation (B-28) 
to improve the accuracy of the calculated pressure. 


P -P 


281,250 = = [161.4 + & (191.2) + 197.4) 





Ping Paty 737.8 psia 
Fo) mec 0U0 737.8 = 2738 psia 


The method may also be used for calculating the flowing sandface 
pressure for a well which produces sour gas by applying corrections to the 
compressibility factors prior to their use in equation (B-27). 

A series of calculations confirm the theoretical soundness of 
the Cullender and Smith method. To really determine which method ia best 
for practical application, would require extensive field measurement of 
both top hole and bottom hole pressures. Such measurements as are pre- 
sently available to the Board lend support to the Cullender and Smith 
method, but adequate information of this nature is not available. Since 
the method is applicable under most conditions and because of its theo- 
retical aptness, it is adopted as a standard for converting top hole 


pressures to bottom hole conditions. 


Annular Flow 


In most cases where the well ia flowing in the annulus between 
the casing and the tubing, it is possible to measure the corresponding 
shut-in tubing pressure. The flowing sandface pressure may then be calcu- 
lated by the previously discussed method (Cullender and Smith) for deter- 
mining the static bottom hole pressure. 

However, it is on occasion necessary to calculate the flowing 
sandface pressure of an annular column from the flowing wellhead pressure, 


Rigorous equations for the determination of the friction loss term for 


PM AOL = 


flow in an annulus are not available. It is necessary to use the equations 
for flow in a circular pipe, incorporating an effective diameter for the 
annular space. This is commonly done through the use of the hydraulic 


radius formula which may be written as 


D si cross-sectional area of flow (Be31) 
eff wetted perimeter 


D =o) SE*h (B32) 


where 


D = inside diameter of casing, ft. 
D s outside diameter of tubing, ft. 


Substituting the effective diameter in the friction loss tern, 


equation (B-11) may be written as 


z 
4F os apy (Bella) 
& 2 uf 


Equation (B-13) for the tubing gas velocity may be written 
for the annular gas velocity 18 
0.4152 TZQ 
va 5 5 (B-13a) 
P (D, — oO, ) 


It should be noted that the Reynolds number is determined 
from 
20,011 G 
N = 
Re apres (B-18a) 
as compared to equation (B-18). 


Substituting equations (B-lla) and (B-1l3a) into the basic 


flow equation gives 


© th 


c 


ZEEE LLY Laer e | 21d (y(t a \" asordod(B-rhad 
G i GeO ey Jak ) 
rquation (Bb-l4a) muy be integrated by any of the methods pre~ 
viously discussed, and results in equations which are then used to calcu-~ 
late flowing bottom hole pressures for the case of annular flow, 
Equations (B-23) and (B-24) as written for annular flow, and 
with the diameter in inches, are as follows, 
a (B-230) 
(a, - d,) ‘ Canta, } 
for effective diameters less than },277 inches, and 


el ee oe (B-2ha) 


1.5 2 
(a, - d,) (4, * a, ) 


for effective diameters greater than 4.277 inches. 


B, TWO PHASES (GAS AND LIQUID) PRESENT 
IN THE WELL BORE 


The preceeding section has dealt with the calculation of 
bottom hole pressures from wellhead measurements, for the case where 
only one phase (gas) is present in the well bore. With the trend to- 
wards the discovery of deep, high pressure, retrograde-type, condensate 
reservoirs, the case where both a gas and a liquid phase exists in the 
well bore becomes increasingly important. Also, where wells are pro- 
ducing significant amounts of free water from the formation, consideration 
should be given to the effect of this liquid phase on the calculation 
of sandface pressures. 

In the static case, if two phases exist in the well bore, the 
liquid phase will exist as a "liquid leg" at the bottom of the flow 


string. In this situation the bottom hole pressure must be measured 


ize = 


mR 


directly with a bottom hole preeeure Gauge, or calculated from a know- 
ledge of the level of the liquid in the well bore and the gradient with- 
{n the liquid phase, In the latter case, the pressure is calculated 
from the wellhead down to the interface by the normal means for 4 gas 
well, and then the bottom hole pressure is calculated from a knowledge 
of the liquid gradient, Since the level of the liquid and its gradient 
are seldom known, the normal means of determining the static bottom hole 
pressure in the two phase case, is by direct measurement. 

For the two phase flowing case, in addition to the possibility 
of direct measurement, there are a number of methods available for calcu- 
lating the flowing sandface pressure from wellhead measurements. A 
brief discussion of some of these methods is included in this section. 

The first thing that must be decided when testing a well ina 
gas-condensate type reservoir is whether or not two phases exist in the 
well bore, This is accomplished by determining the phase condition of the 
well bore fluid through s calculation of the hydrocarbon dewpoint temp- 
erature at the flowing wellhead pressure. The dewpoint is the condition 
at which the first droplet of liquid forms, and if the dewpoint temperature 
is found to be above the flowing wellhead temperature this indicates thet 
liquid exists and two phase flow is taking place in the well bore. If two 
phases are found to exist, equilibrium flash calculations can be made to 
determine the relative amounts and the composition of each phase at any point 
in the well bore, 

To calculate the wellheed dewpoint temperature and composition 
of phases in the well bore, use is made of vapor-liquid equilibrium ratios. 


The vapor- liquid equilibrium ratio for a component in the fluid is defined 


by the equation 


ah 7 37- 


: x (B~ 33) 


where 


Yabual mole fraction of component i 

in the vapor phase, 

x = mole fraction of component i 

in the liquid phase. 

The value of Ky for each component depends on pressure, tempera- 
ture and the composition of the mixture, 

Many correlations have been presented for the determination of Ky 
values (48) (74) (81). All of the recent correlations incorporate some con~ 
position variable to account for the distribution of components in the 
mixture. 

The method of calculating dewpoint temperatures and the com- 
position of phases are presented in many references (49). A convenient 
consolidation of the method, along with illustrative examples have been 
presented by Van Wielingen (92), for calculating the dewpoint temperature 
of a natural gas, the relative amounts of the gas and liquid phases that 
will exist when a mixture of gaseous hydrocarbons is subjected to a 
specific temperature and pressure, and the composition of each of the 
phases at these conditions, 

A study of production from gas-condensate reservoirs in the 
Province of Alberta has been made in an attempt to correlate the existence 
of two-phase flow conditions to the liquid content of the fluid, The 
atudy has resulted in the conclusion that if the products of a high 
pressure separator (operating under average conditions) are such that the 
liquid - gas ratio is less than about 30 barrels of liquid hydrocarbons 
per million cubic feet of gas, there is little likelihood that two-phase 


flow occurs in the well bore. On the other hand, if the ratio is greater 


it) - 


than about 60 barrels per million cubic feet, it 4a very likely that 
two-phase flow will exist. It must be stressed that the foregoing are 
merely generalizations and that many exceptions to the rule may be en- 
countered in practice. For this reason, this rule {ae not intended to 
replace calculations to determine the phase condition. 

When the vertical flow of gas-liquid mixtures does take place, it 
differs materially from one phase flow, in two important respects. 

oh The flow pattern is much more complex, and depends on the 
diameter of the flow tube and the properties and relative amounts of the 
phases. 

fae The composition of the total fluid in the flow tube differs 
from that produced because of the difference in densities of the two 
phases. The less dense phase tends to slip past the heavier phase and 
thus the heavier phase accumulates in the flow tube. This phenomenon is 
generally referred to 4s "slip" or "holdup", and in the case of a well 
producing liquid is represented by the gas slipping past the liquid with 
a resulting higher liquid-gas ratio in the flow tube. 

Many authors have described the flow patterns for the upward 
flow of gas-liquid mixtures using various terminology. Gosline (32) 
described the patterns 48 changing from dispersed gas through gas piston 
to dispersed Jiquid, as increasing quantities of gas are added to a 
flowing liquid column, Martinelli et al (56) described the same range of 
patterns 4s changing from liquid viscous - @a@8 viscous through liquid 
viscous - gas turbulent to liquid turbulent - gas turbulent. Govier, 
Radford and Dunn (33) described the two phase flow patternt based upon 
visual examination and correlated their description to pressure drop 
observations. For convenience, the patterns 48 described by Govier et 


al are adopted and used in this discussion. 


SEAS 


As gas ia introduced into a column of liquid flowing at a con- 
stant rate, the gas disperses as bubbles which increase in number and 
size. This is the bubble flow pattern. With an increase in ¢nas flow, 
the bubbles coalesce and form bullet shaped slugs characteristic of the 
Blug flow pattern. As the gas-liquid ratio is increased further, the 
slugs merge and appear to degenerate into froth. A further increase in 
the gas-liquid ratio results in less relative motion between phases, and 
the liquid usually forms a4 ripply layer moving upwarde on the tube wall. 
This is termed ripple flow. A further increase in the gas flow rate re- 
sults in a smoothing of the waves or ripples on the flow wall with gas 
flowing at very high rates through the inner core. This pattern is 
termed film flow. Eventually, at an even higher gas flow rate, the en- 
tire liquid layer will be captured from the wall, and carried in a dis- 
persed state in the gas. This is known as the mist flow region. 

It should be noted that not all of these patterns can be ob- 
served under all flow conditions. In particular, where the liquid flow 
rate is high, the increase in gas-liquid ratio required to move from 
region to region is considerably less, and the transition from pattern to 
pattern is not so distinguishable. 

The complexity of the flow pattern and the existence of holdup 
greatly alter the relationship between pressure drop and pertinent vari- 
ables for two phase flow when compared to single phase flow. To illus- 
trate, pressure losses in two phase flow do not always increase with an 
increase in the flow rate or with a decrease in the flow tube size. 
Occasionally, increased gas flow rates without significant increase in 
liquid flow rates, will serve to reduce the holdup and the resulting 
portion of the gradient due to the weight of the column without greatly 
increasing the friction loss. This results in an overall reduction in 


pressure drop with the increasing flow rate. Continued increnses in 
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the gas flow rate, eventually result in increasing friction losses which 
more than ofiset the decreasing weight of the column. The pressure loss- 
production rate relationship then reverses and is similar to the single 


phase relationship. 


This phenomenon is illustrated in Figure B-3, taken from the 


work of Govier, Radford and Dunn, 


UNIT PRESSURE DROP, FT OF WATER PER FT 
HOLDUP RATIO 





RESSURE DROP REGIMES 0 
O 100 200 300 400 500 
AIR-WATER VOLUME RATIO, CU. FT. PER CU. FT. 





FIGURE 8-3 — TYPICAL TWO-PHASE PRESSURE DROP CURVE. 
(Redrawn from Govier et a/(33) 
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Many attempts have been nmiade to correlate the flow patterns 
in a two phase system, and to evaluate the pressure drop in each region 
(11)(33)(34)(37)(75). All of these investigations have been based on 
laboratory studies at relatively low pressures. In addition, a number 
of attempts (6)(38)(69) have been made to determine the two phase pressure 
drop based on actual field conditions. Unfortunately, these investigations 
have been concerned primarily with the case of gas-lifting in oil wells 
(under conditions of bubble and slug flow), and have not been extended in- 
to flow regions of particular interest in the producing of gas-condensate 
wells. Use of the laboratory based correlations indicate that flow patterns 
for a typical gas-condensate well, in general are in the range of ripple, 
film or mist, but the different correlations are not in full agreement. 

On the basis of a general review of these wegidis. it has been 
concluded that the Ros (75) approach is probably the most accurate method 
now available for application to gas-condensate field systems, both with 
respect to the correlation of flow patterns and to the determination of 
the pressure drop. The reader is referred to the original paper for de- 


tails of the method. 


Conclusions 

There are several laboratory based methods available for calcu- 
lating flowing bottom hole pressures when two phases exist in the well 
bore. These are not consistent amongst themselves, and none have been 
extensively tested for actual field gas-condensate systems. For this 
reason, and because of the likelihood that for gas-condensate wells, 
two phase flow exists for only a portion of the flow string, it is recom- 
mended that whenever two phase flow is suspected the flowing bottom hole 
pressures be directly measured during initial testing, until accurate pro- 


cedures or correlations are developed for calculating sub-surface pressures, 
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VALUES OF ge FOR VARIOUS TUBING AND CASING SIZES 
(From Cullender and Smith (24)) 


Berge 


TABLE B-1 


(Use only for internal diameters less than TO ie 


Nominal 
Size, in « 


1 
Ved / 
1-1/2 
2 
2-1/2 
3 
3-1/2 
\ 


4o1/2 
3/4 


(Use only for internal diameters greater than 


4a 3/1 


5-3/16 


5-5/8 


6-1/4 


OM Dismeeleng 


1 325 
1,660 
1.990 
O.ats 
2,075 
3.500 
4,000 
4.500 
4,750 
4,750 
5.000 
5.000 


5.000 
5.000 
52900 
5.500 
PP ew) 
5.500 
5.500 
5.500 
6,000 
6,000 
6,000 
6,000 
6,000 
6.625 
6.625 
6.625 
2 63625 
6.625 


0.10 
ee 722 G32 


d 


Lb. Per 
Rte 


1.80 
2.40 
apy fs) 
4.70 
6.50 
9.30 
11.00 
12,70 


LG. 25. 


18.00 
18.00 
21.00 


0,10 
r are? 2 


13.00 
15.00 
14,00 
15,00 
17.00 
20,00 
23.00 
25.00 
15.00 
17.00 
20.00 
23.00 
26.00 
20.00 
22.00 
24.00 
26,00 
28.00 


1,049 
1.380 
1.610 
1.995 
2441 
2.992 
3.476 
3.958 
4.082 
4,000 
4.276 
4.154 


Pee a 


4494 
4.408 
5,012 
4.976 
4,892 
4,778 
4,670 
4,580 
552k 
5.450 
5. dba 
5.2h0 
5.140 
6,049 
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APPENDIX C 


GAS FLOW IN THE RESERVOIR 


Gas flow through porous media in a reservoir may be steady- 
state or unsteady-state. Steady-state flow requires that conditions 
(flow rate and pressure) do not change with time while for unsteady~ 
state flow, these conditions do change. For example, unsteady-state 
producing operations will result in a decline in both the flow rate 
and pressure, with time, 

Katz, et al (49) have shown that even though pressure deple- 
tion of a gas pool is normal in an unsteady=-state process, flow may 
be treated as steady-state in regions near the well bore under 
certain conditions, or after the well has flowed at stable conditions 
for lengthy periods of time, Under these conditions, the rate of decline 
in flow rate and pressure will be constant over the reservoir and will 
be very slow, reflecting the ratio of the withdrawals in a set time 
period to the total volume of the reservoir. This type of flow is often 
referred to as pseudo steady-state and is usually treated as steady-state 


for calculation purposes. 


A. STEADY-STATE FLOW IN THE RESERVOIR 

Steady-state flow in the reservoir may be laminar, turbulent, 
or a mixture of both, This was initially demonstrated in 1933 by Fancher, 
Lewis and Barnes (28), and has since been confirmed by many others (8)(9) 
(12)(36). Fancher, Lewis and Barnes made extensive measurements on the 
pressure drop = flow rate relationships for the flow of gases through 
cores of various permeable materials, They expressed their results in 
terms of a modified Fanning friction factor and a modified Reynolds 


number (both incorporating the grain diameter) and defined respectively 


as 


Se LOG 


@D (AP) 
iO 
BEST 2g (C-1) 
a 
and 
DVop 2 
7 ee eee (C-2) 
ReFLB , 
u 
where 
AP = pressure drop, pounds per square foot. 


: =  32,17h (lb-mass_) ft, 


(lb. force) sec. 


L = length of core, ft. 
re) = fluid density, lb-mass per cubic foot. 
us = absolute viscosity, pounds-mass per 


ft. second, 


Ss = diameter of average grain, ft. 
4 rate of flow \ 
Ais 5 CP epaLeny velocity (rere area | 


ft. per second. 

A comparison of equation (C-1) to the Fanning equation (B-11) 
and of equation (C-2) to the Reynolds number equation (B-17), reveals 
that these equations differ only slightly. These differences result 
from modifications made by Fancher et al to accommodate the flow of 
fluids through permeable media rather than through pipes. Their data 
indicates that the flow of fluids through porous media resembles that 
through pipe in that there is a change from laminar to turbulent flow 
as velocities increase, The correlation of tOLB versus NV eFLB results 
in a straight line relationship with a Slope of minus one, on log co- 
ordinates, for conditions of laminar flow, regardless of the porous 


material through which the flow is taking place. At velocities that 


correspond to a NReFLB of about 1.0, the relationship between t OLB and 


= E916 


NpeFLB begins to flatten out as the transition to turbulent flow takes 


place, 


Laminar Steady-State Flow 


The work of Fancher et al showed that for conditions of laminar 


flow 


fie 
t Ms ¥LB 


FLB Notes (c=3) 





where CoLB is a constant characteristic of the porous medium. 
Substituting equation (C-2) into (C-3), and the resultant equation 


into equation (C-1) and solving for ub 


g& D ‘ LP 
1 geee NOSE, (c-h) 
€ ibe oh o~ 


An independent approach to the laminar flow of fluids through 
porous media is based on the early work of Darcy which has resulted in 


the somcalled Darcy Law, commonly expressed as 


ee oie err (6-5) 
where 
ue = apparent velocity under mean flowing 
conditions, cm, per second. 
q = volume flow rate, cu. cm. per second. | 
A = cross-sectional area of flow, sq. cm. 
u = viscosity of the fluid, cp. 
- = = pressure drop per unit length, atm. 


per cm, 


K = permeability, darcys. 


= LOO 


This is the relationship ordinarily employed in petroleum engineer- 
ing work. The permeability 'K' as defined by equation (c-5) is in fact equal 


to the value of gD, /2 for any particular porous medium. 


FLB 
If the differential form of equation (C-5) is combined with the 
Gas Law, equation (A-1), applied to a radial system, and if customary field 


units are introduced 


— EOS x 107° kh P dP 
ae uTzZ(dr/r) (c-6) 


Assuming constant values of k, T, and h, and that the change in 
pressure is not great so uy may be treated as a constant, equation (C-6) may 


be partly integrated and written as 


= LenLoe ox 107° kh (C-7) 


Me in Tp 2 


2 
5 


Equation (C-7) may be used in this form with tables ot [E a (49) 
but where the change in pressure is not great, Z may also be treated as a con- 
Stant determined at the average pressure condition. In this case equation 
(C+7) may be integrated and expressed as 


2 es 


-6 
Orfuo4 x 10” kh foe ae (c~8) 


Q 4 ih ints Can CCE a ae is 
Weal in re rs 
where 


Q a gas flow rate, millions of cubic feet 


per day at 14,65 psia and 60°F, 


k = permeability, millidarcys. 
h = pay, ft. 
P, = formation pressure (at exterior drain- 


age boundary), psia, 


Ee = flowing sandface pressure, psia. 


a DOL & 


ee = average viscosity, cp. 

Ze = average compressibility factor. 
ne = average flowing temperature, °R. 
a exterior boundary radius, ft. 

r = effective well bore radius, ft. 


Equation (C=8), while sound theoretically is limited in its use 
even under true laminar conditions because of the inevitable variations 
in pay and permeability and because of the difficulties associated in 
determining To and ry 

The following example illustrates the use of the steady-state 


laminar flow equation, 


Example C-1 
4,915 MMcfd of 0.6 gravity gas is being injected into a well. 


The formation is 20 feet thick, the permeability is 50 md., and the tempera- 
ture of the flowing gas is 640°R. If the injection pressure is 2300 psia 
at the sandface of a four inch well bore, calculate the pressure 40 feet 


out from the centre of the well bore, 


Solution 


P 672 psia 


: (Figure A-1) 


ly 358°R 


Cc 


Assume laminar steady-state flow, and use equation (C-8) anda 


trial and error procedure, 


First trial 


Assume Pe = Py = 2300 psia 


r 


i 
%, 2300 psia 


im 92" « 


4.2060 B 


re ors 


640 
r 350 


= 0,887 (Figure A-2) 


S| 
t 
" 
kK 
sc 
=<] 
\O 


i) 
! 


0.0125 (1.35) = 0.0169 (Figure A-10 and A-11) 


= 
" 


Substituting into equation (C-8) : : 
(0.7054 x 107°) (20) (50) (P, - 2300" ) 


4e915 = = —T9~O169)(0.007) (O40) 1n(40/0.33) 
P. = 5,290,000 - 320,100 
Pe = 2230 psia 


Second trial 


Assume P, = 2230 psia 


f 
a 2265 psia 
bee 265 wi 
= 73 32.31 
ae 0. 889 (Figure A-2) 
u = 0.0125 (1.33) = 0.0166 (Figure A-10 and A-11) 


Substituting into equation (C-8) 


(0.7054 x 107°) (20) (50) (P,” = 2300°) 
LG LS ee renner reer netanee ae emementenipenee 
(0.0166) (0.889)(640) 1n(40/0.33) 
Pi = 5,290,000 ~ 315,800 
Ph = 2230 psia 


Hence the answer converges to 2230 psia in two trials. 
A similar approach may be used to determine the permeability- 
psy product where the sandface pressure and the pressure at some distance 


out in the reservoir are known at a particular flow rate. 


= BUS h= 


Turbulent Steady-State Flow 


During the production of gas from a reservoir, the linear velocity 
and the corresponding Reynolds number of the gas increases as it approaches 
the well bore. Laminar conditions will prevail throughout the reservoir 
if the linear velocity of gas at the sandface is such that the corresponding 
modified (l"ancher, Lewis and Barnes) peeraiis number is less than approxi- 
mately one. For higher flow rates, turbulent conditions will exist at the 
sandface, and will prevail some distance back into the reservoir, This 
distance will depend upon the velocity of the gas. Thus, the common situ- 
ation arises where turbulent conditions prevail for some distance out from 
the well bore, and laminar conditions exist from this intermediate point to 
the boundaries of the reservoir, 

Forscheimer (30) has proposed that flow under both laminar and 
turbulent conditions may be handled by a quadratic equation including two 
components; the pressure gradients due to the viscous shearing forces and 
to the inertial forces, The Forscheimer equation (in absolute units as 
used in equation (C-5)), may be expressed as 


V 


AP ua 2 
= Se. |) + i) 
Tiwa) poatabee (C=9) 
where 
B = a constant of proportionality referred 


to as the turbulence factor, ea 
The meaning of 8 can be better illustrated if we consider the 
case of a well producing from the centre of a radial reservoir of a size 
such that turbulent flow exists near the well bore, The total pressure 
drop from the outer boundary to the well bore may be expressed as 


AP, = AP, + AP ap (G30) 
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where 


AP: = total pressure drop. 


th 
AP, = pressure drop in region of laminar flow. 
AP op = pressure drop in regions of transition- 


al and turbulent flow. 
If we further define: 
Lo = length of regions of laminar flow. 
L = length of region of transitional and 
turbulent flow. 
Then a friction factor for the turbulent-transitional region may be 


expressed by a relationship similar to equation (C+1) 





EDGE 
ta See (C#11) 
A = 
ae : Lop 
then by combining equations (C-4), (C-11) and (C-10), 
2G f 
PDBee fuk 2 
AP. =o A 
T BN eh a Vag Slap (c-12) 
gD gD 
eg cg 


this equation is identical to the Forscheimer equation (C-9) where 


2 Cop bly ii 
2 7 a 
Be AL 


(C=13) 


& D AL (C=14) 


The turbulence factor is analogous to a friction factor divided 
by a grain diameter, 
If the differential form of equation (C-9) is combined with 


the Gas Law, equation (A-1) and is applied to a radial system, expressed in 





B1195 - 


customary field units, and integrated assuming constant values of k, T, 


wy, Hh, Q and Z, then it may be expressed as 


6 rs : 
pep 2 . 2418 x 10 Hbat a ole ge 4478 oar? (2. 2) 
aaa a gag \r 


12 
8 kh r h re 
(C=15) 
where 
Pe = formation pressure, psia. 
i = flowing sandface pressure, psia. 


k = permeability, md. 


h = pay, ft. 

ee average viscosity, cp. 

T, = average temperature, rR 

Ze Be average compressibility factor. 
re = exterior boundary radius, ft. 
r, = effective well bore radius, ft. 
8 = turbulence factor. 


Q = gas flow rate, millions of cubic 

feet per day at 14,65 psia and 60°F. 

Before this equation may be used, the turbulence factor 'B' 
must be evaluated, Janicek and Katz (45) have made measurements of the 
turbulence factor for cores of various permeabilities and porosities and 
have correlated their data in the form of a plot of turbulence factor 
versus permeability with porosity as a parameter. Cornell and Katz (17) 
have reviewed these data and proposed a simpler relationship of the 
turbulence factor to the permeability itself, They found that for dolo-= 
mites, limestones and sandstones an approximate relationship exists be- 


tween the turbulence factor and the permeability. This relationship may 
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be expressed as 


hake eddex NOxe (c#16) 
k 3 
where 
k = permeability, millidarcys. 


Cornell and Katz have correlated pressure drop data on consoli- 
dated and unconsolidated cores, in terms of a modified friction factor, 


f..., and a modified Reynolds number, N These factors are really very 


CK 
similar to those employed by Fancher et al, and are defined as 


ReCK* 





64 Bo Ap 
fox =F V2 PL (C~17) 
and 
N e BK ae 
ReCK 10 (C=18) 


6r33ex, 10 u 
Combining equations (C-1) and (C-17), it is apparent that the friction 
factor defined by Fancher et al and by Cornell and Katz are related 


through 





and combining equations (C=2) and (C-18), it is clear that the Boas in the 
Fancher relationship is replaced by 8K/6.33 x 1079, in the Cornell-Katz 
modified Reynolds number, 

If the value of '8' given by equation (c-16) is incorporated 
into equation (C-18), for a radial system, and converted to field units, 


the modified Reynolds number at the sandface becomes 


* = 2209101 GQ 


ReCK | rh ti x i/3 (c-19)) 
where 


P = pressure, psia, 


L 5 length of flow, ft. 
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G me specific gravity of gas. 

Q = flow rate, millions of cubic feet 
per day at 14.65 psia and 60°F. 

k = permeability, millidarcys. 

u = fluid viscosity, centepoises, 

h = pay, ft. 


~~ effective well bore radius, ft. 


The Cornell and Katz friction factor plot (17) (49) indicates that 
the change from laminar to turbulent flow for all types of porous media 


begins at an N of about 0.1. The transition region where the flow is 


ReCK 


partly turbulent corresponds to the N range of from 0.1 to about 1.0. 


ReCK 
The Turbulence Contribution Factor 
Equation (C=15) for laminar-turbulent radial flow may be expressed 


in the form of a quadratic equation 





P me P aes aQ + vQ° (C20) 
it 8 
where 6 
: 1,418 x 10 ui Tin (r,/¥,) 

BS) kh 
and ; 

«od sdahi xohlows $5 .0hentl 96t. +9 

AG ie r 


8 5a 
Govier (35) has shown that if a turbulence contribution factor 
'r ' tg defined as 
L 
b 
= -_ Cx21 
36 an autre ( ) 


then the laminar-turbulent flow equation may be written as 


2 
Pe sa Fie =aQ Ei (Cu22 ) 


and the turbulence contribution factor By, may be evaluated by combining 


ale 


equation (C-15) and the emperical relationship between g and kK, equation 


(C#16), 
If this is done and the small value of 1/r, is neglected 
te Me ct OeOUO NG 1 


ia 1/3 
r hk Me in(r,/r,) 


(G25) 
It may be further shown that for the common case where the 


drainage radius is of the order of one-half mile and the 1n(r,/r,) is of, 


the order of 9, equation (C-23) reduces to 


G 
=. C=24 
Fo AL + 0,01 es ( ) 


r 
lig 
or 


Pe Mei et ne} al (C-25) 


t Nreck? 

The solution of equation (Cc-2h) ig presented in graphical form 
as Figure C-1, 

At a Reynolds number of Ool, where Cornell and Katz found that 
the change from laminar to turbulent flow begins, equation (C=25) shows 
that the turbulence contribution factor would be approximately 1.01, 

To solve the laminar-turbulent equation (C-15) or its equivalent 
(C+22), the turbulence factor can be evaluated from equation (C~23) or 
for certain drainage conditions from Figure C-1, It should be emphasized 
that neither the equation nor the figure are precise, since they incorpor- 
ate the emperical correlation between the turbulence factor 'B' and the 
reservoir permeability 'k' as developed by Cornell and Katz, 

The following example illustrates the use of the turbulence con- 


tribution chart, Figure C~l, in the calculation of sandface pressures, 


Example Cx2 
200 MMcfd of Ua ls gravity gas is being produced from a well come 


pleted in a formation with a 20 foot pay section and a permeability of 
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60 md, The shut-in reservoir pressure is 2500 psia, the well bore radius 
is 4 inches, and the average reservoir temperature is 660°R. A calcue 
lation has revealed that the average compressibility factor is 0.80 and 
the viscosity is 0.02 ep. 

Calculate the sandface pressure assuming that the steady-state 


radius of drainage is one-half mile. 


G ¥ O ° e) 


= . = 28,4 
Wari 0.02)(0.33) (20 


From Figure C-l 


Fy = 1,073 


Substituting into equation (C-22) 


6 
2 2  (1.073)(1.418) (10°) (0.02)(0.80) (660 8.98 
a BLS oe 20) (60 
= 2500° = 601,800 
P = 2377 psia 


8 
Direct use of the laminar flow equation for this problem gives 


a sandface pressure of 2385 psia. 


B, UNSTEADY=-STATE FLOW IN THE RESERVOIR 

Flow in a porous medium is defined as unsteady-state when the 
mass flow rate into an element of the medium is not the same as the mass 
-flow rate out of that element. Hence the fluid content of the element 
changes with time, Such changes in the amount of fluid in the porous 
medium are possible because of the compressibility of fluids. The vari- 
ables in unsteady-state flow, in addition to those already used for 
steady-state flow, are time, the porosity of the formation and the com- 


pressibility of the fluid. 


The basic approach to the solution of an unsteady-state problem, 
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must be consideration of the system as a separate entity. A material 
balance equation known as the continuity equation, is then written to 
account for the fluid injected into and produced out of the system. When 
this equation is combined with the Gas Law (the equation of state) and 
Darcy's Law (the equation of motion), it provides a method of solution of 
unsteady-state problems for known boundary conditions, 

For the radial flow case, the equation of continuity may be 


written as follows 


adel), ov 4 @ (c-26) 
where 
V = velocity in the radial direction. 
ig = radius of the volume being con- 
sidered. 
) = porosity of the volume being con- 
sidered, 
rs) = fluid density, 
t = time, 
Equation (C=5), Darcy's Law for radial flow may be written 
as 
ip) cee, as 
uooor (C-27) 


and the Gas Law, equation (A-1) may be written as 


MP 
ep = ORT (C-28) 


Combining equations (C-26), (C-27) and (C-28) and the 
assumption that viscosity, compressibility factor and permeability 


are constant, the following equation may be written for laminar, 
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radial flow in a horizontal plane. 





we a 
: pe re drow t, ues oe 
ae r or KP at (C-29) 


It should be pointed out that equation (C-29) is not the 
rigorous differential equation, but is only approximate in that the 


right-hand side should be ud 9P. This approximation is not serious 
Ke aoe 


unless drawdowns are very large. 

Equation (C-29) is non-linear (the dependent variable pe appears 
to powers other than one) and insoluble as written. If the 'P'’ in wo/KP 
is assumed to be constant at a value of the average pressure, Pa then 
the equation is linear and can be solved by several methods, These methods 
include the analytical solutions for constant production rate by Hurst and 
Van Everdingen (42), Cornell (20)(21) has applied the Hurst and Van 
Everdingen solution to constant rate gas flow problems. Cornell and Katz 
(17) have adapted the Schmidt graphical method used for unsteady-state 
heat flow problems to the solution of the linearized equation. This method 
is particularly useful in visualizing transient pressure gradients. 

One of the most convenient approaches to solving unsteady~state 
flow problems is through definition of certain dimensionless quantities 
and solution of the linear equation in normalized form, 


If the following dimensionless quantities are defined: 








2 
Ig 
dimensionless pressure Po =|— (c-30) 
Es D Pe 
dimensionless radius ry = — (Ca31) 
8 
= 
2,634 x 10 KEG 
dimensionless time har 3 a deep) 
D thee 2 
as 
Vet Lor 10°u ZT Q 
aaa (C=33) 
dimensionless flow rate Qn = 5 , 


hk Pe 


oy Pe - pe et 
dimensionless drawdown Py = a <3 + peal} 
D Py 
a 
where 

i = pressure, psia. 

be od radius, ft. 

k = permeability, md. 

t = time, hours. 

Lee: average viscosity, cpo 


Z average compressibility factor. 


> = gas filled porosity, fraction. 


[ = average temperature, °R. 


Q = flow rate, millions of cubic feet 
per day at 14.65 psia and 60°F. 
h = pay, ft. 
Pir,t) = pressure at any radius, r, after 
well has been produced at constant 
rate, dps for time, t. 
then equation (C-29) may be normalized (made dimensionless), by dividing 


the numerators by oe and the denominators by es and expressed as 


ee (c-35) 


9 9 9 
2H) 7 eD ty 


and generalized solutions in terms of the defined ratios may be developed 
for various boundary conditions, 
Hurst and Van Everdingen (44) have solved this equation for a 


constant flow rate 'Q' and for an infinite reservoir, and have shown that 


for ty greater than 100, 


Pls eeaig oe 0. 80907 ) (C-36) 


They also solved the equation for a constant flow rate from a finite 


See 


reservoir with a closed circular outer boundary, and have shown that 


0.5147 Ot 3r M - ur (1n ry) = ae - 1 (c-37) 








for values of ty that correspond to times when the boundary effects are 
felt. Hurst and Van Everdingen also solved the equation for the case of 
a finite reservoir with a constant pressure at the outer boundary, but the 
solution is complex and for this reason is not reproduced here in equation 
form, 

These and other solutions may be expressed in tabular or graphi- 
cal form, the latter being perhaps the more convenient for general use. 
Katz et al (49) have presented such graphical representations of the 
solutions for certain boundary conditions, Aziz and Flock (4) have re- 
cently published graphical representations of the solutions of the un- 
steady-state flow equation for the constant production rate case. Their 
graphical presentation is especially convenient because it combines the 
solution for the boundary conditions of greatest interest in a single 
plot. The form of plot is produced as Figure C-2, and includes solutions for 
the following boundary ere 

1. Infinite radial reservoir. 

2a Finite reservoir, sealed outer boundary. 

30 Finite reservoir, constant pressure at the outer boundary. 


The following examples illustrate the use of this figure. 


Example C=3 


A gas well located in a very large reservoir which has been 
shut in for a lengthy period, is produced for five hours at the rate of 


3.01 MMcfd, Calculate the sandface pressure at this time if the following 


reservoir data is given. 


Pe = 2300 psia 
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T = 640°R, 

¢ = 0,15 

k ® 28 md. 

h BLS Tt 

Vet 80 2 att 


8 
The sandface pressure must be assumed to calculate Ze and Wee How= 
ever, for simplification this portion of the calculation will be left out 
of the trial and error procedure. It will be assumed that the following 
values are known, 


es 0,89 


rans 0.0168 


Solution (For five hours producing from a large reservoir, this can be 
considered as an infinite reservoir, ) 
First trial 


Assume P. = 2100 psia 


1 


Re 2200 psia 


From equation (C-32) 


; -h 
t = 20034 x 107" (28)(2200)(5) _ 128, 800 


y (0,0168)(0.15)(0.5)* 
80 from Figure C-2 (infinite case) 
Bas 6,28 
In this case P, may also be calculated from equation (C-36), 


From equation (C=33) 


6 
Del (0. )(0.89) (640) ( ) 
a 8_x_10_ (0.0168) (0.89) (640)(3,01 = 0.0230 


(12) (28) (2300)? 
Substituting into equation (C=34) 


1 2300° . p 2 
Cd on eT 


0.0230 2300° 


ime] 
ul 


2127 psia 


= Pieter 


Second trial 


Assume ie = 2127 psia 


P ae - oe 2214 psia 


a 


ml 
tit (2,634)(10 28)(2214 = 129,600 
(0.0168) (0.15) (0.25) 
From Figure C=2 


4 = 6,29 


6 
— (1.428) (107) (0.0168) (0.89) (640)(3.01) _ 0.0230 


q 
2 (12) (28) (2300)* 
Once again solving equation (C=34) 


P, = 2127 psia 


Hence the sandface pressure converges to 2127 psia in two trials. 


Example C=4 
Assume that the well in Example C-3 is located in a radial 


reservoir only 2000 feet in diameter, and is produced at 3.01 MMcfd 
for 100 hours. Reservoir data are: 


Py = 2300 psia 


ah = 640°R. 
Q = 0.15 
kK = 28 md, 


h = le ft. 


x = 0.5 ft. 
8 
a? a. 0.89 
u = 0.0168 Cpe 


= 206 - 


Solution (This is the case of a finite reservoir with a sealed outer 


boundary, ) 
ry = = Be 2000 
r On 


To avoid numerous trials, we will assume a sandface pressure we 
know to be approximately correct, 


Assume i = 2075 psia 


P = 2300 _+ 2075 = 2187 psia 


a 2 
From equation (C-32) 
2.634 x 1074 28) (2187) (100 
ty = <eSde X30 A28)(2187) (1.00) = 2,560,000 
(0.0168)(0.15)(0.5) 
60 from Figure C-2 (finite case-sealed outer boundary ) 
P, = 8,1 
From equation (C=33) 


6 
a = 1.418 x 10 1020166)(9.69) (640) (3.02) = 0.0230 
(12) (28) (2300) 
Substituting into equation (C=34) 
i 2300° . p 2 
8.1 = ———e Ss 
0.0230 2300 


P 
8 


2075 psia 

So the assumed sandface pressure is correct and no further trial 
is necessary, 

The problem involving the finite reservoir-constant pressure 
at the outer boundary is handled in the same manner ag Example Cm4, except 


that the appropriate portion of Figure C-2 igs used to determine P This 


1 
approach may also be used to calculate the permeability-pay product when 
the sandface pressure and a pressure at some distance out in the reservoir 


are known for a particular flow rate, It may also be used to estimate the 


aetsey fees 


time to reach a certain pressure, or the degree of stabilization. The 
latter matter is dealt with at length in Appendix D. 

In the case where the rate of flow from a reservoir has been 
changed, multiple transients result in the reservoir. This is known as 
the flow after flow case, and the pressure distribution at any time may 
be calculated, The approach involves determining the pressure change 
associated with each change of flow rate, and summing these pressure 
changes for the time over which each change was in effect. The calcue 


lation may be made through use of the following equation. 


e * BS Let 
=11g.) Pee ( f Jape Bot )P + + 
2 D q q (9 qa cee. 
P fete Deer SUN et. D, “GeD3'= t., 
(qa. -@ a 
D, Din) th (c-38) 
where 


t = total elapsed time since the first 


- 
flow rate began. 
t, = total elapsed time since the second 
flow rate began. 
th = total elapsed time since nth flow 
rate began, 
dn = first dimensionless flow rate, 
1 
Tree nth dimensionless flow rate, 
n 
ae = dimensionless pressure change 
number evaluated at time the 
P. = dimensionless pressure change 
n 


number evaluated at time th 


The following example illustrates the application of equation 


(c~+38). 
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Example C=5 


A gas well which has been shut in for a lengthy period, is opened 
and produced at a rate of 3.01 MMcfd for 4,0 hours. The rate is increased 
to 6.02 MMcfd and maintained at this rate for another 4,0 hours, Calewiath 
the sandface pressure at the end of the second production period, assuming 
an infinite reservoir. 


Reservoir data are: 


P = 2300 psia 
iar = 0.0168 Cpo 
h Seleercs 
Ze = 0.89 
= oO 
T,. 640°R. 
k = 28 md. 
ry SAG oe Lin 
> = 0.15 


Solution 
First trial 


Assume P = 2000 psia 


Pa = 2150 psia 
we (2654) 40") (28) (250) WP Ae 
Dy 0.0168) (0.15)(0.25 , 


ct 
i 


ay 
2,634)(10 28)(2150)(4 
D, Sy ee 100 ,700 


From Figure C=-2 


e = 6,49 


I = 6,14 
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6 
af (1.418)(107) (0.89) (0.0268) (640)(3.02) _ 0.0230 


(12) (28) (2300) 
maaan (1.418) (10°) (0, 89) (0.0168) (640) (6.02) ite PON 
Ds (12)(28)(2300)° 


Substituting into equation (C-38) 


2300° - p > 
a = — (00230) (6.49) + (0.0230) (6.14) ] 
2300° 
P = 1937 psia 


Second trial 


Assume Py = 1937 psia 


op eR oa ey 


a 


2 
an 
tt (2.634) (107) (28)(2119)(8) 198,500 
T (0.0168) (0.15) (0.25) 
=i 
ty = (2,634) (107°) (28)(2119)(4) 99 ,230 
2 (0.0168) (0.15) (0.25) 
From Figure C-2 
P = 6,49 
ey 
P = 6,14 
to 
6 
a 1.418)(10- )(0.0168)(0.8 {ens 3-01) | 0.0230 
i (12) (28) (2300) 
and 
q = 0.0460 
Ds 


Substituting into equation (C-38) 


2300° - p* 
—— [(0,.0230)(6.49) + (0.0230)(6.14)] 
2300 













P, = 1937 psia 


Hence the answer converges to 1937 psia in two trials. 7 : Me Li 
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APPENDIX D 


STABILIZATION AND RELATED MATTERS 
LL LT AE a SE Rt i a a ED 


The simplicity of form of the equations describing the steady 
state radial flow of gas in a reservoir as compared to their unsteady- 
state counterpart, has long suggested that well test data obtained under 
steady-state conditions may be more easily interpreted than similar 
unsteady-state data, However, it is clear that the production of a well 
in a finite reservoir is inherently an unsteady~state process, and that 
Steady-state analysis of the situation is an approximation, Nonetheless, 
consideration of the unsteady-state relationship, particularly as the 
graphical solution of Cornell and Katz, indicates that at any instant 
of time the steady-state equation is satisfied if movement of the gas is 
assumed to originate at the apparent radius of drainage - i.e. at the 
radius obtained by extrapolating the linear portion of the P =a. ea 


£ 


versus ln (r/r,) curve, to Beas This means that data points taken at 
times corresponding with equal apparent radii of drainage would simulate 
steady-state data, This approach prompts one to think of a condition 

of apparent stabilization for which the steady-state equations are ade- 
quately descriptive, Clearly, this condition must be close to that which 
would be reached after an appropriately long time, and must be changing 
very slowly with time, The time required to reach such a condition de- 
pends upon the arbitrary definition of the condition, A logical 
definition is that when the transient pressure-radius profile first 
reaches the edge of the spacing unit in which the well is drilled, For 

a fully developed field, this is the point where the profile would inter- 


sect that of the neighboring wells. For a field not completely developed, 


while the profile would eventually extend further, its restriction to the 


mobos 


edge of the spacing unit simulates full development. Janicek and Katz (45) 
have shown that the ratio of the @pparent radius of drainage (as here de- 
fined) to the effective radius (to the edge of the spacing unit), is 
approximately 0.60, Thus, the time to attain this arbitrarily defined 
stabilization is that time required for the apparent radius of drainage 
to move sixetenths of the distance towards the edge of the pattern, 
This time may be calculated from unsteady=-state theory. 

The equation for Steady-state laminar radial flow of @ gas may 


be expressed as 





Hf (P.- ~- cee) 
qi * RT er Re 
D n(x, /r_) pee (D=1) 
where 
qq = the dimensionless flow rate, 


PS = formation pressure, psia, 
P_ = flowing sandface pressure, psia, 
r «= effective well bore radius, ft. 
r, = the apparent radius of drainage, 
ft. (which is equal to 0.60 re 
at the time the transient pressure- 
radius profile reaches the edge of 
the drainage pattern, ) 
If the equation is combined with the unsteady-state flow equation 


in the form of equation (c~34) 


1 Pie - PF r,t 
P,o=— : (C=3 ) 
Ip re 


the following relationship which ig valid for the condition of apparent 


Stabilization results 


P - P 
5 = in Or, aa) al (D-2) 


t 


and if the unsteady-state flow equation is written in terms of the sand- 
face, equation (D-2) may be written as 


P, = ln eure) (D=3) 


The dimensionless time necessary to attain this condition can 
be determined from Figure C-2, for the appropriate boundary conditions. 
The actual time corresponding with this time, may be determined through 
equation (D-4) 

h 


2.634 x 10 kP it. (Deh) 


+ So ere eh oe a ae AE, 
D - 2 
u, or, 


where 
k = permeability, md, 


P = average pressure, psia. 
t = the time of stabilization, hours, 
7 = average viscosity, cp. 


r = gas filled porosity, fraction. 


Py = effective well bore radius, ft. 


This determination requires a knowledge of the average pressure, 
which in turn requires that the sandface pressure is known. The sandface 
pressure aed may be determined from a knowledge of the apparent radius 
of drainage, which coincides with the arbitrarily selected definition of 


stabilization, and the steady-state equation (C-8). 


0.7054 x 10 + kh (E. - Po ) (c-8) 





Q = 
Py In(r,/r,) 
where 


Q = gas flow rate, millions of cubic 
feet per day at 14.65 psia and 60°F, 


mele = 


k = permeability, md. 

h = pay, ft. 

2 = formation pressure, psia. 

Pp «= flowing sandface pressure, psia. 

u = average viscosity, Cpe 

% average compressibility factor. 

[ = average temperature, °R. 

r = apparent radius of drainage, ft. 

r = effective well bore radius, ft. 

The above outlined approach may be simplified by development 
of .an equation which enables one to estimate directly the time of stabi- 
lization. In the case of a finite reservoir with a constant pressure at 
the outer boundary, Figure C-2 (solution of the unsteady-state flow 
equation) shows that the dimensionless pressure drawdown, Pro becomes 

.constant at some tos for each dimensionless radius, Tye This particular 
portion of the figure was taken from work of Hurst et al (44), and the 
relationship between 'r.‘ and 't," at the instant when "P,' becomes con= 


D D 


stant may be expressed as 


0.48858 


Gavan 2.6408 t) (D-5) 
which may be approximated as 

r. = "constant" t ori 6 

D D (D=6) 


where the constant varies depending on the value of tps If equation 


(D-6) is combined which equation (D-4), the following results 
2 


r 
= " " ni a 
ee constant cP (Dat) 


a 


where the constant varies from some 500 to about 900 over the range of 


dimensionless times normally encountered, 


ss aa 


For the infinite reservoir case, Janicek and Katz (45) have 


expressed the relationship between 'ry! and "ty! at the time ne becomes 


constant as 


0.5 (p-8) 


ui 
Lae) 
ct 


"D D 


This expression results in the following equation 


9D 
(a 
3 0s 


= 9g a9 (D=9) 
ies 


t 


Muskat (63) has developed a method for calculating the time 
of stabilization for a radial water-drive oil reservoir, When applied 
to a finite radial gas reservoir with a fixed outer boundary, the approach 


results in the following relationship 


2 
1g 
t. = 950 SUEUR RD (D-10) 
k P 3 
& 


Another criterion of stabilization that is often used is the 
time when the sandface pressure is dropping at a specified low rate per 
unit of time, 

For a known constant production rate, the pressure drop after 
any time may be measured from the slope of the appropriate ,! versus 
"t.' curve, in Figure C-2, From equations (C-34) and (C+32) this rate 


D 


of pressure decline may also be expressed as 








dP 
nh lied 3th (D=11) 
dt, : seibee x lok Pit 
~ = 
wa >? 8 
If we combine equations (D-11) and (C~33) and assume that Me Ze P.? 


Q, k, T,, and >, are constants, and that OBS P. (this is a reasonable 


approximation where the drawdown is not large the following equation 


results, 
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dP, 5.356 x 10> orb (- me | (D=12) 
dt, Zo ,9 dt 


where 


_—& = the pressure decline at the sand- 
face, psi per hour. 
In the most commonly occurring case, that of the finite reservoir with 


a sealed outer boundary, Van Everdingen and Hurst (44) have shown that 


\ 4 2 
: : Oreste en, Epo Ary (1n rp) = 2 Sy or 2 (c-37) 
ttle pkere % 2 2 
To) 
ke = ey 
2 (D=13) 
dt, ln - 1 
2 
or since r => aL 
; D 
dP, 2 ry 
—t az 2 a (D214) 
dat 58) r 
D D ne 


Combining equations (D-12) and (D=14) 


2 -3 2 
() : 5.356 x 10 or, h (- —.) 
re ZT ,2 at (D-15) 


Equation (De15) shows that when the pressure transient reaches the 
edge of the drainage area, the rate of pressure decline will depend upon the 
flow rate, and the reservoir and gas properties, If stabilization is defined 
as a specific rate of pressure decline, the "stabilized radius of drainage” 
calculated from the equation will vary, not only from well to well, but for 
a particular well as the flow rate changes. A different (dP/dt) would be 


required for each flow rate, to ensure equivalent stabilization. Since this 


- P22). - 


situation is undesirable, the rate of pressure decline method for defining 
stabilization is rejected, and an approach more closely related to a set 
radius of drainage is recommended, 

Considering equations (D-7) based on work of Hurst, (D-9) based 
on work of Janicek and Katz, (D-10) based on work of Muskat, and several 
equations by others as summarized by Van Poolen (91), it is concluded that 
the following equation should give a reasonable estimate of the time to 
stabilization, it should be noted that the constant is rounded to 1000, 
and that the reservoir pressure is used in place of the average pressure 


in equations (D-7), (D-9) and (D-10), 


rs pad 
. k P 
where 
t. = time of stabilization, hours. 
2 ee radius to the exterior boundary, ft. 
yu, = average viscosity, cp. 
b = gas filled porosity, fraction. 


k = permeability, md, 

P = formation pressure, psia. 

The following examples illustrate the use of this equation 
for calculating the time to stabilization, and also show that the rate 


of pressure decline method is not a useful approach for estimating 


times to stabilization, 


Example D=1 

Calculate the time of stabilization for a low permeability, 
low pressure well, given the following reservoir data, 

Ee 625 psia. 

T. = 520°R. 
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k . = 10 md, 
hres BeQue? 


h = 20 ft. 


ie malls Sh 
ie = 341°R. 
eC 
Pa = 665 psia. 


r Mi ie duiat tod 


r. = 2640 ft. (1 section spacing) 


Da 
Ze = 0,90 (at an estimated average 
reservoir pressure of 550 psia.) 
eae 0.011 cp. (at an estimated average 


reservoir pressure of 550 psia.) 


The well is flowed at a constant flow rate of 500 Mcfd. 


Solution 
Time of Stabilization 


Substituting into equation (D-16) 


a (2909) (2640) (0,023){0220) 
8 10)( 625 


te = 2450 hours 


Rate of Pressure Decline 


From equation (D-14) 


a, = 2 BGho = lLe (Pex LO 
Combining this with equation (D-12), we can determine the rate 


at which the sandface pressure will be declining when the flow conditions 


are semi-stabilized, 
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iotymuose neeee e408? 40.20} (025") (20) es 
; 0.90)(520)(0.5 dt 


at 0.0031 psi/hour 


So the sandface pressure will be declining at a very low rate 
when stabilization occurs and such a rate would have to be specified to 


result in stabilization equivalent to that attained after 2450 hours. 


Example D-2 
Calculate the time of stabilization for a high pressure, high 


permeability well, given the following reservoir data. 
P = 2700 psia, 
i = 640°R, 


k = 1000 md, 


b = 0.09 
h = 120 ft. 
G = 0,81 
if = 416°R. 


P = 672 psia. 


Re = 0,25 ft. 
sa shoal 0 8 ae ap 
oer (at an estimated average 


reservoir pressure of 2690 psia.) 
u = 0,019 cp. (at an estimated average 
reservoir pressure of 2690 psia.) 
The well is flowed at a constant rate of 30.0 MMcfd, 
Solution 


Time of Stabilization 


Substituting into equation (D-16) 
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3 2 
om (1000) (2640 rte Sys Hours 


Rate of Pressure Decline 


ise 19 


From equation (D-14) 


a 


dP 2 
we. a 2(S3) = 1.79 x 10 Y 


Combining with equation (D-12) 


: ap 
-8 5,356 x 10 3(9,09)(0.25°)(120 mas 
meelOvng 5 0.79 0)(30 dt 
dP, 
rem te 0,075 psi/hour 


It can be seen that if adequate stabilization is to reflect a 


set drainage radius, no one pressure decline rate can be used as a 


criterion. 


Conclusions 


These examples illustrate that the time to achieve this arbi- 


trarily defined stabilization (to the edge of a one section spacing unit), 


as calculated from a theoretical approach may be very long. In fact, 


from Example D-2, one may conclude that if the defined criterion of stabi- 


lization must be attained on each flow point, a conventional back pressure 


test is virtually impractical to run in all but the most permeable re- 


servoirs, 


This calculated time may be adjusted on the basis of experience, 


with the object of determining the time after which the change in the flow 


rate and pressure are inconsequential with respect to the effect on flow 


in the reservoir, It is recognized that attainment of this condition is 


not always practical and an arbitrarily selected "maximum time" to stabi- 


lization is desirable. This Board has chosen a time of fifteen days. This 





means that where the calculated time of stabilization is greater than 


fifteen days, it is considered that by the fifteenth day the changes 

occurring as the radius of drainage continues to move outwards in the ree 
servoir, will have little effect on the interpretation of a flow test, 
bearing in mind the uses to which the test results will be put. This does 

not necessarily mean that a flow point need be extended for the lesser of 

the calculated time to stabilization or fifteen days. However, where the time 
to stabilization is longer than the length of the extended flow period, the 
stabilized conditions should be calculated using the method outlined in 


Section 5 or any other acceptable approach, 
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APPENDIX E 


THE THEORETICAL VALIDITY OF THE ISOCHRONAL 


METHODS OF TESTING GAS WELLS 


The starting point in a discussion of the theoretical validity of 
the isochronal test, is the steady-state equation (C=-15) for combined 


laminar-turbulent flow. 


_ 
oe ™ Eas 818 10° Ha2,7,8 (=) i 
3.14 x ita 2 it 1 
2 BG 2,7, fer - = (E=1) 
where 
Pe = formation pressure, psia. 
fs = flowing sandface pressure, psia. 
k = permeability, md. 
h = pay, ft. 
Cry en average viscosity, cp. 
T, =  sverage temperature, °R. 


Z bad average compressibility factor. 


a 
aes apparent radius of drainage, ft. 
Peles effective well bore radius, ft. 
B = turbulence factor, 

Q = gas flow rate, millions of cubic 


feet per day at 14,65 psia and 60°F. 
Equation (E-1) may be approximated by the following equation 
r 
GQ ie ~ Pa = Q (=) (E-2) 


r 
8 


where 'n' depends upon the total contribution of the two terms on the 


right-hand side of equation (E-1) being unity when the second term makes 
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no contribution, and for example being 0.95 when the maximum contribution 
is some 20 percent. The ‘'n’ in equation (E-2) is represented as being a 
constant, This is an approximation, as a perusal of equation (E-1) clearly 
shows it to be a variable, In fact, 'n’ is a constant (equal to unity) 
only when there is zero turbulence, 

For a series of flow rates, if each of the flow rates is extended 
to a condition of stabilization, the effective radius of drainage is the 
edge of the drainage area, and the apparent radius is a constant. Equation 


(E-2) may then be written as 


= Q (B=3) 
where 

C = the stabilized back pressure 

coefficient, 

Several authors (44)(45)(55) have shown that the effective 
radius of drainage for a particular reservoir is a function only of time. 
Since a definite relationship has been shown to exist between the effective 
radius of drainage and the apparent radius of drainage, it follows that for 
any series of flow rates carried out for a constant duration following a 
prolonged period of shut-in (an isochronal test), the apparent radius of 
drainage will be a constant and equation (E-2) may be expressed as 
ab 


eee 
ae, are = Q (E-4) 


where C'' is some coefficient other than the stabilized coefficient. 

The turbulence term in equation (E-1) is essentially independent 
of the radius of drainage, rae and is the same for a short time isochronal 
flow period as it would be for a longer stabilized flow period, However, 
the laminar term (and thus C'’ as designated in equation (E-h)) differs. 


As a result, the 'n’ in equation (E-4) wii differ from that of equation 


FLOW RATE 





MPSS. = 


(E-3) even though it corrects for the same amount of turbulence, This 
difference in the value of the exponent 'n' will depend on the degree of 
stabilization which exists at the time interval sexectaa for conducting 
the isochronal test. Several calculations have indicated that the change 
in 'n' will be small, even in the extreme case where the drainage radius 
affected during the test is very small compared to the stabilized radius 
of drainage. Experience verifies this, and in fact indicates that in 
most instances, changes in the exponent 'n' as the flow period is extended 
in length, are immeasurable. 

The results of an isochronal test may then be plotted as 
(P - Ba) versus Q and will result in a straight line with a slope such 
that the exponent ‘'n' is essentially the same as would be experienced 
under conditions of stabilized flow. If this isochronal line is positioned 
to reflect stabilized conditions, the results of a fully stabilized con- 


ventional back pressure test are duplicated, 
Modified Isochronal Test 
From unsteady-state flow theory it can be shown as follows that 


the results of a modified isochronal test are almost identical to those 


obtained from an isochronal test. 





AES ain 


Isochronal Test 


For the isochronal test the time between each flow period is 
not equal to T but is sufficient to permit the pressure to return to Pre 
Consider the final flow rate, Qy,° This is plotted against 


Pe - Bh where Po is the sandface pressure at the end of flow rate 


M (ay) 
From equation (C=3)) 
i F =- P g 
2 D 9 
ee Toe, 


P, corresponds to the time since the 7th flow rate (a, > Q),) began. 
T 


T 
From equation (C=-36) iat 1/2 (1n pee 0.809) 
ff T 
Assuming the effects of changes in pressure and viscosity can be neglected 
then P. = 1/2 (1n "constant" T + 0.809), 
T 
Combining the constants into "a" and "b" 
P = 1n bT + a 
t 
i 
2 2 
and Pr = P 
ig ‘i 
3 = Qn (1n bT + a) (1) 
Pe T 


Modified Isochronal 


In this case the time between each flow period is equal to T and 


Poe Py and P¢ are not equal to Pro 


Again consider the final flow rate, Q),° This is plotted against 


q Pp + (q - q ) P + (q. e ) P + ( ae ) P 
2] D qa q a. 


‘ ue . I? ate + (4p, = aay ie 


MAST 


=q = 0 
i 6 


where I * qn 
and P. = time since first flow rate began =6T 
Pe = time since second flow rate began = 5T 
PL = time since third flow rate began = 417 
bse = time since fourth flow rate began = 3 T 


Ee = time since fifth flow rate began = 2T 


E = time since sixth flow rate began = T 





mR Te SekGe iP) todo oP. as eee (PD pep. «| 
P. D, ty t, D, t, t), De t. te 


> Cama be, a Sn ea [2p » (6) + a= anv (or) +a] =a, 1n/6 
eat Hy to oh ap 5 


Similarly q. (P, =P, )=q, In/¥\andq. (P, =P, ) =q nf 2 
D,  t; t), D, (3 Ds * ts te Ds 


< = q, 1n/6\+ q 1n/4\+ q 1n/2 (2) 
Pu ey @, “ayo, He) 


Also from equation (C-38) 





2 2 
P - P 
my isi q id | qa ed, qa ip’ F gi 
P Dy ti) D, D, t, dD D, t3 D), D3 ty, 
(qewumagin Et Leatyuioatlg bs) tbo aft Tp 
oe D), ae De D. te Da to 


where q =q =q = 0 
D D D 
2 4 6 


and P now corresponds to time since first flow rate began or to 


ot: ed ioc Nene ogee edt ret site te | leah peti Weed <P ete a 
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2 2 
then P wolf. 
i ee o (P, = P, ) + P 
=q (P, -P, ) +a, (P pobre tl qQ 
re eat to eels ty Deane Ye “6 hte 
pio hae 
nmel =, In(q\ a 1n(5\+ ap 1n(3) ay (In bT + a) (3) 
Pay Vie his 3 5 2 "i | 
Subtracting equation (2) from (3) 
ee ee 
uml eg (In -in6) +a) (4n 3 - in 4) +q) (In g- in 2) 
Pe 1 5 5 3 i 3 5 2 1 
sak (1n bT + a) 
i; 


= in/35\ +a. 1n/15\ +a, 1n/3\+aq. (in dT + a) (4) 
a en 


Comparing equation (1), the isochronal test, and equation (4), the modified 


isochronal test, the difference in the fourth point, flow rate Q), iss 


py , (33) ap, (ae %p, *(2) | 


Similarly, it can be shown that the third point will differ by 


- 1n/15 \+ i1n/ 3 
? fey ot) 


and the second point will differ by 


q, in/3. 
 |%, y 
The first point will be identical, 
¥ 
The validity of the modified isochronal test may also be il- 
lustrated by a numerical example as follows, showing that the apparent 


radius of drainage will be nearly the same for a series of flow rates of 


equal duration, even though each flow rate does not begin at the same 


pressure, 


Example E-1 
A new well drilled in a reservoir with a virgin pressure of 1000 psia 


(from which no production has taken place), is flowed at a constant rate of 


a 


1.0 MMcfd for a period of three hours, The well is then closed in and flowed 
for alternate three hour periods, The second flow rate is at 1.5 MMcfd, 
the third is at 2,0 MMcfd, and the fourth flow rate ig at 3.0 MMcfd. Cale 
culate the apparent radius of drainage at the end of each of the flow periods, 
The following well data is given, 

a = 550°R, 

k = 30 md. 

) = 0.20 


h = 20 ft. 


G = 0.6 
x (o) 
te 360°R, 
Pe = 672 psia. 


r i (GG) saec 
Zz = 0,86 


U = 0.013 cp. 


Solution 
The flowing sandface pressures would normally be measured, but in 
this example they will be calculated from the unsteady-state flow equation, 


The example calculation assumes laminar flow. 


(a) First Flow Rate 
First trial 


Assume P = 920 psia 


p: = 2000 +920. 06, cate 


a 2 
From equation (C32), the dimensionless time which corresponds 


to three hours is calculated. 


alk 
t= £22634 _x 10) (30) (960) (3) = 1.40 x 10? 


a (0.013) (0.2) (0.25°) 
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For a short flow period such as three hours, the reservoir is 
assumed to act as an infinite reservoir, 80 from Figure C=2 


P, = 6,32 


The dimensionless flow rate which corresponds to 1.0 MMcfd is 


calculated from equation (C-33). 


_ (1.418 x 10 ®)(0 2013) (06 ae (120) . 9. 01h5 


(20) (30) (1000° ) 


Substituting into equation (C-34) 


. 1 (1000° = Poe 


6.3 2 





0.0145 1000° 


PS = 953 psia 


Second trial 


Assume Ps 953 psia 


Pa = cat = O77 psia 


From equation (C32) 


-) 
c= ABsG3h x 20 CSO OTINS) 1.43 x 10° 


4 (0.013) (0.2)(0.25°) 
From Figure C=2 


a ae 6.33 


1 0.0145 


Substituting into equation (C-34) 
1 (1000% = P_*) 

63° 5 = 
0.0145 1000 





P = 953 psia 
So the sandface pressure after three hours will be 953 psia. 


From the steady-state flow equation (C-8), the apparent radius 


of drainage may then be calculated, 


« €3) - 


Low 4027054 x 107°)(20)(30)(1000 ~ 953°) 


(0.013) (0.86)(550) In(r,/0.25) 


ya 138s Ltr 


(b) Second Flow Rate 
The pressure at the end of the second flow rate can be calculated 


from equation (C-38), for the multiple transient case. 


Pol “ Pas 
—— ns P + (q - q ) P + (q = 2 jor 
Pas ea Sy. Ds nabs ute mica ag ice 
2.634 x 1074 30 ) 2 
ty 62.630 x 10 (30) (977) (9) = 4.58 x 26 
1 (0.013)(0.2)(0.25°) 


From Figure Ce2 


P = 6,87 


1.418 x 10° 0.013)(0.86 0)(1.0 


= 0,0145 
1 (20) (30) (1000°) 


t. = 2,86 x 10° 
P = 6,68 


I = 0 
1.43 x 10? 


ct 
a 


P = 6,33 


0.0218 





= = (0,0145)(6.87) = (0.0145)(6.68) + (0.0218) (6.33) 


P = 927 psia 


The sandface pressure at the end of the second flow period is 
927 psia. The radius of drainage at this time can be calculated from 


equation (C-8), if the pressure is known at the start of the flow period. 
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This pressure (to which the well will be built up during the three hour 
shut-in), is also determined from the multiple transient flow equation. 


10007 =» Pp 


——— = (0.0145) (6.68) - (0.0145) (6.33) 
1000 


e = 997 psia. 


So from equation (C-8) 


0.7054 x 107° 20) (30 2 _ 92 2 
0 


1.5 = ““=9°013)(0-60)(550) intr, /0.25 


r = 121 ft. 
a 


(c) Third Flow Rate 
The pressure at the end of the second shut-in period is then 


calculated from equation (C-38). 


ty (for 12 hours) = 5.72 x 10? 
1 
Be ieee Oe 
a 
aq, = 0 
D,, 
1000° « oh 
—————— = (0,0145)(7.02) = (0.0145)(6.87) + (0.0218)(6.68) = 
1000 


(0,0218) (6.33) 
% = 995 psia 
The pressure after the third flow period is also calculated from 


equation (C38). 


ty (for 15 hours) = 7.15 x 10? 
uf 
P = Tart 
‘oh 
an = 0.0290 
p) 
1000° = 7 
—— a (0.0145)(7.14) = (0.0145)(7.02) + (0.0218)(6.87) = 
1000 
(0,.0218)(6.68) + (0.0290) (6.33) 
P = 900 psia 


= 23{ = 


from equation (C-8) 


-6 2 2 
2,0 = 427054 x 107") (20)(30)(995° - 900°) 


(0.013) (0.86)(550) 1n(r./0.25) 
re = 123 ft. 
(d) Fourth Flow Rate 
The pressure at the end of the third shut-in period is calculated 
from equation (C38), 
=D, (for 18 hours) = 8,58 x 10° 


P= 7,23 
Si) 


ee (0.0245)(7.23) = (0.0145)(7.1h) + (0,0218)(7,02)-— 


(0.0218)(6.87) + (0.0290)(6.68) = (0.0290) (6.33) 


E. = 993 psia 


The pressure after the third flow period is also calculated from 


equation (C-38). 


ty (for 21 hours) = 1.00 x 10° 
1 
P = 7.31 
oni 
qa, = 0.0435 
Dy 
1000° = p @ 
= (0,0145)(7.31) ~ (0.0145)(7.23) + (0.0218)(7.14) + 
1000° 


(0.0218)(7.02) + (0.0290)(6.87) = (0.0290)(6.68) + 


(0.0435) (6.33) 


P B4E psia 


From equation (C-8) 


waao = 


3,0 = 4021054 x 107°) (20) (30)(993° ~ 846°) 


(0,013) (0,86) (550) 1n(r/0.25) 


se We Lose foe 


It can be seen for this particular reservoir that the apparent 
radius of drainage after three hours at any of the four flow rates is 
nearly the same, even though the formation pressure prior to each rate is 
not the same. This helps to justify the contention that the modified 


isochronal test yields results close to those of a true isochronal test. 
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